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Abstract 

The "principle of the fermionic projector" provides a new mathematical frame- 
work for the formulation of physical theories and is a promising approach for 
physics beyond the standard model. The book begins with a brief review of 
relativity, relativistic quantum mechanics and classical gauge theories, with the 
emphasis on the basic physical concepts and the mathematical foundations. The 
external field problem and Klein's paradox are discussed and then resolved by 
introducing the so-called fermionic projector, a global object in space-time which 
generalizes the notion of the Dirac sea. The mathematical core of the book is 
to give a precise definition of the fermionic projector and to employ methods of 
hyperbolic differential equations for its detailed analysis. The fermionic projector 
makes it possible to formulate a new type of variational principles in space-time. 
The mathematical tools for the analysis of the corresponding Euler-Lagrange 
equations are developed. A particular variational principle is proposed which 
gives rise to an effective interaction showing many similarities to the interactions 
of the standard model. 

The main chapters of the book are easily accessible for beginning graduate 
students in mathematics or physics. Several appendices provide supplementary 
material which will be useful to the experienced researcher. 
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APPENDIX A 



Connection to the Fock Space Formalism 



In this appendix it is shown that for an observer who is making measurements 
only in a subsystem of the whole physical system, the description of a many-fermion 
system with the fermionic projector is equivalent to the fermionic Fock space formalism, 
provided that the number of fcrmions of the whole system (including the particles of 
the sea) is infinite. The following consideration applies in the same way to either a 
space-time continuum or to discrete space-time. Before beginning we point out that 
the action principle, from which the fundamental physical equations can be deduced, 
involves the fermions only via the Dirac action <^', {i(^+B—m)^>. For the formulation 
of the Dirac action one only needs on the fermionic Fock space the time/position 
operators and the operator ^, which are all one-particle operators. Therefore, we can 
say that many-particle operators (like for example in the four-fermion coupling of the 
Fermi model) are not essential for the formulation of the quantum field theory of the 
standard model. Having this in mind, we may here restrict attention to one-particle 
operators^. 

Let P be a fermionic projector acting on the vector space H. The one-particle 
observables correspond to operators O on H. Our subsystem is described by a non- 
degenerate subspace K C H; we decompose H as & direct sum H = K ® L with 
L = K^. We assume that the observables are localized in N; i.e. they are trivial on 
L, 

Oil = 0,L . (A.l) 

We choose a (properly normalized) basis ^'i, . . . , of the subspace P{H) C H, and 
decompose the states *j in the form 

Substituting into (3.2.1), we obtain for the many-particle wave function the expression 

* = E (-i)''^' ( A ) ( A *i ) ' (A-2) 

■Kevin) yeTT j y^TT / 

where V{n) denotes the set of all subsets of {l,...,n}. For measurements in our 
subsystem, we must calculate the expectation value <^\0\^>f^, where the operators 



^Online version: For the description of entanglement, it is indeed necessary to consider two-particle 
observables; see the paper "Entanglement and second quantization in the framework of the fermionic 
projector" (in preparation). 

^We remark for clarity that this expectation value does not coincide with that of a measurement 
in nonrelativistic quantum mechanics. Namely, in the continuum, the scalar product <.|.> involves 
a time integration. But one can get a connection to nonrelativistic measurements by considering 
operators O with a special time dependence (which, for example, act on the wave functions only in a 
short time interval [t, t + At]). 
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A. CONNECTION TO THE FOCK SPACE FORMALISM 



O act on the Fock space according to 

A • • • A *„) = (O^i) A • • • A + ^1 A {0^2) • • • A + • • • + ^1 A • • • A , 

and where <.|.>/.' is the scalar product on tlic Fock space, induced by the scalar 
product <.|.> on if. It is useful to rewrite the expectation value with the statistical 
operator S, i.e. 

<*|C'|*>F = tTipiSO) with S = |*><*|f, 

where trj? denotes the trace in the Fock space. Using (A.l), we can take the partial 
trace over L and obtain, applying (A. 2), 

<^|0|^>F = ti Fj, {S^ O) with (A.3) 

n 

= Yl Yl I ^^eTT *^^>< Aje^' ^f\FK (A.4) 

k=0 7r,7r' e •P(n), 

where tr^;-^ is the trace in the Fock space Fk = e^io ^'^ ^ generated by K. Thus 
our subsystem is described by a statistical operator on Fk, which is composed 
of mixed states consisting of different numbers of particles. Since the constants c^^^t^' 
depend on the wave functions outside our subsystem, we can consider them as 
arbitrary numbers. 

In the limit when the number n of particles of the whole system tends to infinity, 
(A.4) goes over to a statistical operator of the form 

= E E \^'a'><^?\^K (A-5) 

fc=0 a,/3=0 

(k) (k) h 

with arbitrary complex coefficients c^^ and /c-particle states € Fj^. This statistical 
operator differs from a general statistical operator in that it is diagonal on the 
A;-particle subspaces (i.e. that the wave functions in the "bra" and in the "ket" of 
(A. 5) are both A;-particle states); more precisely, has, compared to (A. 5), the 
more general form 

00 00 

^.^n = E E -S' i^i'^^x*?!.^. • (A.6) 

k,l=0 a,/3=0 

We remark for clarity that a pure state of the Fock space ^ G Fk has a decomposition 
* = Ylk^=o ^k^^'^\ and thus the corresponding statistical operator is 

S = \-f><^FK = E ^^^l*^*^^><*^'Vi^ ■ 
k,l=0 

This statistical operator is a special case of (A.6), but it is not of the form (A. 5). 

The difference between (A. 5) and (A.6) becomes irrelevant if we keep in mind 
that all physically relevant observables commute with the particle number operator. 
Namely in this case, every expectation value reduces to the sum of the expectation 
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values in the /c-particle Fock spaces, 

oo oo 

k,l=0 a,f3=0 
oo oo 

k=0 a,/3=0 

If we choose the coefhcicnts c^jj] in (A. 5) to be c|^^ — \ these expectation values 
are also obtained from the statistical operator , 

tr FAS^onO) = tvpAS^'O). 

We conclude that it is no loss of generality to describe the subsystem by the statistical 
operator . 
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Some Formulas of the Light- Cone Expansion 

This appendix is a compilation of some formulas of the light-cone expansion. More 
precisely, we list the phase-free contribution to the light-cone expansion of the Green's 
functions (cf. Def. 2.5.5). According to Def. 2.5.5 and Theorem 2.5.6, the light-cone 
expansion of the Green's functions is immediately obtained by inserting ordered expo- 
nentials into the line integrals. Furthermore, as explained after (2.5.45), the formulas 
can be applied directly to the fermionic projector; they then describe the singularities 
of P{x, y) on the light cone. Without loss of generality, we restrict attention to the left 
handed component of the Green's functions. We compute precisely those contributions 
which will be of relevance in Appendix G and in Chapters 6-8. The following formulas 
were all generated by a computer program, see [F6] for details. 

We begin with the perturbation by a chiral perturbation to first order. The phase- 
free contribution (denoted by a corresponding superscript on the equal sign) is 

XL i-s {xl4r + Xr4l) s)ix, y) P^"=^''' i 0{{y - xf) + 0{{y - xf) 



+XLS^^\x,y)e 


r dz [0, 1 1 0] i^Au) 

' X 


(B.l) 


-XL>S(°)(x,j/) / 

J X 


y 

dz [0, 1 0] 4^ 


(B.2) 


+XL s(°)(x,i/)4^ 


(x) 


(B.3) 


+\XL 5(°)(x,y)^ 


/ r dz [0,0 \0] {Ul) 

J X 


(B.4) 


-XL S^'\x,y)^ ^ 


["dz [1,0|0] (d'Au) 

' X 


(B.5) 


+lxLS^'\x,y)^ 


e r dz [0,0 \l]{n Au) 

J X 


(B.6) 


+XL S^'\x,y) e 


r dz [0,1 1 1] i^Au) 

Jx 


(B.7) 


+XLS^'\x,y) [ 

Jx 


' dz [0,2 \0]{n4L) 


(B.8) 


-2xL-SW(x,y) J 


^""^^[0,0 11] {0Au), 

X 


(B.9) 



where again ^ = {y — x). The notation ^ 0{{y — x)^) means that we leave out all 
contributions which are of the order 0{{y — x)°) and have a leading factor ^. This 
formula has the disadvantage that it contains partial derivatives of the chiral potential; 
it would be better for physical applications to work instead with the Yang-Mills field 
tensor and the Yang-Mills current. Therefore, we introduce left and right handed 
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gauge-covariant derivatives D^/^, 

and define the corresponding field tensor and current as usual by the commutators 

i^fc = I [D'^^, Dl] , jf = Fil] (c = L or R). (B.IO) 

In the case of an Abelian gauge field, this formula reduces to the familiar formulas for 
the electromagnetic field tensor and current, 

F/fc = djA^k - dkAcj , if = dikA^ -nA^i. 

Notice, however, that in the general case of a system of Dirac seas, (B.IO) involves 
quadratic and cubic terms in the potential. 

By substituting (B.IO) into (B.1-B.9) and manipulating the line integrals with 
integrations by parts, one can rewrite the phase- free contribution in a way where the 
linear terms in the potential are gauge invariant. For example, we can combine (B.l- 
B.3) by transforming the line integrals as 

e r dz [0, 1 I 0] i^ALk) = e r dz [0, 1 I 0] (7^iffc + dk4L) + 0{Al) 

J X J X 

= e r dz [0, 1 I 0] 7^F/fe - + dz [0, | 0] 4^ + 0{Al) . (B.ll) 

Jx J X 

This procedure yields (in the non-Abclian case) quadratic and cubic terms in the 
potential which are not gauge invariant. Fortunately, these gauge-dependent terms are 
all compensated by corresponding contributions to the higher order Feynman diagrams. 
We thus obtain 



XlY.{{-s{xl4r + Xr4l))^ s){x,v) P'^"^^-"'' ^0{{y-xf)+0{{y-x) 
fc=o 

+XL5(0)(x,y)f r dz [0,110]^^,^ 

Jx 

+\xL S^^^ {x, y) i £ dz [0, I 0] 7^^= iffe 

-\xLS^'\x,y)iC [d^ [0,0|l]jf 

-iXL {x,y)i r dzi [0, 1 | 1] F^^ f dz^ [0, 1 | 0] 

Jx J Zl 

+XLS^'Hx,y)e r dz [0,1|1] (^jf) 

J X 

+XLS^^\x,y) r dz [0,2 \0]jtl'' 

Jx 

-iXL S^^\x, y) i r dzi [0, 1 I 1] F^^ f dz^ [0, 1 | 0] 

J X J Zl 

+iXL S^^^ (x, y) ee r dzi [0, 3 I 0] 7' F^- dz2 [0, | 1] 

Jx J Zl 
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+iXL S^'^ {x, y) ee r dzi [0, 2 | 1] f dz^ [0, 1 | 0] 7' Ft, 
-2ixL S^^\x, y) r dz, [0, 2 | 1] T dz2 [0, 2 | 0] {^F, 

J X J Zl 

-2ixLS^'Hx,y)^ie f dz^ [0,2|1] {^F^^) T [0, 1 | 0] Fj 

J X J Zl 

+iXL S^^^ {x, y) f dzi [0, 2 | 1] j'^F^j f dz2 [0, 2 | 0] 

Jx J Zl 

— XL S'-'^ {x, y) e 1^ dz, [0, 2 | 0] Fji J' dz2 [0, | 0] ^^^'F^ 
r dz^ [0,2|0]7Viffe T [0, 1 | 0] Vi^^f 

^ J X J Zl 

+2ixL S^'^ {x,y)^i r dz, [0, 3 | 0] 7^iffc dz2 [0, 1 | 0] 

Jx J Zl 

-2ixL 5(1) (x, y) e r dzi [0, l\l]Ftj f"" dz2 [0, 1 | 0] 7fcF|* 

-2xL S^^\x,y)£,iiie 

X r [0, 4 I 0] y f dz2 [0, 1 I 1] F^. T ^2:3 [0, 1 | 0] Ff 

-2xL S^^\x,y)i,i^e 

X r dzi [0, 3 I 1] 7'F,t r dz2 [0, 3 I 0] F^. r dz3 [0, 1 I 0] Ff 

J X J Zl J Zl 

-2xL S^'^{x,y)i^^e 

X r [0, 3 1 1] F^fc r [0, 3 1 0] 7'f,^ r dz^ % 1 1 o] f^ 

X J Zl J Z2 

-2xL s^^\x,y) e^je 

X r dzi [0, 3 I 1] F^fc r dz2 [0, 3 I 0] F^^' r d^a [0, 1 | 0] Vf, 



<ki 
L 



'Z2 

We call this formulation of the phase-free contributions purely in terms of the Yang- 
Mills field tensor and the Yang-Mills current the gauge invariant form of the light-cone 
expansion. 

It remains to consider the scalar /pseudoscalar perturbation; i.e. we must study 
how the dynamic mass matrices Yj^^j^^x) show up in the light-cone expansion. We 
begin with the case of a single mass matrix. To first order in the external potential, 
the corresponding Feynman diagram has the light-cone expansion 

XL m {-s {-xlYr - XrYl) s){x,y) 

phase-free 1 ^ ^(0) ^ ^ ^ | ^^^^^ 

+XL m {x, y) Yl{x) + 0{{y - xf) . (B.12) 

The higher orders in the chiral potentials yield no phase-free contributions. The next 
orders in the mass parameter are treated similarly. The contributions quadratic in m 
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are 

oo 

XLTT? {{-s4LrsYLs{-4^srYns{-4^sr){x,y) 

ni,n2,rt3=0 

phase-free ^_ ^(o) ^^^^^^ ^ ^^^^^^ 

+iXL 5(1) {x, y) r dz [0, 1 | 0] 7^' (D^Yr) 

J X 

+iXL rr? 5(1) (x, y) f dz [0, 1 | 0] t' I^^.Xl) Yr 

J X 

-iXL rr? S^^\x, y) Yl T dz [0, | 0] T' l^j^fi) 

J X 

+XL 5(1) (x, y) e r dzi [0, 2 | 0] j^F^'i f dz2 [0, 0\0]Yl Yr 

J X J Zl 

+XL rr? {x, y) C f dzi [0, 2\0]Yl Yr f dz2 [0, 1 | 0] 7^F^^ 

+^0{{y-xf) + 0{{y-xf), 

whereas there is only one terms cubic in m, 

00 

XLTu' ^ {{-s4^rsYLs{-4j,srYRs{-4^sr^YLs{-4^sr)ix,y) 

ni,n2,n3,n4=0 
phase-free ^3 ^(1) £ q , 

To the order ~ and higher all contributions are on the light cone of the or- 
der ^0{{y - + 0{iy - xf). 

The above Feynman diagrams completely characterize the Green's functions to the 
order 0{{y — x)^) on the light cone. Notice that in agreement with Theorem 2.5.6 we 
get only a finite number of phase-free contributions. 
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Normalization of Chiral Fermions 



In this appendix we describe a method for normahzing chiral fermions. The main 
difficulty is that for a proper normalization one needs to give the chiral fermions a small 
rest mass; this will be discussed in Section C.l for a single Dirac sea in Minkowski space. 
In Section C.2 wc develop a method for analyzing the normalization of chiral fermions 
with a small generalized "mass," whereas Section C.3 gives the general construction 
including the infrared regularization and the interaction. 

C.l. Massive Chiral Fermions — Preparatory Discussion 

Before introducing the infrared regularization, we need to understand how a chiral 
Dirac sea can be normalized in infinite volume using some kind of "(5-normalization." 
To this end we consider a non-interacting left-handed fermionic projector in Minkowski 
space, 

P{x,y) = XLtm{x,y)\^=o , (C.1.1) 
where we set ^ 

tm = -iPm-k,n) 

with pm and km according to {2.2 A, 2.2.5). Naively, products of this fermionic projec- 
tor vanish due to chiral cancellations, 

P^{x,y) = j (fzP{x,z)P{z,y) = J d'^zxLto{x,z) XLto{z,y) 

= I d^z XL XR to{x, z) to{z, y) '°'^=""^ . (C.1.2) 

However, this formal calculation has no meaning in the formalism of causal perturba- 
tion theory §2.2 because in this formalism we are not allowed to multiply Dirac seas of 
the same fixed mass. Instead, we must treat the masses as variable parameters. Thus 
before we can give a mathematical meaning to products of chiral Dirac seas, we must 
extend the definition of a chiral Dirac sea to non-zero rest mass. 

Giving chiral Dirac particles a mass is a delicate issue which often leads to confusion 
and misunderstandings. Therefore, we discuss the situation in the example (C.1.1) in 
detail. In momentum space, the distribution t„j, m > 0, takes the form 

t^{k) = m + m) 6{k^ - m'^) Q{-k^) . 

On the mass shell, the range of the (4 x 4)-matrix -|- m is two-dimensional; this 
corresponds to a twofold degeneracy of the eigenspaces of the Dirac operator (J^ — m) 
for any fixed k. If m = 0, the Dirac equation splits into two separate equations for the 
left- and right-handed component of the spinor, and this makes it possible to project 
out half of the eigenvectors simply by multiplying by XL^ 

P{k) = XL 1^ Sik^) Qi-k^) . (C.1.3) 
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If m > 0, this method cannot be apphed because the left- and right-handed subspaces 
are no longer invariant. In particular, the product Xi^m for m > is not Hermitian 

and is not a solution of the Dirac equation. Nevertheless, we can project out one of 
the degenerate eigenvectors as follows. For given k on the lower mass shell we choose 
a vector q with 

kq = and ^ -1 . (C.1.4) 

A short calculation shows that 

[tmik), Pi] = and {pif = 1 

(where p is again the pseudoscalar matrix (1.2.13)). This means that the matrix p^ 
has eigenvalues ±1, and that the Dirac equation is invariant on the corresponding 
eigenspaces. Projecting for example onto the eigenspace corresponding to the eigen- 
value —1 gives 

P,n{h) := ^ (1 - p4) {H + m) 5{e - n?) e(-feO) . (C.1.5) 

Thus, similar to the procedure in the massless case (C.1.3), is obtained from by 
projecting out half of the Dirac eigenstates on the lower mass shell. But in contrast 
to (C.1.3), the construction of P„ depends on the vector field g, which apart from 
the conditions (C.1.4) can be chosen arbitrarily. A short calculation shows that Pm is 
idempotent in the sense that 

Pm Pm' = ^("l - fn) Pm ■ (C.1.6) 

The distribution (C.1.5) can be regarded as a generalization of the chiral Dirac 
sea (C.1.3) to the massive case. In order to make this connection clearer, we now show 
that (C.1.5) reduces to (C.1.3) in the limit m \ 0. Thus, for fixed k ^0 and variable 

m > 0, we let k be on the lower mass shell, k{m) = (— y^|fcp -|- m^, k), and choose q{m) 
such that (C.1.4) is satisfied. A simple example for q is 

k ^ 



q{m) = -i-\k\,^\k\^ + m^-^\. (C.1.7) 
m \ \k\J 

In this example, k and mq coincide as m \ 0; more precisely, 

k — mq = 0(m^) . 

This relation holds for a large class of functions q{m). Therefore, it seems general 
enough to concentrate on the situation where 

k — mq = m^ V with. v{m) = 0{m^) . (C.1.8) 

Solving this relation for q and substituting into (C.1.5) gives 

Pm{k) = \(l-p^+ mpi] m + m) S{k'^ - m^) e{-k^) . (C.1.9) 

Using that on the mass shell ]f,{]j, -|- m) = m(]f, + m), wc get 

P„(A;) = ^ (1 - p + mpV) (i^ + m) 5(A;2 - m^) e(-A;°) . (C.1.10) 

If now we take the limit m \ 0, we obtain precisely (C.1.3), i.e. 

lim Pm = P (C.1.11) 

m\0 



C.l. MASSIVE CHIRAL FERMIONS - PREPARATORY DISCUSSION 



209 



with convergence in the sense of distributions. This calculation shows that (C.1.5) 
indeed includes (C.l. 3) as a limiting case and that the dependence on q drops out as 
m\ 0. 

The distribution (C.1.5) gives a possible definition of a massive chiral Dirac sea. 
However, it would be too restrictive to use only (C.1.5) as the basis of our construction, 
because there are other common ways to give chiral Dirac particles a rest mass. These 
alternatives are more general than (C.1.5) in that the wave functions are no longer 
solutions of the Dirac equation. To give a simple example, one could describe a massive 
left-handed Dirac sea for m > by 

Pm{k) = (xL^ + j) S^k" - ^) G(-feO) . (C.1.12) 

This distribution has the advantage over (C.1.5) that it is Lorentz invariant, but it is 
clearly not a solution of the Dirac equation. As m \ 0, we again recover the massless 
chiral Dirac sea (C.1.1). We compute the operator product PmPm' in momentum 
space. 




where in the last step we used that m,m' > 0. Note that in the last line the sum- 
mand m/8 appears (instead of the summand m/4 in (C.1.12)), and therefore P„j is 
not idempotent in the sense (C.l. 6). On the other hand, one can argue that (C.1.6) is 
a too strong normalization condition, because we are interested in the situation when 
the masses of the chiral particles are arbitrarily small, and thus it seems sufficient 
that (C.1.6) should hold in the limit m,m' \ 0. In this limit, the problematic sum- 
mands m/A and m/8 both drop out, and thus we can state the idempotence of Pm as 
follows, 

lim {Pm Pm' - S{m - m') Pm) = . (C.1.13) 

m,m'\0 

The above example shows that, in order to have more flexibility to give the chiral 
Dirac particles a mass, it is preferable to work with the weaker normalization condi- 
tion (C.1.13) instead of (C.1.6). Comparing with the naive calculation (C.l. 2), one sees 
that introducing the mass changes the behavior of the operator products completely, 
even if the masses are arbitrarily small. Therefore, we refer to the limit m, m' \ 
in (C.1.13) as the singular mass limit. 

For the correct understanding of the singular mass limit, it is important to observe 
that, in contrast to operator products as considered in (C.1.13), the formalism of 
the continuum limit is well-behaved as m \ 0. Namely, in the continuum limit we 
consider an expansion in powers of m. The different orders in m have a different 
singular behavior on the light cone. In particular, to every order on the light cone only 
a finite number of orders in m contribute. Thus to every order on the light cone, the 
m-dependence is polynomial and therefore smooth. Expressed in terms of the kernel, 
the limit m \ is singular when we form the product P{x, z) P{z, y) and integrate 
over z (as in (C.l. 2)). But if we take the closed chain P{x,y) P(y,x) and consider 
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the singularities on the hght cone, the hmit m \ is regular and well-behaved. This 
justifies why in Chapters 6-8 it was unnecessary to give the neutrinos a mass and 
take the limit m \ afterwards. We could treat the neutrino sector simply as being 
composed of massless chiral particles. In particular, the chiral cancellations in the 
formalism of the continuum limit are consistent with the singular mass limit. 

Our next goal is to develop the mathematical framework for analyzing the singular 
mass limit for a fermionic projector with interaction. Clearly, this framework should 
be general enough to include the examples (C.1.5) and (C.1.12). Thus we first return 
to (C.1.5). After writing in the form (C.1.10), it seems natural to interpret the 
leading factor as a generalization of the chiral asymmetry matrix X. This is indeed 
convenient in the vacuum, because introducing the operator by 

Xm{k) = ^{l-p-mpiik)), (C.1.14) 
we obtain in analogy to the corresponding formulas for massless chiral particles that 

Unfortunately, the operator X^ does not seem to be useful in the case with interaction. 
The reason is that X^ depends on the momentum k, and this leads to the following 
serious difficulties. First, the /^-dependence of X^ makes it very difficult to satisfy the 
analogue of the causality compatibility condition 

X^{i^ + B-m) = {i^ + B-m)Xrn. 

As a consequence, it is in general not possible to commute the chiral asymmetry matrix 
through the operator products of the causal perturbation expansion; in particular 

tm and tm Xm do in general not coincide (where tm is the interacting Dirac sea 
as defined via the causal perturbation expansion). Even if we assume that there is 
a canonical definition of the fermionic projector Pm obtained by suitably inserting 
factors of X^ and X^ into the operator product expansion for f^, we cannot expect 
that the correspondence to the massless Dirac sea is respected: i.e. in the case with 
interaction, (C. 1.11) will in general be violated. In order to explain how this comes 
about, we point out that our argument leading to (C.1.8) was based on the assumption 
that k converges to the mass cone as m \ 0. More precisely, if limm\o k[ra) is not 
on the mass cone, the function v will diverge like v{'m) ~ m~^, and so Xm{k) will not 
converge to X as m \ 0. Thus lim^^o ^m = ^ only if in this limit all the momenta 
are on the mass cone. But in the causal perturbation expansion off-shell momenta also 
appear (note that the Green's functions are non-zero away from the mass cone). This 
means that in the limit m \ 0, the momenta are in general not on the lower mass 
cone, and so Xm will not converge to X. Because of these problems, we conclude that 
it is not admissible to first perform the perturbation expansion for tm and to multiply 
by Xm afterwards. Instead, the /^-dependence of X^. must be taken into account in 
the perturbation expansion. 

At this point it is very helpful that we stated the normalization condition for a chi- 
ral Dirac sea in the form (C.1.13). The key observation is that if we substitute (C.1.10) 
into (C.1.13), compute the operator product and take the limits m,m' \ 0, all con- 
tributions to (C.1.10) which are at least quadratic in m drop out. More precisely, if 
we expand Pm in the form 

= (xL ^ + ^ (1 - p) + I Pil^ + 0(m2)) 5ik^ - m") Q{-k^) , (C.1.15) 
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the error term is of no relevance for the normahzation condition (C.l. 13). Taking the 
inner product of (C.l. 8) with k and using the first part of (C.l. 4) together with the 
relation k"^ = m^, one sees that vk = 1. We use this identity in (C.l. 15) to obtain 

= (xL ^ + f + J P [y, ^]) - rr?) ei-k') . (C.1.16) 

Writing in this form has the advantage that we can pull out the chiral projectors 
by setting 

Pm = ^{Xim + imX*) (C.1.17) 

with X = XL and 

irn = + m + ^ p [y, ^]) Si^ - m^) e(-fc°) . 

Again neglecting terms quadratic in m, tm is a solution of the Dirac equation, 

{i0 + Bo-m)im = 0, (C.l. 18) 

where 

777 

Mk) = -^p\i.H\- (c.1.19) 

The formulation of the vacuum (C.1.17) and (C.l. 18, C.1.19) has the advantage that 
the interaction can easily be introduced. Namely, in order to describe the interaction 
we simply insert the external potentials into the Dirac equation (C.l. 18). In this way, 
the problems mentioned after (C.l. 14) get resolved. Instead of working with a k- 
dependent chiral asymmetry matrix X^, the /c-dependent vector field v in (C.l. 10) is 
now taken into account by a perturbation Bq of the Dirac equation, making it possible 
to apply perturbative methods in the spirit of §2.2. 

An obvious technical problem in this approach is that the perturbation operator 
Ho, (C.1.19), is not of a form previously considered in that it is nonlocal, is not causality 
compatible and does not decay at infinity. This problem will be analyzed in detail in 
Section C.2. What makes the problem tractible is that i3o tends to zero as m \ and 
is homogeneous, meaning that its kernel Bo{x, y) depends only on the difference x — y. 

Let us verify in which generality the above method (C.1.17, C.l. 18) applies. In the 
example (C.l. 12), we can write the chiral Dirac sea in the form (C.1.17) with 

im = + f ) - ^) Qi-k') , (C.1.20) 

and tjn is a solution of the Dirac equation (C.l. 18) with Bo = m/2. Thus in this case, 
Bq is a homogeneous local operator. More generally, the method of pulling out the 
chiral asymmetry (C.1.17) applies to any distribution of the form 

Pm{k) = (xL (odd) + (even) + 0(m^)) S{k^ -cm^) e(-A;°) , 

where "(odd)" and "(even)" refer to a product of an odd and even number of Dirac 
matrices, respectively (and c is a constant). Namely, the corresponding tm is 

im{k) = ((odd) + 2 (even) + 0{m^)) S{k^ -cm^) @{-k°) . 

Hence the only restriction of the method (C.1.17, C.l. 18) is that the right-handed odd 
contribution to Pm should be of the order 0{m?'). For example, our method does not 
apply to 

Pm{k) = {xLH + mxRfH + m + 0{m')) d{k^ - m^) Qi-k"^) 
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with a scalar function f{k), although in this case the normalization condition (C.1.13) 
is satisfied. Dropping this restriction would make it necessary to give up (C.1.17) and 

thus to treat the trace compatibility on a level which goes far beyond what we can 
accomplish here. It is our view that assuming that the right-handed odd contribution 
to Pm is of the order 0{m'^) is a reasonable technical simplification. 

We close our discussion with a comment on the example (C.1.12). Wc saw above 
can be written in the form (C.1.17) with tm according to (C.1.20), and that 
tm is a solution of the Dirac equation (C.1.18) with the perturbation Bq = m/2. An 
alternative point of view is that im is a solution of the free Dirac equation of half the 
mass, 

777 

{i^-M)tm = ^ with ^ = y- (C.1.21) 

We refer to the method of considering a Dirac equation in which the mass parameter 
is multiplied by a constant as the modified mass scaling. The modified mass scaling 
has the advantage that one can satisfy the normalization conditions for chiral Dirac 
seas (C.1.13) with according to (C.1.17) and im a solution of the free Dirac equa- 
tion. 

C.2. The Homogeneous Perturbation Expansion 

In the above examples we saw that there are different methods for giving a chiral 
Dirac sea a rest mass, which all correspond to inserting a suitable homogeneous oper- 
ator Bq into the Dirac equation. Furthermore, we found that the terms quadratic in 
the mass were irrelevant for the normalization of the Dirac sea, and this suggests that 
it should be possible to treat Bq perturbatively. This is indeed possible, as we shall 
now show for a general class of operators Bq. 

For simplicity, we again consider a single Dirac sea. We let Bq be a homogeneous 
operator, whose further properties will be specified below. In order to keep track of 
the different orders in perturbation theory, we multiply Bq by a small parameter e > 0. 
Exactly as in (2.6.10), we insert a parameter /Lt into the Dirac equation, which then 
reads 

{i^ + eBo- nl)^ = 0. 

Here the Dirac operator is homogeneous and is therefore diagonal in momentum space. 
Thus for given momentum k, the Dirac equation reduces to the 4x4 matrix equation 

{}ji + £Bo{k)-n)^{k) = 0. (C.2.1) 

Our aim is to introduce and analyze the spectral projectors and Green's functions of 
the Dirac operator i^ + eBq, where we regard /j, as the eigenvalue. In preparation, we 
shall now analyze the matrix equation (C.2.1) for fixed A; in a perturbation expansion 
to first order in e. If A;^ 7^ 0, the matrix ]^ is diagonalizable with eigenvalues and 
spectral projectors 

/x± = ±fik, E± = Ul±-^) , (C.2.2) 

where we set /x^ = Vk^ (if fc^ < 0, our sign convention is such that f^k lies in the 
upper complex half plane). The cigcnspaces lmE± are two-dimensional. The spectral 

projectors E± become singular as k'^ — > 0. The reason is that on the mass cone 
C := {k \ k'^ = 0}, the matrix ^ is not diagonalizable. We will address this problem 
later and for the moment simply assume that fc^ 7^ 0. We next consider the Dirac 
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operator ^ + eBq for small e. Perturbing the eigenspaces Im E± gives rise to two- 
dimensional invariant subspaces, and a standard calculation shows that the projectors 
Ej^ onto these subspaces are given by 

El = Es + {Es Bo E-s + E-s Bo E,) + ©(e^) (C.2.3) 

with s = ± and s = —s = ^. It remains to diagonalize the operator ^ + eBo on the 
invariant subspaces ImEg. This is carried out in the next lemma. We choose three 
(possibly complex) Lorentz vectors (q'j)i=i,2,3 such that 

{qi, k)=0 and qj) = -6ij . (C.2.4) 

More precisely, if k is timc-likc, wc choose the (qi) as a real orthonormal basis of the 
space-like hypersurface <A:>-'-. If on the other hand q is space-like, we choose qi and q2 
real and space-like, whereas ^3 is time-like and imaginary. We use the vector notation 
Q = {QIiQ2-,Qz) and introduce the matrices ^i,2,z by 

t = p^. (C.2.5) 

Lemma C.2.1. Suppose that k"^ ^ and that for small e, the matrix ^ + eBo is 
diagonalizable. Then its eigenvalues (fJ-s) s=±,a=y2 o,f^ given by 

tif = fik + £ ± T+) + 0{e^) (C.2.6) 

/x^ = -f^k + £ {i^- T r_) + , (C.2.7) 

where 

Us = \i:^{EsBo) (C.2.8) 

Ts = ^TY{tEsBo) (C.2.9) 

rs = V{r}? + irjf + {rlf . (C.2.10) 
The corresponding spectral projectors can he written as 

E^ = WEs + (n" E, Bo E-s + E-s Bo Es) + 0{e^) (C.2.11) 

2/^fe 

with 

= ^ (l ± ;^ Tsf^j . (C.2.12) 
If fs = 0, the invariant subspace ImE^ is an eigenspace to first order in e; i.e. 
{H + eBo)\i^Ei = {^l^k + ei^s)ni^Et +0{e''). 

Proof. We restrict attention to the invariant subspace Im E^ ; for Ei_ the proof is 
similar. A short calculation using (C.2.5) and (C.2.2, C.2.4) shows that 

[Si, S+] = , S2 = 1 , Tr(S, Sj- E+) = 2 5,j . 

This means that the matrices Sj are invariant on ImEj^, have the eigenvalues ±1 on 
this subspace and are orthogonal. Thus by choosing a suitable basis (and possibly 
after changing the orientation of S by exchanging Si with S2), we can arrange that 
the matrices S|iin£'_|_ coincide with the Pauli matrices ex. To first order in e, the eigen- 
values are obtained by diagonalizing Bq on the unperturbed invariant subspace Im Ej^ . 
A short calculation shows that the 2x2 matrix vl-^-ra has the eigenvalues vzizT and 
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corresponding spectral projectors IIi^ = ^(1 + 7 rff) with r = ^ {jiY + {T2Y + {j^Y- 
This gives (C.2.6, C.2.7). Finally, (C.2.11) follows from standard perturbation theory 
without degeneracies. ■ 



To avoid confusion, we point out that in general 7^ \ fs\ because (C.2.10) involves 
ordinary squares instead of absolute squares. In particular, it is possible that = 
although Tg 7^ 0. However, in this case the 2x2 matrix e^Bolim is not diagonalizable, 
and thus the above lemma does not apply. 

Remark C.2.2. In the proof of the previous lemma we used that the three matrices 

Silim E+ can be represented as the Pauli matrices cjj. It is instructive to verify explicitly 
that these matrices satisfy the correct commutation relations, for example 

2 Pl) ^2]|lmE+ = 5]3|iniE_,_ . 

We now give this calculation in detail. By a choice of coordinates, we can arrange 
that k = {u},p) and = {0,qi^). The standard identity between the Dirac matrices 

io'jk = 2 ^jklm <7^™' yields that (possibly after changing the orientation of E), 

Prom the definition of S, (C.2.5), one sees that [Si, II2] = — 2^i^2, and using (C.2.13) 
as well as the identity [/ifc,7°] = 0, we conclude that 

S2] = -^[^-^,, 7V (C.2.14) 

In order to simplify the rhs of (C.2.14) on Im E_^_, we use that satisfies the 
Dirac equation 

{li-lik)E+ = 0. (C.2.15) 
This identity allows us to replace the commutator with — n^hy an anti-commutator, 

W - fik, 7°] E+ = {]^- fik, 7°} E+ = {2u - 2fikl°) E+ . (C.2.16) 
Multiplying (C.2.15) by 2uj/iJ.k and adding (C.2.16) gives 

[^-^,,7°]i?+ = -U-^^Ae^. 

fJ-k \ J 

Using this identity in (C.2.14) gives 

2 Pl' ^2]|lmE+ = PMlmE+ 

with 

and a short calculation shows that this vector has indeed all the properties listed 
after (C.2.4). 

We shall now define the spectral projectors and Green's functions corresponding 
to the Dirac operator + eBq. We denote the spectrum of the matrix in (C.2.1) by 

a%k) = a{]^ + £Bo{k)) . 
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It is natural to define the spectrum cr^ of the Dirac operator + eBo as the union of 
the a^{k)s, 

= [j a%k) . 

As we saw above, the matrix ^ + eBo{k) in general is not diagonalizable, and thus 
we cannot introduce the spectral projectors for all k pointwise. But since the diag- 
onalizable matrices are dense in G1(C*^), it is reasonable to assume that the matrix 
}^ + £Bo{k) is diagonalizable for almost all (a.a.) k. Our formalism will involve mo- 
mentum integrals where sets of measure zero are irrelevant. Therefore, we may in what 
follows restrict attention to those k for which the matrix }j: + sBQ{k) is diagonalizable. 
Moreover, we shall assume that Bq is smooth and hounded. According to (C.2.2), the 
spectrum of the unperturbed Dirac operator is a^^^ = M U iM. The next lemma shows 
that the real part of the spectrum is stable under perturbations. 

Lemma C.2.3. Suppose that k'^ > 0. Then for e sufficiently small, cF^{k) C M. 

Proof. Choosing coordinates such that k = (w, 0), it is obvious that the eigenspaces 
of 1^ are definite, i.e. 

I 7^ for all eigenvectors 

By continuity, the eigenspaces oi 1^ + eBq (k) will also be definite for sufficiently small 
e. As a consequence, the corresponding eigenvalues are real, because 



Unfortunately, we have a-priori no control of how the imaginary part of the spec- 
trum changes with e. For this reason, it is most convenient to introduce the spectral 
projectors for all fi £ C, such that they vanish identically for ji ^ . For the nor- 
malization, we work with (^-distributions supported at one point in the complex plane. 
More precisely, we set 

5'^{z) = J(Re z) (5(Im z) 

I d^z ■■■ = I c?(Re z) d(Im z) ■ ■ ■ . 
Jc ilR2 

Dep. C.2.4. For h e C and keR'^ we set 

p^k) = E Kik)S'{l^-l^tik)) (C.2.17) 

s=±, a=l/2 

kf,{k) = e{k')f^{k) (C.2.18) 
r pp 

slik) = / pl{k). (C.2.19) 

JC IJ. — 

We also consider p^^, k^, and as multiplication operators in momentum space. 
In formal calculations, the operators and fc^ are solutions of the Dirac equation, 

{i^ + eBo - ^i) p'l, = = {i^ + eBo - ii) k% 
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and satisfy in analogy to (2.2.26-2.2.28) the multiplication rules 

Pi Pi' = KK' = '^'(/^-/^Op^ (C.2.20) 

as well as the "completeness relation" 

/ Pld'f^ = 1. 
Jc 

Using these identities in (C.2.19) yields that 

(i^ + £i3o-/x)4 = 1. 

Thus on a formal level, the operators p^, and are the spectral projectors and 
Green's functions of the Dirac operator, respectively. In order to give these operators 
a mathematical meaning, we can proceed as follows. Let k be such that the matrix 
^ + £Bo{k) can be diagonalizcd. Then the functional calculus for finite matrices (as 
defined e.g. via the approximation by polynomials) allows us to introduce for / G C^{C) 
the matrix f{l^ + eBo{k)). Formally, we can write the functional calculus with the 
spectral projectors, 

^ /(/i) Plik) d^fj. = fm + sBoik)) . (C.2.22) 

We can use this relation to give the integral in (C.2.22) a rigorous sense for a. a. k. 
The same argument applies to k^. For s^, we can similarly use the formal identity 

^ /(/x) s^ik) d'ti ^''='^ g{ix) pl{k) (C.2.23) 

with 

/" PP 

= / /(M)dV- 

Jc l^ — ^ 

In this way, one sees that the operators p^, A;^ and are well-defined when evaluated 
weakly in and k. 

Under additional assumptions, we can make sense of the operators in Def. C.2.4 
even for fixed real fi. We first justify the 5-distribution and the principal part. 

Lemma C.2.5. Suppose that for a given interval / C M, the spectral projectors 
in (C.2.17) are hounded uniformly in ^ ^ I and /c G M^. Then for a.a. (j, E I, the 
operators p^, fe^, and are well-defined distributions in momentum space. 

Proof. We write the Dirac equation {1^ + eBo{k))'^ = in the Hamiltonian form 

= H{LO,fi)-^ with H{oj,n) = --f'^ + eBo{uj,k) - nl) 

and k = {uj,k). In what follows we keep k fixed and consider this equation for variable 
parameters uj,n gM. The matrix H{uj,fi) is Hermitian with respect to the positive 
scalar product (.|.) = -<. 17*^1.;^. Thus it can be diagonalized; we denote its eigenvalues 
(counting multiplicities) by fii < • • • < ^4. The min-max principle (see [RS]) allows 
us to write 0„ as 

Qn = mill max , 

U,dimU=n u&U, \\u\\=l 
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where ||.|| is the norm induced by (.|.) and U denotes a subspace of C^. It fohows from 
this representation that the 0„ depend Lipschitz-continuously on lo and fx. Namely, 

f2„(a;) = min max ||ii'(a;) u\\ 

U,dimU=n u(^U, \\u\\ = l 

= min max \\H{uj') u + (H(uj) — H{uj')) u\\ 

U, dim. U=n u&U, \\u\\=l 

^ rrJ^^r ^ax ^11 + II^M - ^(^0 II 11^11) 

C7, dimC7=n ||u||=l 

= nn{uj') + \\H{Lo)-H{u;')\\. 
Using that Bo{k) is with bounded derivatives, we obtain the estimate 

\\H{lv) - H{oj')\\ < ||e7° (^o(c^) -^o('^'))ll < ec\LO-Lo'\ 

and thus r2n('^) — ^n{^') < ec|a; — uj'\. Exchanging the roles of u and u' gives the 
bound 

\nn{u}) - nn{uj')\ < ec\u-u'\. (C.2.24) 
A similar calculation shows that 

\nnifl) - nn{fl')\ < \n-fi'\. (C.2.25) 

We next consider for given n the equation 

w = nniio,n). (C.2.26) 

The following argument shows that for sufficiently small e, this equation has a unique 
solution ujn, which depends Lipschitz-continuously on fj,. Let (p (for fixed ji and n) be 
the mapping 

: R ^ M : a; HH- ri„(a;,/x) . 

According to (C.2.24), 

\<p{uj) — = \Cln{u>) — Cln{uj')\ < £c\uj — Uj'\ . 

Thus if wc choose e small enough, is a contraction. The Banach fixed point theorem 
yields a unique fixed point Un- The dependence on the parameter /x is controlled 
by (C.2.24) and (C.2.25). Namely, 

|a;n(M) - Wn(MOl = |f^n(^n(/i), /i) - ^^n('^n(/^')> /"')l 

< ec |Wn(/^) - t^n(/^')l + IA* - /^'l 

and thus 

\ujn{fi) - LOn{fi')\ < {1 - ec:)-^ \n - . (C.2.27) 

If we regard the spectral projector (C.2.17) as a distribution in lo, it is supported 
at those uj for which the Dirac equation + eBo — /x)^' = has a non-trivial solution. 
These are precisely the solutions Un of the equation (C.2.26). Thus we can write as 

4 



n=l 



duj{ii) 



(C.2.28) 



where the parameters a = a(n) and s = s{n) must be chosen such that Hg{ijJn) = At- 
Since a;„(/x) is Lipschitz (C.2.27), the factor |9^a;„(//)| in (C.2.28) is well-defined for 
a.a. fi and is uniformly bounded. Thus p^(a;) is a well-defined distribution for a.a. /x. 
The same argument applies to k^. 
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It remains to justify the Green's function s^. We can write it in the Hamiltonian 
framework as 

e ^ PP ^ PP 

l^ + eBo- fil cj - H{uj, ix) ^ ' 
Thus denoting the spectral projectors of H by (i^n)n=i,...,4) we have 

4 PP 

According to (C.2.24), ^ni'-^) is Lipschitz and thus differentiable almost everywhere 
with |c?a)^^n| < ec. The spectral projectors Fn{oj) can also be chosen to be Lipschitz. 
As a consequence, the principal part in (C.2.29) is well-defined for a.a. ji. ■ 

This lemma involves the strong assumption that the spectral projectors Eg must be 
uniformly bounded. We shall now analyze this assumption in detail. As one sees 
from (C.2.2) in the limit /i 0, the spectral projectors can have poles and thus in 
general are not uniformly bounded. Thus we need to impose an extra condition, which 
we will formulate using the following notion. 

Dep. C.2.6. Let A be a Ax 4 matrix, which is Hermitian (with respect to -<.\.y). 
A point ji G (y{A) is called e-definite if there is a subset a+ C (t{A) such that 

(i) The invariant subspace /_(_ corresponding to 0"+ is definite. 

(ii) dist{a+, o-{A) \a+) > e. 

Lemma C.2.7. If fJ- & ^^i-^) e-definite, the matrix A is diagonalizable on /+, and 
its spectral projectors Ea are bounded by 

Pall < , (C.2.30) 

where \\.\\ is a matrix norm and c is a constant which depends only on the choice of 
the norm \\.\\. 

Proof. It clearly suffices to consider a particular matrix norm. We introduce the 
positive scalar product (.|.) = -<;.|7''|.;^, let ||.|| = (.|.)2 be the corresponding norm, 
and set 

mil = sup \\A'^\\ . 

1' with |[1'||=1 

We denote the projector onto by E. E can be constructed with a functional calculus. 
Namely, let V{z) be a complex polynomial satisfying the conditions 

V\a+ = 1 and V\a_ = 0. 
Since these are at most four conditions, V can be chosen to be of degree three, 

3 

V{z) = ^c„z". 

n=0 

Furthermore, the fact that A is e-definite can be used to bound the coefficients Cn by 

|Cn| < ^ (C.2.31) 
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with a suitable constant C (this is easily seen from a scaling argument). The projector 
E is given hy E = V{A), and (C.2.31) gives the estimate 

' C „ , ^f\\A\^^ 



E\\ <Y.^\\Ar <Cr^] , (C.2.32) 



n=0 

where we used e < \\A\\ in the last step. 

By definition, Im = /+ is a definite subspace. We can assume without loss of 
generality that it is positive, i.e. 

-<^\E'^'^ > for all 

The matrix A\j_^ is Hermitian with respect to the positive scalar product -<.|.;^|/+. 
Thus it has a spectral decomposition with eigenvalues and corresponding spectral 
projectors Ea, a = 1, . . . ,N, 

n 

A\l+ = ^HaEa\l+ . 
a=l 

Extending the Ea by zero to the invariant subspace corresponding to cr{A) \ a+, the 
spectral projectors satisfy the relations 

N 

El = Ea = El, = , <^>\Ea^'^>Q ioi ah ^, 

a=l 

where the star denotes the adjoint with respect to ^.|.;^. 
We introduce the operators F and Fa by 

F = 7° E , F„ = 7° F„ . 

It is straightforward to check that these operators have the following properties, 

F+ = Fa, (^-IF^*) > (C.2.33) 

J^Fa = F, (C.2.34) 

a 

where denotes the adjoint with respect to (.|.). The relations (C.2.33) mean that 
the Fa are positive self-adjoint operators on a Hilbert space. This makes it possible to 
estimate the norm of the spectral projectors as follows, 

ll^all = ||7°i^a|| < ||7°ll lli^all < lli^all = SUp (* | F„ *) 

<!> with ||>I'||=1 

N 

< sup ^(^'IF;,*) = sup (*|F*) 

with ||>I'||=1 * with ||'I'||=1 

= ||F|| = ||7° E\\ < \\E\\ . 
We now apply (C.2.32) to obtain (C.2.30). ■ 



Dep. C.2.8. The Dirac operator + eBq has an e-definite kernel if for all 
fi G (— e^e) and all k with /x G (T^{k), /x is in the e-definite spectrum of the matrix 



Combining Lemma C.2.5 and Lemma C.2.7 gives the following result. 



220 C. NORMALIZATION OF CHIRAL FERMIONS 

Theorem C.2.9. // the Dirac operator + eBq has an e-definite kernel, then 
its spectral projectors and Green's functions (as given in Def. C.2.4) are well-defined 
distributions in momentum space for a.a. G {—e,e) . 

It remains to specify under which assumptions on Bq the Dirac operator has an 
£-definite kernel. We decompose Bq as 

Bo{k) = al + iPp + iZ + pcjl+^-^Wija'^ . (C.2.35) 

Here a, f3, v, a, and w are real potentials (namely the scalar, pseudoscalar, vector, 
axial, and bilinear potentials, respectively; clearly we assume w to be anti-symmetric). 
We introduce the function A(k) as the following combination of the axial and bilinear 
potentials, 

= -k'^ {a,a) + {a,kf-Wijk^ w'^k. (C.2.36) 
The first two summands can also be written as 

-k'^ {a,a) + {a,kf = -k'^ (^a- ^ {a,k) . (C.2.37) 

For timelike k, the vector inside the round brackets is spacelike, and thus (C.2.37) > 0. 
Similarly, the vector Wijk^ is spacelike for k timelike. We conclude that 

A{k) > if A;2 > 0. (C.2.38) 

Furthermore, A(g) vanishes on the mass cone C = {q^ = 0} if and only if q is collinear 
to the vector a and is an eigenvector of w, 

a = vq and UHjq> = Xqi {u,XeR, q€C). (C.2.39) 

Expanding (C.2.36), one sees that in this case, A is finite to the next order on the 
light cone, i.e. 

A{q) = ^ Z = hm A(A;) exists. (C.2.40) 

Qualitatively speaking, the next theorem states that the Dirac operator has an e- 
definite kernel if and only if the scalar potential is non-zero and dominates the axial 
and bilinear potentials. 



Theorem C.2.10. Suppose that for all qeC, 

Wij{q)a''q^ 



Hq)\ > I 



^ (i + e(i-2^7M) Vm\ */A(g) = o. 



i/A(g)/0 

(C.2.41) 



Then for sufficiently small e, the Dirac operator i^ + eBq has an e-definite kernel. If 

conversely there is q & C for which the opposite inequality holds (i.e. (C.2.41) with 
"> " replaced by "< "), then the Dirac operator has no e-definite kernel. 

Proof. A short calculation using (C.2.8, C.2.2, C.2.35) gives 

i^± = a± — {v,k). (C.2.42) 

In the special case 1j = Hk^^ and ^ = 7, we obtain furthermore from (C.2.9) that 

(r±)r = ar±Wro (r = 1,2,3). 
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Thus, according to (C.2.10), 

3 

{T±f = ^(Or)^ ±2ar WrO + {Wrof , 

r=l 

and this can be written eovariantly as 
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1 1 2 • • 

— {a, a) -\ — ^ {a, k)"^ ^ Wijk^ w^^ki ^ — WijO^k-' 



(C.2.43) 



This tensor equation is vahd for any time-hke k, and it is easy to check that it holds 
for spaceUke k as well. 

Let q & C. We first consider the case A{q) / 0. By continuity, A 7^ in a neigh- 
borhood U of q, and according to (C.2.38), A is positive in U. We substitute (C.2.42) 
and (C.2.43) into (C.2.6) and (C.2.7). In order to remove the singularities at /x^ = 0, 
we write the eigenvalues /i" in the form 

y/p + 2eS^ + £ (a ± K+) + 0(£2) 

= -y/A;2 + 2eSi^ + e (a t + 0{e^) ^ 

where we set ^ 

6^ = {v,k) ± A , K± = T± - — A. 



(C.2.44) 



The functions k± have the following expansion. 



A2 IF 2/Ltfe Wija^'W - A 



(C.2.45) 



In particular, one sees that these functions are bounded locally uniformly in //fc. Let 
us verify under which conditions the Dirac operator restricted to [7 has an e-definite 
kernel. Suppose that G (— e, e). Then, due to the square root in (C.2.44), 

fc2 + 2£(5i = 0{e^). 

It follows from (C.2.44) that 



Vfc2:f2^+0(e) = VO(£2) + 2£((52 - 5i) + 0{e) 
= V-4£ A + 0(e) - ^/s 

and therefore 

Moreover, we obtain from (C.2.44) and (C.2.45) that 

2/X+ + 2ea - e{K+ - /t_) + ©(e^) 



2iJ,l + 2e{a + ^ + 0{f^k 



Thus the condition |/x^ — /xL| > £ is satisfied if 



a + 



WiiO^k^ 



A 



> 



As is proved in Lemma C.2.11 below, the eigenspace corresponding to fi^ is definite. 
We conclude that ^\ is an e-definite eigenvalue of A. Repeating the above argument 
in the three other cases iJ?i,ii\. G (— e,e), one obtains that for sufficiently small e, the 
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kernel of the Dirac operator is e-definite in U. If conversely (C.2.41) holds with ">" 
replaced by "<", it is straightforward to check that the Dirac operator for small e has 
no £-definite kernel. 

It remains to consider the case A(gf) = 0. We write the eigenvalues as 



III = V^2 + 2£ {v,k) +e{a±T+) 



(C.2.46) 



According to (C.2.40), the first three summands in (C.2.43) have a finite limit at q. 
Furthermore, (C.2.39) yields that 

2 ■ ■ 

We conclude that the functions t± in a neighborhood of q have the expansion 

r± = ^f\l\ + 0{^\). (C.2.47) 

For small e, ~ ^/e, and so the term 0{yJjLk) is of higher order in e and can be 
omitted. Furthermore, the following continuity argument varying I shows that the 
eigenvalues /x" and /x" correspond to positive and negative eigenvectors, respectively: 

If / = 0, only the scalar and vector potentials enter the perturbation calculation to 
first order in e (sec (C.2.46, C.2.47)). If only scalar and vector potentials are present, 
the spectral decomposition of the matrix ^ + eBq is easily obtained from the identity 

[m + eif + ea) ~ eaf = {k + evf 1 . 

One sees that the eigenvalues are twofold degenerate, cr^{k) = and that if 

they are real, the corresponding eigenspaces are definite. The parameter I removes the 
degeneracy of these eigenspaces, but the resulting invariant subspaces remain definite. 

Suppose that /v,^ € (— e, e). We consider the two subcases 2-y/p| > e and 2iVjI| < e 
separately. In the first case, \iJ,\ — > e, and thus we must arrange that 

1/4 - > £ (2 VI > e). (C.2.48) 

In the second case, — /7,^| < e. Thus we must combine the eigenvalues to pairs and 
consider the definite eigenspaces corresponding to the sets as = s = ±, and 

must satisfy the condition 

dist (cr+,cr_) > e {2^fl<e). (C.2.49) 

Evaluating (C.2.48) and (C.2.49) using (C.2.46, C.2.47) and analyzing similarly the 
three other cases fiL,fi± G (—£,£) gives the condition (C.2.41). ■ 



Lemma C.2.11. Let A be a Hermitian matrix (with respect to -<.\.y). If fj- & cr(^) 
is real and the corresponding invariant eigenspace I is one- dimensional, then I is a 
definite eigenspace. 

Proof. Since each invariant subspace contains at least one eigenvector, / is clearly 
an eigenspace. We must show that / is definite. Assume to the contrary that / = <*> 
is null, i.e. 

^* = A* with A G M and -<;*|*:^ = 0. 
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We denote the invariant subspaces of A by iIfj,)^£a{A)- Since Ix = <^> is one- 
dimensional and null, there must be an invariant subspace 1^, /x 7^ A, which is not 
orthogonal to 

/^n<*>^ / 0. 

We choose on a basis (ei, . . . , e„) such that A is in the Jordan form, i.e. 

/ /i 1 ••• \ 
fi ■■■ 

1 

V ••• /X / 

Let G {1, . . . , n} be the smallest index such that ~<ek\^>- / 0. Then 

A-<efe|^'^ = ^EklA'i'y = ^Aetl'^y 

= /Li -<efc I + -<efc_i I = I2^ek\ *^ 

This is a contradiction. 



Suppose that the homogeneous operator Bq satisfies the condition (C.2.41) in The- 
orem C.2.10. Then the Dirac operator has an e-definite kernel. As a consequence, 
the distributions = ^(p^ — /cp are well-defined (see Def. C.2.4 and Theorem C.2.9). 
Following (C.1.17), we introduce the fermionic projector by 

P'^ = liXf^ + t^X*) (C.2.50) 
with X = XL- In order to analyze the normalization of P^, we consider the product 

P'^ P;, = \{Xt'^ tl, X*+X tl X f^, + X* tl, X*) . (C.2.51) 

According to (C.2.20) and (C.2.21), 

t^t^, = S^{fi-fi')f^. (C.2.52) 

Thus the only problem is to compute the products t^^Xf^^, and Using the 

relations XL / R = ^(1 T p) together with (C.2.52), this problem reduces to making 
mathematical sense of the operator product 

P ■ 

It seems impossible to give this expression a meaning without making additional as- 
sumptions on Bq. For simplicity, we shall impose a quite strong condition, which is 
motivated as follows. The spectral projectors corresponding to the unperturbed 
Dirac operator — fj, satisfy the relations pp^ p = P-^ and thus Pfj, pPfj, = (/x > 0) . 
It is natural to demand that the last identity should also hold in the presence of the 
homogeneous perturbation for small e. 

Dep. C.2.12. The kernel of the homogeneous Dirac operator + B{e,k) is e- 
orthogonal to p if for all p^p! G a"^(fe)n (— |, |), the corresponding spectral projectors 
E^{k) and Efj_i{k) satisfy the condition 

E^pE^, = 0. (C.2.53) 
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If the kernel of the Dirac operator is e-definite and £-orthogonal to p, it follows 
immediately that for all G (— |, |), 

tlptl, = 0. (C.2.54) 

Using (C.2.52) and (C.2.54) in (C.2.51), we see that 

Now we can take the limits £, /x \ to obtain 

lim lim ( P^P":, -- S^in- p') P^] =0. (C.2.55) 

In analogy to (C.1.13), this relation states that the fermionic projector is idempotent 
(apart from the factor \ which will be treated in Section C.3 using the modified mass 
scaling) . 

In the remainder of this section, wc analyze under which assumptions on So the 
kernel of the Dirac operator is £-orthogonal to p. We begin with a simple calculation 
in first order perturbation theory. 

Lemma C.2.13. Suppose that the Dirac operator i(p + ei3o has an e-definite kernel 
and that the homogeneous potentials in (C.2.35) satisfy for all k eM!^ the relations 

(3{k) = and eijim w'^ {k) k^ = 0. (C.2.56) 

Then for all k and p, p' G a^{k) fl (— |, |), 

E^{k)pE^{k) = Oie"). (C.2.57) 

Proof. Choose k and p,p' G a^{k) H (— |, |). Since the Dirac operator has an e- 
definite kernel, the invariant subspace / corresponding to the set {p, p'} C cr^{k) is def- 
inite (notice that p G (—£,£) and \p — p'\ < e). We saw in the proof of Theorem C.2.10 
that the invariant subspaces Im E^ and Im E^ (with E^ according to (C.2.3)) are 
definite. Thus / C Im E^ or / C Im Ef_. Therefore, it suffices to show that for all 
s = ±, 

ElpEl = 0{e'^) . (C.2.58) 

Substituting (C.2.3) and using the relations pE±p = E^:, we obtain the equivalent 
condition 

Es{Bo,p}Es = 0. (C.2.59) 
This equation means that the matrix {Bo,p} must vanish on the two-dimensional 
subspace ImEg. Since on this subspace, the matrices E, (C.2.5), have a representation 
as the Pauli matrices, we can restate (C.2.59) as the four conditions 

Tr{Es{Bo,p}) = = Tr(tE,{Bo,p}) . 

Evaluating these relations using (C.2.2, C.2.5, C.2.35) gives (C.2.56). ■ 

This lemma is not satisfactory because it gives no information on how the error term 
in (C.2.57) depends on k. More specifically, the error term may have poles on the mass 
cone (and explicit calculations show that such poles /c"^" indeed occur for n = 1 
and n = 2). Since in the limit e \ the kernel of the Dirac operator is the mass 
cone, it is far from obvious how to control the error term in this limit. In other words. 
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(C.2.57) cannot be interpreted as "the kernel of the Dirac operator is £-orthogonal to 
/9 up to a small error term." 

In order to resolve this difficulty, we must proceed non-pcrtTirbativcly. In general- 
ization of our previous ansatz ilp + eBq, we shall consider the Dirac operator + , 
where we assume that B^{k) is a homogeneous potential which is smooth in both 
arguments and has the power expansion 

B^k) = e Bo{k) + Bi{k) + B2{k) + ■■■ . (C.2.60) 

The higher order potentials Bi,B2,--- arc irrelevant for Def. C.2.8 because they arc 
negligible for small e. In particular, the statement of Theorem C.2.10 remains valid 
without changes. Furthermore, the potentials Bi,B2, ■ ■ ■ should be irrelevant for the 
statement of idempotence (C.2.55) because (C.2.55) involves a limit e \ 0. Therefore, 
it seems unnecessary to enter a detailed study of these potentials. The only point of 
interest is under which assumptions on Bq there exist smooth potentials Bi, B2,- ■ ■ such 
that the spectral projectors corresponding to the Dirac operator + B^ satisfy the 
conditions (C.2.53) exactly. 

Theorem C.2.14. Suppose that the Dirac operator i^+eBo has an e-definite kernel 
and that the homogeneous potentials in (C.2.35) satisfy for all k the relations (C.2.56). 
Then there is an e > and a smooth potential B^{k) having the expansion (C.2.60) 
such that the kernel of the Dirac operator i^ + B^ is e-orthogonal to p. 

Proof. Choose k and G '^'^{k) n (— Similar to what was described 
before (C.2.58), we know from the proof of Theorem C.2.10 that the matrix A = 
^ + B^{k) has a positive and a negative definite invariant subspacc, one of which 
contains Im£^^ U Im E^/. Again denoting the projectors onto these subspaces by 
and Et, respectively, it thus suffices to show that for s = it, 

£;fpEf = 0. (C.2.61) 

We first evaluate these conditions in a special spinor basis. Namely, we let ei and 62 
be an orthonormal basis of ImE^ and set 63 = pei, 64 = pe2. The conditions (C.2.61) 
imply that 63 and 64 span Im E^. Using the relation = 1 as well as that the 
subspaces <{ei,e2}> and <{e3,e4}> are invariant under A, we conclude that the 
matrices p and A are of the form 

where we used a block matrix notation corresponding to the splitting 

= <{ei,e2}>e<{e3,e4}> , 

and "*" denotes an arbitrary block matrix entry. Furthermore, the relation p* = — 1 
yields that 

-<e3 I 63^ = — 1 = -<e4 I e4>- , 
and thus the basis (ea) is pseudo-orthonormal, 
4 

^^r I $^ = ^ ^ with si = S2 = 1, S3 = S4 = -1- (C.2.63) 

a=l 

We see that the matrix p and the spin scalar product arc in the usual Dirac repre- 
sentation. In this representation, the fact that A is block diagonal (C.2.62) can be 
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expressed by saying that A must be a real linear combination of the 8 matrices 

1, 7°, p7, (C.2.64) 

We next express this result in a general basis, but again in the Dirac representa- 
tion. Since the representations of the matrix p, (C.2.62), and of the scalar product, 
(C.2.63), are fixed, the freedom in choosing the basis is described by even U{2, 2) trans- 
formations. This group, which wc denote by J7(2,2)''™", contains the normal Abelian 
subgroup U = {exp('!9/9/2) : -i? G M}. Acting by U on (C.2.64) gives the matrices 

1, ((coshi?-|-p sinh??)7°, ((cosh^9 + p sinh^9) p7, p7°7 . (C.2.65) 

When the factor group J7(2, 2)'^™YC/ acts on (C.2.65), the resulting transformations 
correspond precisely to Lorentz transformations of the tensor indices (see Lemma 1.2.1 
for details). Thus the conditions (C.2.61) are satisfied if and only if A is of the form 

A = al + {{cosh -& + p sinh '&)'!/ + {{p cosh -3 + sinh '3)4 + pt/p (C.2.66) 

with a time-like vector field u and two vector fields a and b, which are orthogonal to 
u, 

{u,a) = = {u,b) . (C.2.67) 

We substitute the identity A = ^ + B^{k) into (C.2.66) and solve for B^{k). Expanding 
in powers of £ gives the result. ■ 



C.3. The General Construction, Proof of Idempotence 

In this section we shall make precise what "idempotence" means for a fermionic 
projector with chiral asymmetry in the presence of a general interaction. We proceed 
in several steps. We begin with a straightforward extension of the results of Section C.2 
to systems of Dirac seas. Then we introduce the interaction and perform the causal 
perturbation expansion. After putting in an infrared regularization, we can define the 
fermionic projector. Finally, idempotence is established via a singular mass limit. 

We begin with a system of Dirac seas in the vacuum, described by the mass matrix 
Y and the chiral asymmetry matrix X (see §2.3). In order to give the chiral fermions 
a "small generalized mass," we introduce a homogeneous operator Bq and consider for 
£ > the Dirac operator iI^ + eBq — mY. For simplicity, we assume that Bq is diagonal 
on the sectors and is non-trivial only in the chiral blocks, i.e. 

{Bo)\i;] = 51 51 with 4'^") = if X, = 1. 

Then on each sector the methods of Section C.2 apply; let us collect the assumptions 
on Bq and the main results: For every index {aa) with ^ 1 we assume that 

(1) B^'^\k) depends smoothly on /c G M'^ and grows at most polynomially at 
infinity. 

(2) The (4 X 4)-matrix + £B^''°'\k) is diagonalizable for a.a. k. 

(3) bI)^°'^ has the decomposition into scalar, vector, axial, and bilinear potentials, 

i3("")(A;) = al + i + p4+'-^Wija'^ , 
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such that for all k G 

eijimw'^{k) k^ = 



and g G C the following conditions are satisfied, 



\a{q)\ > 



+ < 



Wij{q)a 



1 + 6(1 



{k I k 



if A(g)^0 

if A(g) = 0. 
^ = 0} is the mass 



A(9) 

2v1M) Vm\ 

(with A and I defined by (C.2.36) and (C.2.40), C 
cone). 

Then for sufficiently small e, the Dirac operator + eB^^ has an £-definite kernel 
(see Def. C.2.8 and Theorem C.2.10). Thus for a.a. ji G (— £,e), the spectral projec- 
tors pIi^'^"\ fc^'''""'' and the Green's functions s^'''""^ are well-defined distributions in 
momentum space (see Def. C.2.4 and Theorem C.2.9). Furthermore, the kernel of the 
Dirac operator is e-orthogonal to p (see Def. C.2.12 and Theorem C.2.14; for simplicity 
we here omit the higher order potentials 81,82, ■■ ■ in (C.2.60), this is justified because 
these potentials obviously drop out in the singular mass limit), and this can be stated 
in the form (cf. (C.2.17)) 



e,(aa) 







for all fi, fj' G 



'2' 2^ 



We build up the spectral projectors p^^, k^^ and the Green's function of the 
whole system by taking direct sums; namely, 




if = 1 



(C.3.1) 



where A stands for p, k, or s. Note that in the chiral blocks the mass parameter ^ 
(and not /i) is used. The purpose of this modified mass scaling is to get rid of the 
factor ^ in the normalization of a chiral Dirac sea (C.2.55) (also see the paragraph 
after (C.1.21)). The corresponding Dirac operator is 

+ e8o — mY — fj,Z , 

where the matrix Z = ^(X + X*) takes into account the modified mass scaling. The 
spectral projectors satisfy the multiplication rules 







S^{p-p') z-^p%^ 



1 ue 



-- 5\p-p')Z 
for 11,11' G (-§,§ 



where C stands for k or p. The Green's functions satisfy the relations 

PP 

PP 



p-ll- 



z-^ {g 



6 



Sj-i,') 



> 



(C.3.2) 
(C.3.3) 

(C.3.4) 



/i — /i' 

These multiplication rules differ from those in §2.2 only by the additional factor Z^^. 

To describe the interaction, we insert a potential 8 into the Dirac operator, which 
then reads 

i^ + B + eBo-niY - fiZ . (C.3.5) 
We assume that Y and 8 have the following properties: 



pe,S 
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(a) Only the chiral particles are massless, i.e. 

y{'i")>0 ifX„ = l. 

(b) B is the operator of multiplication with the Schwartz function B{x). 

(c) Y and B are causality compatible, i.e. 

X* {i0 + B -niY) = {i0 + B -niY) X . (C.3.6) 

In order to introduce the spectral projectors with interaction and k^^, we take 
the operator expansion of causal perturbation theory §2.3 and replace the operators 
according to A ^ A^_^^ (with A = p, k, or s). All the operator products of the resulting 

expansion are well-defined for a. a. /i (note that B{k) has rapid decay and >1^^(A;) grows 
at most polynomially at infinity). 

For the infrared regularization, we proceed exactly as in §2.6 and replace space by 
the three-dimensional torus (2.6.2). Furthermore, we "average" the mass parameter 
/X. More precisely, combining (2.6.13) with (C.2.50), the auxiliary fermionic projector 
is defined by 

If {Xi%^ + i%^X*)d'i,, (C.3.7) 

where as usual = |(p+^ — k^^l)■ Finally, the regularized fermionic projector is 
obtained by taking the partial trace, 

^pe,S^a ^ J]^(pM)(-). (C.3.8) 

a=l 13=1 

Before we can prove idempotence, we need to impose the following extension of 
the non-degeneracy assumption (2.6.1). We set 

^taa) = <^{H + eB^r^). 

Def. C.3.1. The Dirac operator i^ + eBq — mY has £-non-degenerate masses 

if for all a and /3 7^ 7, 

'^wn(-f,|)7^0 =^ 4^)n(-£,6) = 0. (C.3.9) 

Roughly speaking, the next theorem states that the masses are e-non-degenerate 

if they arc non-degenerate in the massive sectors and if the homogeneous potentials in 
the chiral sectors are sufficiently different from each other. 

Theorem C.3.2. Suppose that Y and Bq have the following properties: 

(i) In the massive blocks (i.e. Xa = Ij, the masses are non-degenerate, 

ym ^ y(b7) iff3j^j. 

(ii) In the chiral blocks ( i. e. Xa 7^ 1^, for all P ^ j and all q E C either 
{v^^^\q) + s A^^P\q) 7^ {v^^\q) + s' A^^\q) for all s, s' G {±1} (C.3.10) 

or else 



\a 



{q) - a^''^\q)\ > 2 + 2 \d^''^\q) + d^''^\q)\ (C.3.11) 
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with 



Wij{q)a 



z/A(g)/0 



d{q) = { A{q) 

. Vm\ ^/A(g) = 

(and A, I according to (C.2.36) and (C.2.40)). 
Then for sufficiently small e, the Dirac operator has e -non- degenerate masses. 

Proof. The condition in (i) follows immediately from the fact that the eigenvalues 
in the two sectors differ precisely by m{Y^^^^ — Y^^^^). For part (ii) we consider the 
formulas for the eigenvalues (C.2.44) and (C.2.46). If (C.3.10) holds, the eigenvalues 
in the two sectors all differ by contributions of the order y^, and so (C.3.9) is satisfied 
for small e. If on the other hand (C.3.10) is violated, there are eigenvalues in two 
different sectors such that the square roots in (C.2.44) and/or (C.2.46) coincide. Thus 
these cigcm -alucs differ by (a + a)^'"^^ - (a + a)^''''\ where each a is an element of the 
set {ibK+, ibK_, ibT+, zbr-}. The condition (C.3.11) guarantees that this difference is 
greater than 2e, and so (C.3.9) is again satisfied. ■ 



We can now state the main result of this chapter. 



Theorem C.3.3. (Idempotence) Consider the Dirac operator (C.3.5) under the 
above assumptions (l)-(3) and (a)-(c). Assume furthermore that the masses are e- 
non-degenera,te (see Def. C.3.1 and Theorem C.3.2). Then the corresponding fermionic 
projector (C.3.7, C.3.8) satisfies the identity 



lim lim 5 
£\0 6\o 







(C.3.12) 



with convergence as a distribution to every order in perturbation theory. 



Proof. Similar to (C.2.19), the Green's function s^^ has a spectral representation 

in a mass parameter v. We want to decompose ,s^^ into contributions and s"^^ 
where — /i| is small and large, respectively. To this end, we introduce in each sector 
the operator 

■p^^d'v 



.6,{aa) 



L 
L 



pp 



B./2(m) 



/i 



d^u 



if Xi 7^ 1 
if Xi = 1 



and define s^^ by taking as in (C.3.1) the direct sum. Setting 
obtain the decomposition 

= ^+11 + • (C.3.13) 
Our first step is to show that for small /i, the matrix s+^(A;) is bounded, more 
precisely that 



\s%,m < ^ for^e(-i,i) 



(C.3.14) 



with C{k) a smooth function with at most polynomial growth at infinity (the exponent 
7 is probably not optimal, but (C.3.14) is sufficient for our purpose). It clearly suffices 
to prove (C.3.14) in a given sector (aa); for simplicity the sector index will be omitted 
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(i.e. Bq = ;Sq""^). Furthermore, we only consider the case Xa / 1; the other case is 
analogous (and even simpler, because in the massive sectors no homogeneous potentials 
are present). We introduce the projector E{k) by 

E = J2 with S={ia, s) with - < £} 

(o,s)e5 

According to Lemma C.2.7, 

ll^ll < (C.3.15) 

with C\ (k) smooth with at most polynomial growth at infinity. The matrix has a 
simple spectral representation. 



Unfortunately, this representation is not suitable for estimates, because we have no 
control of H-E" || for (a, s) ^ S. To avoid this problem, we rewrite as follows. 



I Y K + — ^ — E<i I {i-E) 

\{a,s)^S ^ {a,s)£S ^ ' / 



= + sBo - + sE)-^ {1 - E) . 

Introducing the "Hamiltonian" H = —7° (^ + £^0(0;, fc) — /x + eE), we obtain 

= {u-H)-^j^{l-E). (C.3.16) 

The matrix H{k) is Hermitian with respect to the positive scalar product (.|.) = 
^.|7°|.;^ and can thus be diagonalized, i.e. 

4 

H = ^ Fn 

n=l 

with real eigenvalues 0„ and spectral projectors Substituting into (C.3.16) gives 

4 



si.. = j:-;^Fnl'il-E), 



n=l " 



and (C.3.15) yields the bound 

pi„|| < 2max^ ^-—\\1-E\\ < ^^max-^ ^— . (C.3.17) 

It remains to estimate the factors \uj — from below. We use that the determinant 
is multiplicative to obtain 

4 

JJ(a;-Q„) = det(a;-/f) = det(7° (a; - iJ)) 

n=l 

= detm + eB^-ix + eE) = (^u^ - ^) [] (/i^ - + e) . 

(a,c)^S {a,c)£S 
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Taking the absolute value, the factors — //| and 1//" — + e| are all greater than |, 
and thus 

^|a;-^^n| > (I)'. (C.3.18) 

n=l 

Since the eigenvalues of the Hermitian matrix oj — H can be estimated by the sup-norm 
of the matrix, 

we can deduce from (C.3.18) that each factor \u! — ri„| is bounded by 

\LO-nn\ > \\iO-H\\-^ (^^y . 

Substituting this inequality into (C.3.17) gives (C.3.14). 

The causal perturbation expansion expresses as a sum of operator products of 
the form 

i%, X Bo ■ ■ ■ Bo , (C.3.19) 

where each factor A stands for p, k, or s. Since Bo{k) has rapid decay and A^^{k) grows 
at most polynomially, these operator products are wcU-dcfincd. According to (C.3.7) 
and (C.3.8), the first summand inside the brackets in (C.3.12) can be written as 
^5 s 

- dfJ, dfl'^^ + (*+M)(fe^) (*+/i')(cI) 

b /3,7 

+ (*+m)(6^) -^b (*V)(c]j + (^+;i)(6^) ^b ^b (i+;i')(c]j) • (C.3.20) 

When we substitute (C.3.19) into (C.3.20), the difficult point is to multiply the right- 
most factor A of the first factor t to the leftmost factor A of the second factor t. More 
precisely, we must analyze the following operator products, 

(•••^;.)&1p(^V---)S (C.3.22) 

with A = p, k, or s. 

If one of the factors A in (C.3.21) or (C.3.22) is the Green's function, we substi- 
tute (C.3.13) and expand. Since is bounded (C.3.14), the products involving 
have a finite limit as 5 \ 0. Since the two integrals in (C.3.20) give a factor 5^, these 
products all drop out when the limit 5 \ is taken in (C.3.12). Thus it suffices to 
consider the case when the factors A in (C.3.21) and (C.3.22) stand for p, k, or s. 

Since the Dirac operator has e-non-degenerate masses, the distributions A^^(fe) 
have disjoint supports in different sectors. More precisely, for all /it, |Lt' G (— |, |), 

supp {A'/^^^) n supp (A'f^^) =0 if /? / 7 and / 1, 

where each factor A stands for p, k, or s. A similar relation holds in the massive 
blocks. Therefore, (C.3.21) and (C.3.22) vanish if /3 7^ 7. 

In the case /3 = 7, (C.3.22) is zero because the Dirac operator is e-orthogonal to 
p (C.3.3). Thus, using a matrix notation in the sectors, we only need to take into 
account the operator products 

(••• A-+^)(AV •••) 
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with A = p, k, oi s (here we may again consider s instead of s because, as we saw 
above, all factors s drop out in the limit 6 \0). Now we can apply the multiplication 
rules (C.3.2) and (C.3.4). Applying (C.3.2) gives a factor (5^(/7, — //'), and we can carry 
out the /i'-integral. After dividing by 6, we can take the limits \ and e \ 0. 
Using that in this limit the Dirac operator is causality compatible (C.3.6), we can 
"commute X through" the resulting operator products (see §2.3). In this way, one 
recovers precisely the unregularizcd fcrmionic projector P = limg^^x^o P^'^- If (C.3.4) 
is applied, the resulting principal part is bounded after the integrals over /j, and /x' are 
carried out, and we can take the limits 6 \ and e \ 0. After commuting X through 
the resulting operator products we find that all terms cancel. ■ 

For understanding better what the above results mean physically, it is instructive 
to consider a cosmological situation where the 4-volume of space-time is finite. In this 
case, the limits e, (5 \ in (C.3.12) are a merely mathematical idealization correspond- 
ing to the fact that the size of the universe is very large compared to the usual length 
scales on earth. We can extrapolate from (C.3.7, C.3.8) to get some information on 
how the properly normalized physical fermionic projector should look like: The pa- 
rameter 6 is to be chosen of the order with T the lifetime of the universe (also 
see §2.6). Then, due to the /x-integral in (C.3.7), the Dirac seas are built up from those 
fermionic states whose momenta lie in a thin strip around the mass cone. Naively, the 
modified mass scaling implies that for neutrinos this strip must be thinner. However, 
this naive picture is misleading because the detailed form of the chiral Dirac seas de- 
pend strongly on the homogeneous operator Bq, which is unknown. We point out that 
in (C.3.12) the order of limits is essential: we must first take the infinite volume limit 
and then the limit e \ 0. This means for our cosmology in finite 4-volume that the 
homogeneous perturbation eBq must be large compared to T~^. One possibility for 
realizing this is to give the neutrinos a small rest mass. But, as shown above, the same 
can be achieved by more general, possibly nonlocal potentials which do not decay at 
infinity. 
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The Regularized Causal Perturbation Theory 

In §4.5 we gave a procedure for regularizing the formulas of the light-cone expan- 
sion (4.5.5-4.5.10). We shall now derive this regularization procedure. The basic idea 
is to extend the causal perturbation expansion of §2.2 to the case with regularization, 
in such a way that the causality and gauge symmetry arc preserved for macroscopic 
perturbations. Using the methods of §2.5 one can then analyze the behavior of the 
so-regularized Feynman diagrams near the light cone. For simplicity, we will restrict 
attention to the first order in perturbation theory. But our methods could be applied 
also to higher order Feynman diagrams, and the required gauge symmetry suggests 
that our main result, Theorem D.2, should hold to higher order in perturbation theory 
as well. 

We first state our assumptions on the fcrmionic projector of the vacuum. As in 
Chapter 4 we describe the vacuum by a fermionic projector P{x, y) of the form (4.1.3) 
with vector-scalar structure (4.1.5). For small energy-momentum, P should coincide 
with the unregularized fermionic projector of the vacuum, i.e. 



P(fe) = (^ + m)(5(/c2-m2)e(-A;°) il \k^\ <^ Ep and \k\ ^ Ep. (D.l) 



Furthermore, we assume that the vector component is null on the light cone (i.e. that 
(4.4.21) holds with Sghear ^ 1)) and that P satisfies all the regularity assumptions 
considered in §4.3 and §4.4. For simplicity, we finally assume that P is supported 
inside the lower mass cone. 



(with according to (5.6.15)). This last condition is quite strong, but nevertheless 
reasonable. In particular, it is satisfied when P is composed of one-particle states 
which are small perturbations of the Dirac eigenstates on the lower mass shell. 

In this appendix we shall address the question of how one can introduce a classical 
external field into the system. For clarity, we will develop our methods mainly in the 
example of an electromagnetic field. As described in §4.1, we consider the regularized 
fermionic projector as a model for the fermionic projector of discrete space-time. In 
this sense, the regularization specifies the microscopic structure of space-time. Follow- 
ing the concept of macroscopic potentials and wave functions introduced in §4.1, the 
electromagnetic field should modify the fermionic projector only on length scales which 
are large compared to the Planck length, but should leave the microscopic structure 
of space-time unchanged. In order to fulfill this requirement, we impose the following 
conditions. First of all, we assume that the electromagnetic field be "macroscopic" in 
the sense that it can be described by an electromagnetic potential A which vanishes 
outside the low-energy region, i.e. 



supp P d C 



(D.2) 



A{k) = 



unless \k^\ < Ep and \k\ < Ep, 



(D.3) 
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where A is the Fourier transform of A. We denote the fermionic projector in the pres- 
ence of the electromagnetic field by P[4\- In order to prevent that the electromagnetic 
potential might influence the microscopic structure of space-time locally, we demand 
that A can locally be made to zero by a gauge transformation. Thus we impose that 
the usual behavior under U{1) gauge transformations 

P[4 + (^A)] (x, y) = e^^(^) P[4] (x, y) e'^^^^) (D.4) 

(with a real function A) should hold also for the regularized fermionic projector, as- 
suming that the involved potentials A and {A + 9A) arc both macroscopic (D.3). We 
point out that, because of the gauge symmetry in discrete space-time (following from 
the freedom in choosing the gauge (3.3.8)), the local phase transformations in (D.4) 
are irrelevant in the equations of discrete space-time, and thus the transformation 
law (D.4) implies the freedom to transform the electromagnetic potential according to 
J]L ^ JjL + ^A. Finally, we must rule out the possibility that the electromagnetic po- 
tential might influence the microscopic structure of space-time in a nonlocal way. For 
this purpose, we impose that the perturbation expansion for the regularized fermionic 
projector be causal, in the sense introduced in §2.2. 

Let us consider how these conditions can be implemented in the perturbation the- 
ory to first order. We first recall the standard perturbation theory for Dirac eigenstates. 
For a solution ^ of the free Dirac equation (i^ — m)^ = 0, the perturbation to first 
order in A, which we denote by A^'[^], is given by 

A^4]{x) = - J d%sm{x,y) 4{y) ny) , (D-5) 

where Sm{x,y) is the Dirac Green's function (2.2.18), 

sm{x, y) = 1-0^, i¥ + m) e-^'(^-y^ . (D.6) 

If we consider Sm{x,y) as the integral kernel of an operator Sm and the potentials as 
multiplication operators, we can calculate in the case 4 — 9^ to be 

A^'[^A] = (^A) ^ = isra [i^ - m. A] ^' 

= - m)sm) A"^ - is„ A ((i^- m)^') = iA 'I' . (D.7) 

Thus in this case, A^{x) = iA(x) ^(x) is simply the contribution linear in A to 
the phase transformed wave function exp(iA(x)) ^(x); this shows exphcitly that the 
perturbation calculation is gauge invariant. 

As a consequence of the rcgularization, the fermionic projector P{x, y) is in general 
not composed of Dirac eigenstates. Therefore, we next consider a wave function 
which is not necessarily a solution of the free Dirac equation. But according to (D.2), 
we may assume that its Fourier transform 4f has its support in the interior of the mass 
cone, 

supp C {k\k^ >0} . (D.8) 

In this case we can introduce A^'[4] as follows. The spectral projector of the free 
Dirac operator corresponding to the eigenvalue /x G M has the form 

P,{x, y) = 1^ + ^) ^i^' - ^"^'^'"'^ (D-9) 

(see (2.2.4); notice that wc added the step function e(^) to allow for the case < 0). 
Since the real axis is only part of the spectrum of the free Dirac operator (namely, the 
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free Dirac operator has also an imaginary spectrum), the spectral projectors (p/i)/ieiR 
are clearly not complete, i.e. J^^Pudfx ^ 1. By integrating (D.9) over /x, 

p,{x,y)di, = |^G(fe2)e-'=M, (D.IO) 

one sees more precisely that the operator p^dii is the projector on all the momenta 
in the mass cone. But according to (D.8), ^ lies in the image of this projector, and 
we can thus use the spectral projectors p^ to decompose ^' into eigenstates of the free 
Dirac operator. Each eigenstate can then be perturbed using (D.5). This leads us to 
introduce A^'f^] according to 

/oo 
diis^4p^^. (D.ll) 
-oo 

This definition of A^* shows the correct behavior under gauge transformations; namely, 
similar to (D.7), 

/oo 
dn Sf, [i^ - n, A] Pf, ^ 
-oo 

= p.d^y ^""-''J"-'^ (D.12) 

Thinking in terms of the decomposition (2.2.1) of the fermionic projector into the 
one-particle states, it seems natural to introduce the perturbation of the fermionic 
projector AP[4] by perturbing each one-particle state according to (D.ll). This leads 
to the formula 

/oo 
dii{s^4p^P + Pp^JfiS^) . (D.13) 
-oo 

The gauge symmetry can again be verified explicitly. Namely, a calculation similar to 
(D.12) using (D.2) yields that 

AP[^A](x,y) = iK{x) P{x,y) - iP{x,y) K{y) , 

and this is the contribution linear in A to (D.4). The perturbation calculation (D.13) 
is immediately extended to a general perturbation operator B by setting 

/oo 
df^is^Bp^P + Pp^Bs^) . (D.14) 
-oo 

Let us verify if the perturbation calculation (D.14) is causal in the sense of §2.2. 
Since it seems impossible to write (D.14) in a manifestly causal form, we apply here 
a different method, which allows us to analyze the causality of the perturbation ex- 
pansion in momentum space. As mentioned in §2.5, the causality of the perturbation 
expansion can be understood via the causality of the line integrals over the exter- 
nal potentials and fields which appear in the light cone expansion. More precisely, 
causality means that the light-cone expansion of AP(x, y) should involve only line 
integrals along the line segment xy, but no unbounded line integrals like for example 
Jq°° dX B{Xy + (1 — A)a;). This way of understanding the causality of the perturbation 
expansion yields a simple condition in momentum space. Namely, if B has the form 
of a plane wave of momentum g, i.e. B{x) = BqCxp{—iqx), then the unbounded line 
integrals become infinite when q goes to zero (for By fixed), whereas integrals along 
the line segment xy are clearly bounded in this limit. Hence we can say that the per- 
turbation calculation (D.14) is causal only if it is regular in the limit q ^ 0. In order 
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to analyze this condition, we substitute the exphcit formulas (D.6, D.9) into (D.14) 
and obtain 



/CO r 
-oo J 



(27r) 



4 



PP 

(^^^^2_^2 {^ + i + Pi)B,{li + fi) 5{e - fi') P{k) e-^(^+«)-+^'=s' 

We set q = sq with a fixed vector q and consider the behavior for e \ 0. Taking only 
the leading order in e, one can easily carry out the /^-integration and gets 

J4i 



AP[B]{x,y) = -\ j ^e-^'^M 



^^{^B,+B,^) p{k) 



2kq + eq^ 



PP \ 

+ P{k) 0^ B, + B, z^^^q— 2 J + 0{e^) . (D.15) 



Since 



PP PP PP 

lim — ^ = lim 



2kq + eg2 2kq — eq^ 2kq 

in the sense of distributions in the argument kq (notice that this kind of convergence 
is sufficient using the regularity of P), the leading singularity of (D.15) for e \ has 
the form 

Taking the Fourier transform in the variable (x — y), it is clear that (D.16) vanishes 
only if the commutator /anti-commutator combination [{Bg, I^}, P{k)] is zero for all 
A;. Since the perturbation Bq can be arbitrary, one sees (for example by considering 
a scalar perturbation, Bq ~ 1) that it is a necessary condition for the perturbation 
calculation (D.14) to be regular in the limit q ^ that 

[}^, P{k)] = for all k. (D.17) 

This commutator vanishes only if the vector field v{k) in (4.1.5) is a multiple of k, or, 
using the notation of §4.4, if the surface states have no shear. We conclude that the 
perturbation calculation (D.14) is in general not causal. 

Before resolving this causality problem, we briefly discuss how this problem comes 
about. The condition (D.17) can be stated equivalently that the operator P must com- 
mute with the free Dirac operator. In other words, the perturbation calculation (D.14) 
is causal only if the fermionic projector of the vacuum is composed of eigenstates of 
the free Dirac operator. In this formulation, our causality problem can be understood 
directly. Namely, since our perturbation method is based on the perturbation calcula- 
tion (D.5) for Dirac eigenstates, it is not astonishing that the method is inappropriate 
for non-eigenstatcs, because the perturbation expansion is then performed around the 
wrong unperturbed states. It is interesting to see that this shortcoming leads to a 
breakdown of causality in the perturbation expansion. 
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In order to comply with causality, we must modify the perturbation calculation 
(D.14). Our idea is to deduce the perturbation calculation for the fermionic projector 
from that for a modified fermionic projector, which satisfies the causality condition 
(D.17). The simplest idea for modifying the fermionic projector would be to introduce 
a unitary transformation U{k) e U (2, 2) which makes the vector v{k) in (4.1.5) parallel 
to k, more precisely 

U{k)-^ Vj{k)-f^ U{k) = \{k)li withA(fc)GM. 

However, a unitary transformation is too restrictive because it keeps the Lorentzian 
scalar product v{k)'^ invariant, and thus it cannot be used for example in the case 
when v{k) is space-like, but k is time-like. Therefore, we shall consider here a linear 
combination of unitary transformations. More precisely, we introduce for L > 1 and 
/ = 1, . . . , L unitary operators Ui{k) G U{2, 2) and real coefficients q such that 

L L 

J2ci{k) = l and vj{k)j^ = ^ci{k)Ui{k)X{k)Ij^Ui{k)-' (D.18) 
1=1 1=1 

with A(A;) G R. The existence of {Ui,ci) is guaranteed by the fact that the U{2,2) 
transformations comprise Lorentzian transformations §1.5. Clearly, the representation 
(D.18) is not unique. According to (D.l), we can choose the transformation (D.18) to 
be the identity in the low-energy region, and can thus assume that 

Ui{k) = 1 if \k^\ < Ep and \k\ < Ep. (D.19) 

Furthermore, the regularity assumptions and the particular properties of the fermionic 
projector mentioned before (D.2) give rise to corresponding properties of the operators 
[/;; this will be specified below (see (D.31, D.52)). The operators obtained by multipli- 
cation with Ui{k) in momentum space are denoted by Ui; they have in position space 
the kernels 

Ui{x,y) = j -0-^U,{k)e-'^^--y\ (D.20) 

Finally, we introduce the "modified fermionic projector" Q by replacing the vector 
field v{k) in (4.1.5) by A(A:) ^, i.e. 

Q{k) = (A(/c) ^ + m 1) f{k) . (D.21) 

According to (D.18), the fermionic projector P is obtained from Q by the transforma- 
tion 

L 

P = ^ciUiQUr'. (D.22) 
1=1 

The modified fermionic projector (D.21) satisfies the condition [Q[k), Jj(] = 0. 
Hence the perturbation calculation for Q does not suffer from our above causality 
problem, and we can introduce AQ[B] in analogy to (D.14) by 

/oo 
d/x {s^Bp^Q + Qp^B s^) . (D.23) 
-oo 
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We now deduce the perturbation of P by applying to (D.23) a transformation analogous 
to that in (D.22), namely 

L 

AP[B] := J2ciUiAQ[B]Ui-' (D.24) 
1=1 

^ /"OO 

= - 51 / d^iUiis^Bp^Q + Qpf^B s^) Uf^ . (D.25) 

1=1 

This last transformation should not affect the causality (in the sense of §2.2) because if 
(D.23) is regular when the momentum q of the bosonic potential goes to zero, then the 
transformed operator (D.24) will clearly also be regular in this limit. We call (D.25) 
the regularized causal perturbation of the fermionic projector to first order. 

The perturbation calculation (D.25) requires a detailed explanation. Qualitatively 
speaking, the difference between (D.14) and (D.25) is that the spectral projectors p^, 
the Green's functions s^, and the perturbation operator B have been replaced by the 
unitarily transformed operators 

UiP^,U-\ Uis^U-^ and UiBU'^ (D.26) 

and that a linear combination is taken. According to (D.19), the unitary transforma- 
tions in (D.26) have no influence on the macroscopic properties of these operators, i.e. 
on the behavior when these operators are applied to wave functions with support in 
the low-energy region. But the transformation (D.26) changes the operators on the 
microscopic scale, in such a way that causality is fulfilled in the perturbation expan- 
sion. We point out that in the case where B is the usual operator of multiplication 
with the external potentials, the transformed operator UiBU^^ is in general no longer 
a multiplication operator in position space; thus one can say that the classical poten- 
tials have become nonlocal on the microscopic scale. In order to better understand 
why the causality problem of (D.14) has disappeared in (D.25), it is useful to observe 
that Q commutes with the spectral projectors p^. This means that Q is composed of 
eigenstates of the Dirac operator, so that the perturbation expansion is now performed 
around the correct unperturbed states. 

Let us consider a gauge transformation. In the case B = ^A, the perturbation 
(D.25) is computed to be 

^ j'OO 

AP[^A] = ij^ci dixUi {Sf, [i^ - fi, A]p^Q + Qp^ [i^ - /x, A] s^) Uf^ 
1=1 

= i^ci diiUiiAp^Q - Qp^A) Uf^ 
1=1 

= 5^ q (^UiA p^ di^ QUf^ - iUiQ (^J°° p^ df^j AUf^^ . (D.27) 

By construction oi Q, we can assume that the distributions P and Q have the same 
support, and thus (D.2) holds for Q as well, 

supp Q C . (D.28) 
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Hence, according to (D.IO), the projectors J^^Pudju, in (D.27) can be omitted, and 
we conclude that 

L 

AP[^A] = Y,ci{iUihU^^UiQU^^ - iUiQUf^ UiAUf^) . (D.29) 
1=1 

If in this formula we were allowed to replace the factors UihU^^ by A, we could 
substitute in (D.22) and would obtain the contribution linear in A to the required 
transformation law (D.4). Indeed, the difference between A and UiAUf is irrelevant, 
as one sees in detail as follows. We consider one summand in (D.29) and set for 
ease in notation U = Ui. According to (D.19), the operators A and UAU~^ coincide 
macroscopically (i.e. when applied to functions with support in the low-energy region), 
and thus (D.29) yields gauge symmetry on the macroscopic scale. However, such a 
macroscopic gauge symmetry is not sufficient for us: to ensure that the microscopic 
structure of space-time is not influenced by the electromagnetic field, it is essential that 
(D.4) holds even on the Planck scale. In order to show microscopic gauge invariance, 
we consider the operator UAU~^ in momentum space, 

{UAU-' f){q) = / U{q) Hq ' P) Uip)'' f{p) , (D-30) 

where A is the Fourier transform of A, and / is a test function in momentum space. 
Since we assume that the electromagnetic potential ^ = ^Ais macroscopic (D.3), the 
integrand in (D.30) vanishes unless g — p is in the low-energy region. More precisely, 
we can say that 

k° \q-p\ ~ CLro > 

where /macro denotes a typical length scale of macroscopic physics. Since the vector 
g — p is in this sense small, it is reasonable to expand the factor U (q) in (D.30) in a 
Taylor series around p. As the operators Ui are characterized via (D.18), we can assume 
that they have similar regularity properties as P. In particular, we may assume that 
the partial derivatives of Ui{p) scale in powers of Ep^, in the sense that there should 
be a constant c ^ Zmacro-^P such that 

/ c N'"' 

\d'^Ui{p)\ < I -— — I for any multi-index n. (D.31) 

\EpJ 

Prom this we conclude that the Taylor expansion of U{q) around p is an expansion in 
powers of (Zmacro-^'p)"^) and thus 

{UAU-'f){q) = I^U{p)A{q-p)U{p)-'f{p) 

-\- (higher orders in (/macro 

Ep)-^). (D.32) 

Using that A{q — p) is a multiple of the identity matrix, the factors U{p) and U{p)~^ 
in (D.32) cancel each other. We conclude that the operators UAU~^ and A coincide 
up to higher order in {ImacioEp)"^ . For the integral kernels in position space, we thus 
have 

{UAU-'^){x,y) = A{x)5'^{x-y) + (higher orders in (/macro^'p)"^)- (D.33) 

We point out that this statement is much stronger than the equality of the operators 
UAU~^ and A on the macroscopic scale that was mentioned at the beginning of this 



240 



D. THE REGULARIZED CAUSAL PERTURBATION THEORY 



paragraph. Namely, (D.33) shows that these operators coincide even microscopically, 
up to a very small error term. Notice that it was essential for the derivation that A 

is a scalar function (for example, (D.33) would in general be false if we replaced A by 
jf). Using (D.33) in each summand of (D.29) and applying (D.22), we conclude that 

AP[^A]{x,y) = iA{x) Pix,y) - iP{x,y) A{y) 

+ (higher orders in (/macro-Ep)"^). (D.34) 

This shows gauge symmetry of the perturbation calculation (D.25). 

It is interesting that, according to (D.34), gauge symmetry holds only up to an error 
term. This is unproblematic as long as the length scales of macroscopic physics are 
large compared to the Planck length. But (D.34) indicates that the regularized causal 
perturbation theory fails when energy or momentum of the perturbation B are of the 
order of the Planck energy. In this case, the distinction between the "macroscopic" and 
"microscopic" length scales, on which our constructions relied from the very beginning 
(cf. (D.3)), can no longer be made, and it becomes impossible to introduce a causal 
and gauge invariant perturbation theory. 

We conclude the discussion of the regularized causal perturbation expansion by 
pointing out that our construction was based on condition (D.17), which is only a 
necessary condition for causality. Hence the causality of (D.25) has not yet been 
proved. We shall now perform the light-cone expansion of (D.25). This will show 
explicitly that the light-cone expansion involves, to leading orders in (Zmacro-^'p)"^ 
and {lEp)~^, no unbounded line integrals, thereby establishing causaUty in the sense 
of §2.2. 

In the remainder of this appendix, we will analyze the regularized causal pertur- 
bation calculation (D.25) near the light cone. Our method is to first perform the 
light-cone expansion of AQ, and then to transform the resulting formulas according to 
(D.24) to finally obtain the light-cone expansion of AP. In preparation, we describe 
how a decomposition into Dirac eigenstates can be used for an analysis of the operator 
Q near the hght cone. A short computation using (D.21, D.28) yields that Q can be 
represented in the form 

/oo 
di2 w^{k) ein) (^ + ii) dik^ - n^) e(-fc°) (D.35) 
-oo 

with the real-valued distribution 

w^{k) = {cl){k)+tiX{k)) f{k) and kik) = {-\J\k\^ + n^,k). (D.36) 

This representation can be understood as follows. According to (D.9), the distributions 
e{ii) (j^ + jj) S{k'^ — //^) in the integrand of (D.35) are the spectral projectors of the 
free Dirac operator in momentum space. The factor Q{—k^) projects out all states 
on the upper mass cone, and the function w^j_{k) multiplies the states on the lower 

mass shell k = + /W^, k) with a scalar weight factor. In this sense, (D.35) can 

be regarded as the spectral decomposition of the operator Q into Dirac eigenstates. 
Notice that the factor S{k'^ — fi"^) @{—k^) in (D.35) is the Fourier transform of the 
distribution Ta, (2.5.40). Exactly as described for the scalar component in §4.3, we are 
here interested only in the regularization effects for large energy or momentum and 
may thus disregard the logarithmic mass problem (see §2.5 for details). Therefore, we 
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"regularize" Ta according to (2.5.42) and consider instead of (D.35) the operator 

/oo 
d/xe(/x) w^ik)il^ + ^^)Tp{k), 
-oo 

where Ta^^{k) is the Fourier transform of (2.5.42). We expand the distribution T^j^ in 
a power series in 

/oo oo w 

d/x e(/x) w^{k) {^ + ^,)Y,-, ik) 
n=0 

with r(") according to (2.5.43). Commuting the integral and the sum, we obtain 

oo 

Q-g(^) = 327r3 (5[nl(^) ^ + hnp)) T^'^^ik) (D-37) 

n=0 

with 

9[n]ik) 
h[n]{k) 

The representation (D.37) is very useful because it reveals the behavior of the operator 
Q near the light cone. To see this, wc consider the Fourier transform of (D.37) in light- 
cone coordinates (s, I, X2, xa). For the Fourier transform of the factor T^"'\k), we have 
the representation (4.5.2). This representation can immediately be extended to the 
Fourier transform of )j( T^'^\k) by acting on (4.5.2) with the differential operator i^; 
more precisely in light-cone coordinates y — a; = (s, /, X2, X3), 




(D.40) 



In order to treat the factors g^^] ^[n] '^^ (D.37), we note that the Fourier transform of 
(D.37) can be computed similar as described in §4.3 by integrating out the transversal 
momenta according to (4.3.4) and analyzing the remaining two-dimensional Fourier 
integral (4.3.7) with the integration-by-parts method (4.3.14). If this is done, the 
functions g^^] ^[n] appear in the integrand of (4.3.7). Our regularity assumption 
on the fermionic projector of the vacuum (see §4.3 and §4.4) imply that g^^j and are 
smooth functions, whose partial derivatives scale in powers of Ep^. Hence all derivative 
terms of the functions g^^] and h[n] which arise in the integration-by-parts procedure 
(4.3.14) are of higher order in {lEp)~^. Taking into account only the leading order 
in {lEp)~^, we thus obtain a representation of the fermionic projector of the vacuum 
involving only and at the boundary v = a„. Comparing this representation 
with (4.5.2) and (D.40), wc conclude that the Fourier transform of (D.37) is obtained, 
to leading order in (lEp)^^, simply by inserting the functions g^^] and into the 
integrands of (4.5.2) and (D.40), evaluated along the line k = {kx = 2u, ky = Q,kz = 0). 



1 /■ ^ 
= ^ J_ dfie{,,)w^{k)fi-'- (D.38) 

/oo 
dii e{tx) w^{k) i^^'+K (D.39) 
-00 



327r3 
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Thus 



reg 

e~*"* hin] (u) du 



+ (higher orders in {lEpy^), (D.41) 



reg . / T \ ^^S-\ 



where and are the functions (D.38, D.39) with k = {—2u,0,0). 

The decomposition of the operator Q into Dirac eigenstates (D.35) is also useful 
for analyzing its perturbation AQ. 

Lemma D.l. Let B{x) G C^(M^) n L^(M^) he a matrix potential which decays so 
fast at infinity that the functions XiB{x) and XiXjB{x) are also . Then the light- 
cone expansion of the operator /S.Q[B\, (D.23), is obtained by regularizing the light- 
cone expansion of the Dirac sea to first order in the external potential as follows. A 
summand of the light- cone expansion of the Dirac sea which is proportional to , 

(iterated line integrals in bosonic potentials and fields) T^'^\s,l) , 

must be replaced by 

(—1) (iterated line integrals in bosonic potentials and fields) 

J-oo yW"/ [ 9[E] forp even 

+ (rapid decay in I) + (higher orders in (lEp)^^, {ImacroEp)^^ )■ (D.42) 

A contribution ~ mP which contains a factor {y — x)j^^ , 

mF (iterated line integrals in bosonic potentials and fields) {y — x)jj^ r*-"^(s,Z) , 

is to be replaced by 

(—1) (iterated line integrals in bosonic potentials and fields) 



/oo 
du 
-oo 



I \reg / 1 ^ 

21-,' 1-) +a«7'bsT 



h,p-i-, forp odd 
^ 2 J -|- (contributions ^ ^ , J ) 

I g^E] for p even 

+ (rapid decay in I) -\- (higher orders in {lEp) ^, {Imacro 

Ep)-^) . (D.43) 

In these formulas, g\^n] (^''^d h^^] '^f^ ^h,e functions (D.38, D.39) with k = (—2^,0,0). 
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Proof. By substituting (D.6) and (D.35) into (D.23), we obtain the following 
representation for AQ in momentum space, 



/ pp ^ q q 

q q PF 

= £J,,e{^Ji){li + 'L+,,)B,{^-'L+,,) 

X ^ + f) + |) - Mk-l) T,.{k - |)) . (D.44) 

Using the methods developed in [F5], we now perform the light-cone expansion in 
momentum space and then transform back to position space. Since we are here in- 
terested in the regularization effects for large energy or momentum, we may disregard 
the logarithmic mass problem and work on the level of the formal light-cone expansion 
of [F5, Section 3] (our constructions could be made rigorous using the resummation 
method of [F5, Section 4]). As in [F5, Section 3], we expand the distributions r^2 in 
a Taylor series in q and rewrite the resulting ^-derivatives as derivatives with respect 
to /x^. This gives 

°° / 2 \ 

T^.{k±l) = c,r i±kqy (^) Ti;r\k) (D.45) 

with combinatorial factors Cj> whose detailed form is not needed in what follows. Next, 
we expand (D.45) in a Taylor series in /^^ and obtain 

oo / 2 \ 

T^2{k±\) = Cn,.M'"(±A;9)M^J T^-+^+^){k) (D.46) 

n,j,r=Q ^ ' 

with new combinatorial factors c„j>- We substitute the expansions (D.46) into (D.44) 
and write the even and odd terms in kq together, 

_ ..2n /7._\i I 'i \ rn(n+i+r^fi.\ I ... til , ?^ 

2kq 



X 



oo / 2 \ ' / -> -> 

cnjr t^''' (kqy (^] T(-+^+^){k)(w^{k+p-w,{k~^) 
=0, j even 

E 'njrfi''- {kqy (^^J T(-+^+^\k) (w,{k + ^)+W,{k-^)^ 



n,j,r=0, j odd 

(D.47) 

We first consider the contributions to (D.47) for even j. These terms contain the 
factor {w^{k + f ) — w^{k — |)). If the distribution t/;^ were a smooth function and its 
derivatives had the natural scaling behavior in powers of the Planck length, we could 
immediately conclude that \Wfj,{k + |) — 'w^{k — |)| ~ \q[ \dw^\ ~ {ImacmEp)^^ , and 
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9{k + l)-9{k + l] 



(D.48) 



thus all the terms for even j would be negligible. Unfortunately, the situation is more 
difficult because W/j, is in general not a smooth function (cf. (D.36)), and we obtain 
the desired regularity in k only after the yU-intcgration has been carried out. This 
makes it necessary to use the following argument. Consider one summand in (D.47) 
for even j. After carrying out the //-integration, this summand yields a finite number 
of contributions to AQ{k + |, A; — |) of the following form, 

kq 

where each symbol "• • • " stands for a possible factor ^ or ^, and where 5 is a scalar func- 
tion, which coincides with one of the functions g^^] ^[n] (^^^ (D.38) and (D.39)). As 
already mentioned after (D.40), our regularity assumptions on the fermionic projector 
of the vacuum imply that the functions ^j^j and , and thus also g, are smooth, and 
that their derivatives scale in powers of the Planck length. Applying the fundamental 
theorem of calculus, we rewrite the square bracket in (D.48) as a line integral, 

(Z?.48) = J_^d\ — {kqy i^^j • • • • • • r("+^+^)(A;) {qV)g{k + Xq) . (D.49) 

We now transform (D.49) to position space. Our regularity assumptions on B mean 
in momentum space that B{q) G fl L°°. Using furthermore the regularity of Vg, we 
can carry out the g-integration in the Fourier integral. Carrying out also the integral 
over A, we end up with a contribution to AQ(x,y) of the form 

/ (0 T^^'^'^'Hk) F{k, x + y) e-'^^-~y^ (D.50) 

with a (matrix- valued) function F which is differentiable in k and of the order (Zmacro-^'p) 
In the low-energy region, the function g in (D.48) is constant and thus F is homoge- 
neous in k of degree at most j + 1. After transforming to light-cone coordinates, this 
implies that (D.50) is close to the light cone dominated by the fermionic projector of 
the vacuum, in the sense that in light-cone coordinates, |(£).50)| < const(^) \P{s,l)\. 
We conclude that all summands in (D.47) for even j are of higher order in (/macro-^p)"^- 
It remains to consider the summands in (D.47) for odd j. In this case, one factor 
kq cancels the principal value, and we obtain 

AQ[B](A; + |, = -|"d/.6(/x)(^ + | + M)H,(^-|+/x) 

X Y: C^^r f^'- ikqfU\) T(-^'^^'^'\k)Uik+l)+w,{k-l)) 

n,j,r=0 ^ ^ ^ ^ 

-|- (higher orders in (/macro 

Ep)-^) (D.51) 

with some combinatorial factors Cnjr- This formula has similarities to the light-cone 
expansion of the Dirac sea in momentum space [F5, equation (3.15)]. In [F5, Section 
3], we proceeded by rewriting the factors kq as /c-derivatives acting on T^ K When 
taking the Fourier transform, these fc-derivatives were integrated by parts onto the 
exponential factor exjp{—ik{x — y)) to yield factors (y — x). After collecting and rear- 
ranging all resulting terms, we obtained the line-integrals of the light-cone expansion. 
This method can be applied also to the integrand of (D.51), and we can carry out 
the /Li-integration afterwards. We shall not go through all these constructions steps in 
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detail here, but merely consider what happens in principle. Whenever a A;-derivative 
d]^i acts on the factors Wfj, in the integration-by-parts procedure, we get instead of a 
factor (y — x)i (which is obtained when the ^-derivative acts on the exponential 
exp(— i/c(x — y))) a factor (JjW^. After carrying out the /^-integration, one sees that the 
resulting term is of higher order in {lEp)~^. Thus we can, to leading order in {lEp)~^, 
neglect all derivatives of the factors w^. But then, the intcgration-by-parts procedure 
reduces to the construction in [F5, Section 3], and we thus obtain precisely the line 
integrals of the light-cone expansion [F5]. Furthermore, we can replace the factor 
{w^{k+^) + w^{k — |)) in (D.51) by 2w^{k), because a Taylor expansion of this factor 
around ^ = amounts, again after carrying out the ;U-integration, to an expansion 
in powers of (^macro-E'p)""'^, and it thus suffices to take into account the leading term 
of this expansion. These considerations show that the light-cone expansion of (D.51) 
differs from that in [F5] merely by the additional /^-integration and the factor w^(k). 
Hence the light-cone expansion of (D.51) is obtained from that of the Dirac sea by the 
following replacements, 

^ j d^i e(/x) r(-) {k) e-^'^^-y^ w^{k) 

J^J Jf^ (-2^^) r("+^)(A;) e-^(--^) w,{k) 

(where we used the identity (y — xyT^'^\x,y) = 2d^iT^'"'^^\x,y); see [F5, equa- 
tion (3.5)]). The lemma follows by carrying out the yU-integrals applying (D.38, D.39) 
and by analyzing the behavior near the light cone as explained before (D.41). 

■ 

From this lemma we can deduce the light-cone expansion of the regularized fermionic 
projector. 

Theorem D.2. The light-cone expansion of the regularized causal perturbation 
(D.25) is obtained by regularizing the light-cone expansion of the Dirac sea to first 
order in the external potential as follows. A summand of the light-cone expansion of 
the Dirac sea which is proportional to mP , 

(iterated line integrals in bosonic potentials and fields) T^^\s,l) , 

must be replaced by (4-5.7). A contribution ~ which contains a factor (y — x)j^^ , 

mP (iterated line integrals in bosonic potentials and fields) (y — x)j^^ T^^\s,l) , 

is to be replaced by (4-5.9). In these formulas, g, h, a and b are the regularization 
functions introduced in ^4-3 and ^4-4 (see (4-4-^'^> 4-3-25, 4-5-10, 4-4-19))- 

Proof. As mentioned at the beginning of this appendix, we assume here that the 
vector component is null on the light cone (4.4.21). Let us consider what this condition 
tells us about the operators Ui. According to (D.19), the operators [7/ are trivial in 
the low-energy region. Conversely, for large energy or momentum, (4.4.21) yields that 
the vector field v{k) is parallel to k, up to a perturbation of the order ffghear- Hence 
we can assume that the transformation (D.22) is a small perturbation of the identity, 
in the sense that 

ci\Ui{k)-l\ ~ e^hear for ah k. (D.52) 
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We next derive the light-cone expansion of AP by transforming the result of 
Lemma D.l according to (D.24). Since the transformation (D.24) is small in the 
sense of (D.52), it leaves the iterated line integrals in (D.42) and (D.43) unchanged 
to leading order in eghcar- Hence it suffices to consider the transformation of the u- 
integrals in (D.42) and (D.43). The n-integral in (D.42) is as a homogeneous scalar 
operator invariant under the unitary transformations. In the tt-integral in (D.43), on 
the other hand, only the Dirac matrices 7^ and 7* are modified. More precisely, we 
have to leading order in eshean 

^ q {Ua^U^^)(u, v = ay) = 7* + V + (contributions ~ 7^> 7^) 



1=1 

L 



^ ci {Un^Ui ^)(u, v = au) = 7' + 7* + (contributions ~ 7^> 7^) 



=1 



with suitable regularization functions bg and 6; which are small in the following sense, 

9 ^shear • 

Notice that in the high-energy region u ~ Ep, the contribution 7' in the integrand 
of (D.43) is smaller than the contribution ~ 7* by a relative factor of {lEp)~^. Hence 
we can neglect 62, whereas 61 must be taken into account. We conclude that the 
transformation (D.24) of the contributions (D.42) and (D.43) is carried out simply by 
the replacement 

+ ,0.53) 

It remains to derive relations between the regularization functions /?.[„], and bs, 
which appear in the transformed contributions (D.42) and (D.43), and the regulariza- 
tion functions g, h, a, and b in (4.5.7) and (4.5.9). For this, we apply the transformation 
(D.22) to Q'^^^, (D.41). Exactly as described above, this transformation reduces to the 
replacement (D.53), and we obtain the following expansion of the fermionic projector 
near the light cone, 

°° 1 foo / 1 \ reg 

P-^s, = - E i-'ir-' I /.[„](u) du 

°° 1 POO 

n=0 

/ 1 \ reg / , X reg 

-(^-1)^\^J + ill' bin 



X 



du 



+ (higher orders in Sshcar, (lEp) 



Comparing this result with the formulas for the fermionic projector derived in §4.3 and 
§4.4 (see (4.3.27, 4.3.28) and (4.4.6, 4.4.7)), one gets the following identities between 
the regularization functions, 

9in]{u) = g{u)a{uy- , ^[n](") = h{u)a{u)'', bi{u) = b{u) . g 
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We finally explain in which sense the regularized causal perturbation theory is 
unique. In order to ensure regularity of the perturbation theory in the limit when the 

momentum q of the external field goes to zero, one must satisfy a causality condition 
similar to (D.17), and to this end one has to work with a modified fermionic projector 
Q. Since we must modify the direction of the vector field v, it is natural to describe 
the transformation from Q to P by linear combinations of unitary transformations 
(D.22). Nevertheless, we remark that one could just as well work with a different 
or more general transformation Q ^ P. The reason is that the particular form of 
this transformation enters only in the proof of Theorem D.2, and we use merely that 
this transformation is close to the identity, in the sense similar to (D.52). Hence the 
restriction to transformations of type (D.22) is no loss in generality. Furthermore, we 
point out that the gauge symmetry (D.34) uniquely determines the precise form of 
how the potential B enters into the perturbation calculation (e.g. one may not replace 
B in (D.25) by Uj^-^BUi). We conclude that our construction of the regularized causal 
perturbation theory is canonical up to the freedom in choosing the coefficients ci{k) and 
the unitary transformations Ui{k). By assuming that the unitary transformations are 
regular (D.31) and small (D.52), the arbitrariness in choosing {ci,Ui) was constrained 
so much that it has no influence on the regularization of the light-cone expansion. 
Indeed, the c; and Ui do not enter the statement of Theorem D.2. Thus we can say 
that the regularized causal perturbation expansion is unique up to contributions of 
higher order in {lEp)-'^, (ImacToEp)"^ and £shear- 
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Linear Independence of the Basic Fractions 



In this appendix we consider simple fractions of degree L > 2 of the form 




(E.l) 




with integer parameters ct, P > 1, j,S > which satisfy the additional conditions 



We prove the following theorem which makes precise that the only relations between 
the simple fractions are given by the integration-by-parts rules. 

Theorem E.l. Assume that a linear combination of simple fractions (E.1-E.3) 
vanishes when evaluated weakly near the light cone (4-5.29) to leading order in {lEp)~^ 
and (ImacroEp)^^ , for any choice ofrj and the regularization functions. Then the linear 
combination is trivial after suitably applying the integration-by-parts rules. 

The condition (E.3) ensures that the simple fractions are asymmetric under com- 
plex conjugations. Such an asymmetry is essential for our proof. However, (E.3) could 
easily be weakened or replaced by other asymmetry conditions. Also, (E.2) and the 
fact that the denominator involves only the square indices [0] is mainly a matter of 
convenience. The reason why we are content with (E.1-E.3) is that all EL equations 
in this book can be expressed in terms of simple fractions of this form. 

We point out that the above theorem does not imply that the basic monomials are 
independent in the sense that, by choosing suitable regularization functions, the ba- 
sic regularization parameters can be given arbitrary values. Theorem E.l states that 
there are no identities between the basic fractions, but the basic regularization pa- 
rameters might nevertheless be constrained by inequalities between them (e.g. certain 
regularization parameters might be always positive). Furthermore, we remind that the 
assumptions of positivity of the scalar component and of half occupied surface states 
(see the last paragraph of §4.4) yield relations between the regularization functions 
which might give additional constraints for the regularization parameters. For these 
reasons, one should in applications always verify that the values for the basic regular- 
ization parameters obtained in the effective continuum theory can actually be realized 
by suitable regularization functions. 

In the proof we will work with a class of regularization functions for which the 
Fourier integrals and the weak evaluation integral can be computed explicitly; then we 
will analyze in detail how the resulting formulas depend on the regularization. More 



cj,dj e {-1,0} 
a-7 > (3 -5 > \ . 



(E.2) 
(E.3) 
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precisely, we choose the regularization functions in (4.5.12, 4.5.14) as follows, 

g(u) = u"-^ {l + eu'')e~^ e{u) , h{u) = u^" g{u) (E.4) 
a{u) = u^" , h{u) = (E.5) 

with real parameters e, c, z/, Ep > 0. These regularization functions have all the prop- 

u 

erties required in §2.5 if a 1 and e <C <C 1; note that the factor e gives the 
desired decay on the scale of the Planck energy. Using the decay of the integrand for 
large Re u, we can deform the integration contours to obtain for any p > 0, 



Jo V ^EpJ Jo 

where in the last step we set 

1 



= (E.6) 



IS + 



2Ep 

and used the definition of the gamma function. Here the power z^^ is understood as 
exp(— plog z) with the logarithm defined on the complex plane cut along the negative 
real axis. By analytic continuation we we can extend (E.6) to p in the complex plane 
with the exception of the poles of the gamma function. 



f 

Jo 



uP-^ e ^'^p e-'^' du = , for pGC\ {0,-1,-2,...}. (E.7) 



This Fourier integral is also useful for computing the L^-scalar product of the Fourier 
transform via Plancherel. Namely, under the conditions p,p'>^ we obtain 



poo 1 I" 

Jo 27r J_ 

and computing the integral on the left gives 



(E.8) 



The fact that the integral (E.7) diverges when p tends to a negative integer corre- 
sponds precisely to the logarithmic mass problem as discussed after (4.5.1). A short 
calculation shows that the infrared regularization can be introduced here simply by 
subtracting the pole, i.e. for n G No 

r (u-^^-T^ e" W du = hm _ izll! . (e.9) 

Jo ^ ' P—n\ zP n\ p + n) ^ ' 

Using this formula in (4.5.12, 4.5.14), we obtain 



fe=0,l 
fc=0,l 

where "(IR-reg)" means that we subtract a counterterm as in (E.9). 
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For clarity we disregard the infrared regularization for the moment and consider 

(n) 

only the zeroth order in e. Substituting the obtained formulas for To and their 
complex conjugates into (E.l), we obtain 

^_-^^A-C+p-S _ ai + ozy) • • • r{a - a„ + ou) 1 



(E.l) 



{-il)^ T{a - Ci) • • • r(cr - C^) ^{a-i)a-A+C+»v 

T{a — 5i + ov) ■ ■ ■ r((T — 6/3 + oi>) 1 



(E.12) 



r(c7 - di) • • • T{a - ds) ■^(p~&)a-B+D+*u ' 

where A = Yl^=i ^i' ^ — ^j=i ^i' ^ ~ ^i=i '^i' ~ ^j=i ^-J' -'^^^^ ^'^^ parameter o 

takes into account the lower indices of the corresponding factors ri"^-* or To^^^; more 
precisely for an index [p] and {p} it is equal to 4p and 4p + 2, respectively. The 
parameter • stands for the sum of the parameters o in the same line. Integrating 
over s, 



-i- o ' ' ' o o ' ' ' o 



ds , (E.13) 



(E.14) 



[0] -[0] -[0] -[0] 
using (E.8) and leaving out irrelevant prefactors, we obtain the expression 

r((7 — ai + ov) ■ ■ ■ r((T — Ga + Ol') r((T — 6l + ou) • • • r((T — hf) + ov) 

r((T — ci) • • • T{a — c^) r((7 — di) • • • T{a — ds) 

E^-^ r(A - 1) 

^ r((a - 7)a - A + C + r((/? - - 5 + £» + ' 

where A is the sum of the arguments of the two gamma functions in the denominator of 
the second line. Since the £^p-dependence tells us about A, we are led to a combination 
of gamma functions as considered in the next lemma. Although the statement of the 
lemma is not surprising, the proof is a bit delicate, and we give it in detail. 

Lemma E.2. Consider for given parameters N,M G N quotients of gamma func- 
tions of the form 

T{a - ai + vbi) ■ ■ ■ T{a - aj + vhj) 1 

T{a — ci + udi) ■ ■ ■ T{a — ck + I'dx) T{nia — + vmi) T{n2(J — h + i'm2) 

with integers J,K>0 and aj, bj, Cj,dj,nj, lj,mj G Z, which satisfy the relations 

ni + n2 = N , m > n2 > 1 (E.16) 

mi+m2 = M. (E.17) 

// a linear combination of expressions of the form (E.15-E.17) vanishes for all {a, u) 
in an open set o/M^, then the linear combination is trivial after suitably applying the 
identity 

xr{x) = r(x + 1) . (E.18) 

Proof. Assume that a linear combination of terms of the form (E.15-E.17) vanishes 
for all {o;u) in an open set of R^. By analytic continuation we can assume that the 
linear combination vanishes for a and /x in the whole complex plane with the exception 
of the poles of the gamma functions. 

We first consider the asymptotics for large a and i^. If we fix and choose a large, 
we can approximate the gamma functions with the Stirling formula 

r(x) = V2^a;^+^ + (E.19) 
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to obtain 

log(£.15) = {J-K-ni-n2)a (log(a) - 1) - ^ m log(ni) a + C>(log a) . 

1=1,2 

Terms with a different asymptotics cannot compensate each other in the hnear combi- 
nation and must therefore vanish separately. Thus we can restrict attention to a hnear 
combination with fixed values of the parameters 

J — K — ni — n2 and ^nilog(ni). (E.20) 

1=1,2 

More generally, we can choose u = ea for small fixed £ > 0. Then for large a, 
log(E.15) = I ^(1 + ebi) - ^(1 + edk) - ^ (ni + erm) ) a (log(a) - 1) 

\j=l k=l i=l,2 J 

J K 

+ ^(1 + ehj) log(l + ehj) a - ^(1 + edk) log(l + ed^) a 
j=i k=i 

- ^ (nj + errii) log(ni + errii) a + 0{loga) . 

1=1,2 

Expanding in powers of e, we see that the asymptotics also determines the parameter 

m log(ni) . (E.21) 

1=1,2 

We can assume that the parameters (E.20, E.21) are the same for all summands of our 
linear combination. Combining (E.16) with the right of (E.20), we can compute n\ 
and 77-2. Furthermore, (E.17) and (E.21) uniquely determine mi and m2- 
By iterativcly applying (E.18) we can write each term (E.15) as 

V{(T,u) T{a + ubi)---T{a + vbj) 1 



Q(o-, v) T{(T + vdi) ■ ■ ■ r((T + udx) r(nicr + vmi) T{n2(y + vm2) ' 

where V and Q arc polynomials in a and v (which clearly depend also on all the integer 
parameters). After bringing the summands of the linear combination on a common 
denominator, the numerator must vanish identically. Thus it suffices to consider a sum 
of expressions of the form 

l[ Via + vbi) (E.22) 

1=1 

with new polynomials V and parameters L, 6;. According to the left of (E.20), the 
parameter L is the same for all summands. We denote the number of such summands 
by P. For notational convenience we also write (E.22) in the more general form 

L 

V{a,u) Wv{nia + uai) (E.23) 

1=1 
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with new parameters ni and a;. We evaluate (E.23) at successive points a+p with p = 
0, . . . ,P — 1. Again using (E.18), we obtain the expression 

L / nip-1 \ 

7^(0-+^, I/) Y\ inia + q + uai)\r{nia + i'ai). 

1=1 \ q=0 J 

We now consider the asymptotic regime 

CT > > P. 

Then (E.23) simplifies to 

L 

1=1 

It is convenient to divide by IT^iC'^^'^)"'^ (this is possible in view of the fact that the 
parameters and L in (E.22) are known). This gives 

L r -1 mp 

V{a,u)W 1 + ^ V{nia + vai)i\ + 0{pa-^) + 0{pv-^)). 

Since the parameters and in (E.22) are also known, we can finally divide by the 
factors for I = 1, 2 to obtain the expressions 

L-2 

Fp{a,u) := r{a,u)ll [l + k r{nia + ubi) . 

1=1 

These functions satisfy the simple relations 

Fp{a,u) = Fo{a,iy)G{^y (E.24) 

with 

L-2 

G{X) := llil + X bi] . (E.25) 

1=1 

Let us verify that the function G(A) determines all the parameters bi in (E.22): 
Assume that the functions G and G corresponding to two choices of the parameters 6/ 
are equal. Collecting and counting common factors, their quotient G{X)/G{X) can be 
written as 

i=l 

for some parameters G Z satisfying the conditions 

I 

J2qi = 0. (E.26) 

1=1 

Here the bi are (with a slight abuse of notation) a selection of the parameters bi and bi , 
which are all different from each other, and / denotes the number of such parameters. 
We must show that the powers qi are all zero. To this end we take the logarithm, 

I 

logG-logG = 5^ log (1 + A 6,). 

i=l 
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Expanding in powers of A up to the order J, we obtain the conditions 

7 

^qib\ = for all / = 1, 



i=l 



We write these equations in the matrix form 

= 

with 



(E.27) 



A 



( 1 

hi 



1 

hi 



1 \ 
bi 



An elementary consideration shows that 

det^l = J]^ 



I 



( qih \ 

^2^2 



bi) ^ 



(E.28) 



l<i<j<I 



because the hj are all different. We conclude that the matrix A is invertible and 
thus 5' = 0. Hence all the powers qi vanish whenever hi / 0. In the remaining 
case 6j = the corresponding power qi is zero because of (E.26). 

We now return to (E.24). By assumption the sum of the functions Fp{a, v) vanishes, 

p 

= forp = 0,...,P-l, 

a=l 

where the index {a) labels the summands of the linear combination. Using (E.24) and 
keeping a fixed, we can again write these equations in matrix form (E.27) with 



A{v) 



( 1 


1 


G(i) 




G(2) 



p-i 
(1) 



^(2) 



1 \ 

G(P) 
G 



•{P) 



i^p-i 



4'' 



Suppose that "^(vq) / 0. Then there is e > such that ^'(i^) / for all v G Bs{i'q). 
Computing the determinant of A again using the formula (E.28) we conclude that for 
each ly G B^{uq), at least two of the functions coincide. Since there is only a finite 
number of combinations to choose the indices, there must be two indices (a) / (/3) 
such that the function G^"^ — G^^^ has an infinite number of zeros on Bs{uq)Bs{i'q). 
Due to analyticity, it follows that G^") = G^^\ in contradiction to our above result 
that the functions G^"^ are all different. We conclude that *(fo) = and thus 



Fria^i^) = 



for all a = 1, ... P . 



This means that the terms (E.22) all vanish identically. ■ 

Proof of Theorem E.l. Consider a linear combination Z of simple fractions which 

satisfies the assumptions of the theorem. We regularize according to (E.4, E.5, E.IO, 
E.ll), expand in powers of e and evaluate weakly near the light cone (E.13). This 
gives a series of terms of the general form (E.14) (note that due to (E.2) no infrared 
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regularization is necessary in the denominator, and so the counter terms appear only 
in the numerator). The scahng in Ep distinguishes between the contributions with 
different values of A. For every basic fractions (E.l) we introduce the parameter 

Ar = a- 7 + /3-(5 (E.29) 

and let A^max be the maximum which this parameter attains for the basic fractions 
in £. We restrict attention to those contributions (E.14) where A can be written as 

A = Na-l + mv (E.30) 

with integer parameters l,m. These contributions all come from those simple fractions 
for which the parameter (E.29) is equal to A'^. Furthermore, since the counter terms 

in (E.IO, E.ll) involve no factor z~'^ (see (E.9)), they contribute only for A = naH 

with n < N^aax and thus do not show up in our analysis. By considering the con- 
tributions (E.14) with A of the form (E.30), wc will show that all simple fractions 
with N = A^max vanish after suitably applying the integration-by-parts rules. Then 
the corresponding counter terms also drop out, because the infrared regularization 
is compatible with the integration- by-parts rules. Hence these simple fractions com- 
pletely drop out of £, and we can proceed inductively to the analysis of the simple 
fractions with < A^max- This argument allows us to completely ignore the infrared 
regularization in what follows. 

Since the summands with different value of A scale differently in Ep, wc can as- 
sume that the parameter A is the same for all simple fractions in Dividing (E.14) 
through Ep~^ T{X — 1) and using (E.3), we obtain terms which are precisely of the 
form as considered in Lemma E.2. Therefore, our representation of as a linear com- 
bination of quotients of gamma functions (E.14) is unique up to applying (E.18). We 
will consider this arbitrariness later and for the moment consider a fixed choice of 
summands of the form (E.14). 

Our goal is to get a one-to-one connection between quotients of gamma functions of 
the form (E.14) and our original simple fractions. Unfortunately, to zeroth order in e 
one cannot reconstruct the simple fraction from the expression (E.14), because (E.14) 
is symmetric in the parameters aj and bk, and thus it is impossible to tell which of these 

parameters came from a factor tP or T^. This is the reason why we need to consider 
the higher orders in e as well. Out of the many terms of the general form (E.13) we 
select a few terms according to the following rules, which we apply one after the other: 

(i) No factors T{a + v) or T{a -|- 1 + i^) appear. 

(ii) For the factor V{{q. — — n + mv) the parameter m is maximal. 

(iii) The number of factors r(a" — n + (2m + l)i/) in the numerator is minimal. 

By "maximal" (and similarly "minimal" ) we mean that there is no summand for which 
the corresponding parameter is larger (no matter how all other parameters look like). 
Note that the factor r((Q; — j)a — n + mv) in (E.14) is uniquely determined because 
of (E.3). 

Assuming that £ is non-trivial, the above procedure gives us at least one term of 
the form (E.13) (note that the zero-order term in e clearly satisfies (i)). The point is 
that we can uniquely construct from this term a corresponding simple fraction from £, 
as the following consideration shows. According to (i), the factors Tj^j and 
are taken into account only to lowest order in e, because otherwise a factor of the 
form T{a — n + v), n = —1,0, would appear. In particular, one sees from (E.2) that we 
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do not get £-terms of the denominator. Hence the higher orders in e are obtained simply 
by expanding the numerator in (E.l) in powers of e. With the rule (ii) we arranged 
that all factors T^"^ with n > or o 7^ [0] are taken into account linearly in e. On the 

other hand, (iii) ensures that none of the factors tJ""^ is taken into account linearly 
in e. Therefore, all gamma functions in the numerator whose argument contains an 

in] 

odd number times v belong to a factor To . Conversely, the gamma functions of the 
form r((7 — n + 2mv) belong to a factor To"'\ at least when ra > or o ^ [0]. In 
this way, the gamma functions determine the simple fraction up to factors of T^^^\ 



T^Q^ and T^^j '^'^ , T^^^ . But the factors Tj^j and T^^^ can easily be determined from 
the argument of the factor r((a — 7)0- — n + mu), because a — 7 gives us how many 
factors T^^T^ we muse use, whereas n tells us about how many factors ^pj^^ we must use. 
Since A is known, we also know the arguments of the factor T{{f3 — 5)a — B + D + mv) 

in (E.13), and this determines in turn the factors T^j and ^[Qp 

We conclude that the above construction allows us to determine one summand 
of £. Subtracting this summand, we can proceed iteratively to determine all other 
summands of £. This construction is unique up to the transformation of the gamma 
functions with (E.18). 

We conclude the proof by establishing a one-to-one correspondence between the 
transformation (E.18) of the gamma functions and the integration- by-parts rule for the 
simple fraction. To every simple fraction (E.l) we can associate a contribution of the 
form (E.14) which satisfies the rules (i)-(iii) with the following symbolic replacements. 



^[oT'^ ' ^[0]'^ r(a + 1) , r(°) , r(°) Via) 



Tjjf) T{a-n + {4p + iy) 

tI^^ r{a-n + Api^) 



(if n > or p > 0) 



T^^l r(a - n + (4p + 3)i.) , T^^l ^ T{a-n+ (4p + 2)^/) . 

These replacement rules determine the first line in (E.14), whereas the arguments of 

the gamma functions in the second line are obtained as explained above by adding the 
arguments of the gamma functions in the first line. If (E.18) is applied to the gamma 
functions in the first line of (E.14), 

r(cr -n + ou) — > (cr - (n + 1) + ou) T{a - {n + 1) + ou) , 

we take this into account with the following symbolic transformation inside the simple 
fraction. 

Here V is the derivation as introduced in (4.5.34). Using the Leibniz rule this corre- 
spondence can be extended to composite expressions; for example, 

_ / 1 \ vri"^ 



(a — n + oa) T{a ~ n + oa) {a — n + oa) 

r((T — n + 0(t)^ r((T — n + oa) 
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< — > (2(7 — (ni + n2 - 2) + (oi + 02)1/) T{a - rii + oiu) T{a - ni + oiu) 

and similarly for other composite expressions. It remains to consider the transforma- 
tions of the gamma functions in the second line of (E.14). Since A is fixed, wc can only 
increment the argument of one gamma function in the denominator and at the same 
time decrement the argument of the other, for example 

(a - 7)0- — (n + 1) + •lU 
r((a -j)a-n + •lu) T{{P -d)a-m + •2^) 

{(3 — S)a — m + •2v) 

' r((a - 7)c7 - (n + 1) + •ii^) T{{f3 - 5)a - (m - 1) + •21^) ' 
This transformation can be related to the integration by parts rule 

V(...)(...)d. - / i...)Vi---)ds, 

'OO J —00 

where (• • • ) and (• • • ) stand for simple fractions composed of Ti"-* and ri"\ respec- 
tively. As is easily verified, these replacements rules are all compatible with each other 
and with the Leibniz rule. They allow us to identify the transformation (E.18) with 
the integration-by-parts rules. ■ 



APPENDIX F 



The Commutator [P, Q] 



The Euler-Lagrange equations corresponding to our variational principles involve 
the commutator [P, Q] (5.2.9), where Q is a composite expression in the fermionic 
projector. In Chapter 4 we developed a method with which composite expressions in 
the fermionic projector can be evaluated weakly near the light cone. In this appendix 
we shall describe in detail how these methods can be used to evaluate the commuta- 
tor [P, Q] in the continuum limit. We begin by collecting a few formulas from §4.5 and 
bring them into a form convenient for what follows. The kernel Q{x, y) can be written 
as a linear combination of terms of the form (4.5.22) 



rp{ci) rp{Cj) rp(dl) rp{ds) 

where / is a smooth function composed of the bosonic fields and fermionic wave func- 
tions. Here the factors ri"^-* and ri''*^ are the regularized distributions of the light-cone 
expansion. The quotient of monomials is called a simple fraction, and its degree L is de- 
fined by (4.5.27). If L > 1, the monomial becomes singular on the light cone when the 
regularization is removed by letting Ep — > oo. In light-cone coordinates {s,l,X2,xs), 
this singular behavior on the light cone is quantified by a weak integral over s for fixed 
X2, X3, and I S> Ep^. More precisely (4.5.29), 



J — ( 



o o o ' ' ' ^ o 

-,(ci) rp(C'y) rpidl) rp(ds) (U)^ 

+ (higher orders in (lEp)-^ and (/macro^p)"^) , (F.2) 



rp{ci) (c-,) rp{dl) rpid''' < ■> ' 

■J- o ' ' ' ^ o o ***-^o 



where g is an integer, Creg is the so-called regularization parameter, and 77 is a test 
function, which must be macroscopic in the sense that its derivatives scale in powers 
of or /macro- The asymptotic formula (F.2) applies on the upper light cone s = 0, 
but by taking the adjoint and using that Q is Hermitian, Q{x,y)* = Q{y,x), it is 
immediately extended to the lower light cone. Furthermore, we can integrate (F.2) over 
/, X2, and X3, provided that I 3> Ep^ . In polar coordinates (y — x) = {t, r, O = (1?, ip)), 
we thus have 



/oo poo p 

dt / r'^dr dnri{t,r,n) f{x,y) 
-00 Jrn Js^ 



rp{ai) rpiaa) rji{bl) rpibp) 

-to ' ' ' J-0 J-o ■■■-to 



= log^(Ep) E^p-' [ fdtf dn (r?/)|,=|,| 

JM\[-ro,ro] Js^ m 

+ (higher orders in {rEp)~^ and (/macro-Ep)"^) , 
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valid for every ro ^ Ep^. Next we expand the function / in a Taylor series around 
(y - x) = e = 0, 

f{x,y) = Y.fj{x)e (F.3) 

J 

with J a multi-index, and write the Taylor coefficients together with the regularization 
parameter. Collecting all contributions, we obtain for Q the weak evaluation formula 



/OO /"OO 
dt I r'^dr dQ r]{t,r,n) Q{x,y) 
-OO Jrn Js^ 



'ro JS^ 

LraavL 9max „ ~ 

+ (higher orders in ((r + \t\) Ep)-^ and (Zmacro^p)"^) + o{Ep) (F.4) 

with suitable functions hj{0,), which depend on L and g. The integrand on the right 
side of (F.4) is evaluated on the light cone r = \t\. The maximal degree of the mono- 
mials -Lmax as well as ^max are clearly finite parameters. Notice that the monomials of 
degree L < 2 are omitted in (F.4); this is justified as follows. For L < 2, the integral 
(F.2) diverges at most logarithmically as Ep oo, and furthermore has a pole in I of 
order at most one. Thus the corresponding contribution to (F.4) is at most logarith- 
mically divergent as Ep oo, with bounds uniform in ro. This is what we mean by 
o{Ep). 

We point out that the asymptotic expansion near the light cone (F.4) does not 
give any information on the behavior of Q{x,y) near the origin, i.e. when x and y are 
so close that r,\t\ ~ Ep^. Namely, due to the restriction ro S> Ep^, the region near 
the origin is excluded from the integration domain. Also, near the origin the terms of 
higher order in ((r + \t\)Ep)~^, which are left out in (F.4), cannot be neglected. As 
explained in detail in Appendix D, the reason for this limitation is that near the origin, 
Q depends essentially on the detailed form of the fermionic projector on the Planck 
scale and thus remains undetermined within the method of variable regularization. 

Our aim is to evaluate the commutator [P, Q] using the expansion (F.4). The main 
difficulty is that products of the operators P and Q, like for example 

{QP){x,y) = J d''zQix,z)P{z,y), (F.5) 

involve Q{x, z) near the origin x = z, where (F.4) does not apply. In order to explain 
our strategy for dealing with this so-called problem at the origin, we briefiy discuss a 
simple one-dimensional example. Assume that we are given a function f{x), a:; G M, 
and a positive integer n such that for all xq S> Ep^ and test functions ry, 



f{x) ri{x) dx 

iMXl—xojXo] 

= [ dx + (higher orders in (xE^p)"^). (F.6) 

In analogy to (F.4), this formula does not give any information on the behavior of 
f{x) near the origin x = 0. Thus there are many different functions satisfying (F.6), 
a typical example is 
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The question is if (F.6) is useful for analyzing the weak integral 

/oo 
f{x) r}{x) dx . (F.8) 
-oo 

The answer to this question depends very much on the properties of rj. If rj is an 
arbitrary test function with compact support, we can restrict attention to test functions 
with support away from the origin, supp C M \ [—xo,xo]. Then (F.6) applies, and 
we find that f{x) ~ x~'". Thus by evaluating (F.8) for suitable test functions, we 
can find out that, as long as \x\ 3> Ep^, f{x) behaves like the function which 
has a pole of order n at the origin. We refer to this statement as we can detect the 
pole of f by testing with rj. Unfortunately, the situation becomes more difficult if we 
assume that rj belongs to a more restricted class of functions. Assume for example 
that r]{x) is rational, goes to zero at infinity and has all its poles in the upper half 
plane {Im x > 0}. Then for / as in (F.7), the integral (F.8) can be closed to a contour 
integral in the lower complex plane, and we get zero, independent of 77. This shows 
that when testing only with rational functions with poles in the upper half plane, the 
formula (F.6) is of no use, and we cannot detect the pole of /. Indeed, the problem 
in (F.5) can be understood in a similar way. If we apply the operator product QP to a 
test function rj and write the result as {QP)r] = Q{Pr]), the problem of making sense 
out of the integral (F.5) can be restated by saying that Q may be tested only with the 
functions Pr]. In other words, the test functions must lie in the image of P, i.e. they 
must be negative-energy solutions of the Dirac equation. Thus the question is if by 
evaluating only with such special functions, can we nevertheless detect the poles of Q, 
and if yes, how can this be done? Once these questions are settled, we can compute 
the operator products PQ, QP and take their difference. 

For clarity we begin the analysis with the simplified situation where both P and 
Q are homogeneous, i.e. 

Pix,y) = P{y-x), Qix,y) = Q{y - x) . (F.9) 

Under this assumption, the operators P and Q are diagonal in momentum space, 

nx,y) = J ^ P{k) e-''^^-y\ Q{x,y) = J ^ Q(k) e'^'^^-y^ 

and their products can be taken "pointwise" in k, i.e. in the example (F.5), 

{QP){x,y) = J ^ Q{k) P{k) e-^'^^-yK (F.IO) 

Due to this simplification, it is preferable to work in momentum space. Now the 
problem at the origin becomes apparent in the Fourier integral 

Q{k) = jd^iQiOe-'^^ ii^y-x), (F.ll) 

where we must integrate over a neighborhood of ^ = 0. In order to handle this problem, 
we must carefully keep track of how the unknown behavior of Q near the origin effects 
the Fourier integral: If we consider the integral in position space 

/ d'imQ{^) (F.12) 

with a smooth, macroscopic function 77, we can make the unknown contribution near 
the origin to the integral small by assuming that rj goes to zero sufficiently fast near 
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^ = 0. Thus we expect that (F.12) is well-defined provided that the partial derivatives 
of at ^ = vanish up to the order n, 

d^r]{0) = for all / with \I\ < n, (F.13) 

where n is a sufficiently large parameter (which we shall specify below) . In momentum 
space, the conditions (F.13) take the form 



d^'k rj{k) = , |i"| < n. 

This means that the Fourier transform of Q is well-defined by (F.ll), as long as it is 
evaluated weakly only with test functions fi{k) which are orthogonal to the polynomials 
k^ . Equivalently, we can say that Q{k) is defined only up to the polynomials k^ , 

Q{k) = J Q(0 e-'''^ mod V'ik) , (F.14) 

where V"'{k) denotes the polynomials in k of degree at most n. Let us compute the 
Fourier integral (F.14) using the expansion (F.4), and at the same time determine the 
parameter n. Our method is to consider (F.4) for 77 = exp(— ■iA;^) and to choose n so 
large that we can take the limit tq \ to obtain the Fourier integral (F.14). For each 
summand in (F.4), the resulting t-integral is of the form 

lim / dtt^-^+l-^le-^'^^mod 7'"(a;). 

'•oXO 7lR\[-ro,r-o] 

Here in the integrand one may distinguish between the two regions \t\ > lo^^-. where 
the factor t2-L+|J| jg j-gg^i^r and e"*"^* is oscillating, and \t\ < uj^'^, where the pole 
of t^-i+l./! niust be taken into account and the exponential is well-approximated by a 
Taylor polynomial. Since we are interested in the scaling behavior of the integral over 
the pole, it suffices to consider the region \t\ < OJ^^, and calculating modulo V^{uj), 
the leading contribution to the integral is 

lim / ^^^2-L+|J| (7^^)""^' . (F.15) 



ro\oJ[^^-i^^-i] \ [-ro,ro] (n + 1)\ 

If n > L — I J| — 3, the integrand is bounded near t = 0. In the case n = L — \ J\ — 4, 
the limit in (F.15) may be defined as a principal value, whereas for n < L — | J| — 4, 
(F.15) is ill-defined. Thus we need to assume that n > L — \ J\ — A. Moreover, we must 
ensure that the terms of higher order in ((r -|- \t\)Ep)~^, which are omitted in (F.4), 
are negligible in the Fourier integral. Since these terms are regularized on the Planck 
scale, the scaling behavior of these higher order terms is in analogy to (F.15) given by 
the integrals 



I 



-11 \ hE-\E-^] itEpY 



^3-L+|J|+n 

dt — — r-- , n> 1. 



A simple calculation shows that these integrals are negligible compared to (F.15) if 
and only if n > L — | J| — 3, and under this assumption, they are of higher order in 
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io/Ep. We conclude that the Fourier transform of Q has the expansion 

L=2 g=0 J 

/oo f 
dt t-^+\J\+^ / hj{i) i'^ e-^*^«* mod V^-\^\-^{k) 
-oo Js^ 

+ (higher orders in k/Ep and (/macro-E'p)^^) + o{Ep) , (F.16) 

where ^ is the "unit null vector" ^ = (1, O G C M^) and hj{i = (1,0)) = hj{Q). 
Carrying out the t-integration gives the following result. 

Lemma F.l. Suppose that the operator Q is homogeneous. Then its Fourier trans- 
form Q is of the form 

Lmax gmax 

Qik) = J2 Y.log<^iEp)E^-'Y.Q^/ik)modV''-\'\-\k) 

L=2 g=0 J 

+ (higher orders in k/Ep and {Imacro 

Ep)-^) + o{Ep) (F.17) 

with 

Q^fik) = -2m {-iuj)^-\''\-^ 
/ d^ hj{0 X < 



X 



Q{ujU) if\J\<L-2 



{L-\J\- 3)1 ' ^' ' (F.18) 

^ e(a;)^-|-^|-3 6(^+\'^\-^\ki) if\J\>L-2, 

where uj = k^ is the energy and k = k/u (e is again the step function e{x) = 1 for 
a; > and e{x) = —1 otherwise). 

Proof. The t-integral in (F.16) is of the form 



/. 



00 

t-'^e-'^^dt mod V^'^ik) 



00 



with n = L — \ J\ — 2 and A = kS,. For n = 0, we have 

e-'^^dt = 2'k5{\). (F.19) 



/ 

J — < 



The case n < follows by differentiating this equation (— n) times with respect to A. 
In order to treat the case n > 0, we integrate (F.19) n times in the variable A. The 
integration constant is a polynomial in A of degree n — 1 and can thus be omitted. ■ 



Let us briefly discuss the above expansion. The parameters L and g give the scaling 
behavior in the Planck energy. The multi-index J enters at two different points: it 
determines via the factor a;~l"^l the dependence on the energy, and it also influences the 
5^-integral. This integral gives detailed information on the behavior of Q{k) in k, but 
it is independent of |a;|. Integrating over 5^ takes into account the angular dependence 
of the regularization functions and of the macroscopic physical objects and tells about 
how the different angles contribute to Q. In the case | J| > L — 2, the integrand has 
a 5- like singularity localized at fc^ = 0, and so the S'^-integral reduces to integrating 
over the intersection of the hyperplane {£, \ k£, = 0} with the two-sphere t = 1 = r. 
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This intersection is empty for time-like A;, and is a one-sphere for space-hke k. As a 
consequence, Qj^{k) is zero inside the mass cone C and in general non-zero outside, 
without being regular on its boundary {k"^ = 0}. If on the other hand \J\ < L — 2, 
the factor Q(ujkS^) is essential, because without this factor, we would simply have a 
polynomial in k of degree L — | J| — 3, being zero modulo V^'~^'^^~^ {k) . Note that for any 
^, the factor @{u!k$,) vanishes inside the lower mass cone {k^ > 0, u < 0}, whereas it is 
in general non-zero otherwise. This means that Qj^ik) again vanishes in the interior 
of the lower mass cone and is not regular on the mass cone {/c^ = 0}. The singular 
behavior on the lower mass cone is made more explicit in the following lemma. 

Lemma F.2. The operators Qj^{k), (F.18), vanish inside the lower mass cone 
{k^ > 0, oj < 0}. Near the lower mass cone, they have the asymptotic form 



ILU 



L-|J|-3 



X < 



hj{k) F(1 + 0(P)) 
^ (_l)l+|J|-L^(l+|J|-L)^P) if\J\>L-2. 



p(L-|J|-2) 



(F.20) 



Proof. Without loss of generality, we can assume that k points in the tx-direction 
of our Cartesian coordinate system, i.e. k = (1,A, 0, 0) with A > 0. Then k$, = 
1 — A cos I?, and we can write (F.18) in the region a; < as 

/I /'27r 
dcos?9 / d(phj{C)l'' 

Acos^?)^-l-^l-3 



Qj 



X < 



(1 



e(Acosi?-l) if|J|<L-2 



(L-lJl-3)! 

^ (-l)^-|-^|-3 5(2+kl--'^)(l- Acosi?) if|J|>L-2. 

Inside the lower mass cone, the parameter A < 1, and the integrand is identically equal 
to zero. Outside and near the lower mass cone, 1 < A ~ 1, and the integrand vanishes 
unless COST? « 1. Hence to leading order in k^, we may replace /ij(^) ^'^ by its value 
at the coordinate pole = and carry out the (/^-integration. 



-47r^ i {-iui) 



i: 



du X < 



"i-^i-^' hjik) F(i + o(P)) 

(^-^^)'-'''-^(A.-l) 



(L-lJl-3)! 

^ (^_-^^L-\J\-3g{2+\J\-L)(^l_Xu) 

In the case \J\ < L — 2, the remaining integral is of the form 



if |J| < L-2 
if |J| >L-2. 



whereas in the case \ J\ > L — 2, 
^(")(1-Au)dtx 



1) du 



A(n + 1) 



(i-Ar+^e(A-i) 



-) y__„-«(i-A„)<i„ 

N n 

J (A"-ie(A-i)). 



d_ 
dX 
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Finally, we use that A = 1 + 0(P) and P = 1 - = 2(1 - A) + 0(A;'^). 
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According to Lemma F.l and Lemma F.2, all the information on the behavior of Q 
on the light cone contained in (F.4) is in momentum space encoded in a neighborhood 
of the lower mass cone. More precisely, this information can be retrieved as follows. 
Due to the factor a;"!'^' in (F.18) and (F.20), the sum over the multi-index J in (F.17) 
is an expansion in powers of uj~^. Thus by considering Q for large energies, more 
precisely for lo in the range 

Wcro < 1^1 < Ep, (F.21) 

we can make the contributions for large \J\ small (the restriction |u;| <C Ep is clearly 
necessary because the terms of higher order in k/Ep are omitted in (F.17)). Hence in 
this energy range, the series in (F.17) converges fast, and the scaling behavior in Ep 
and as well as the dependence on k given explicitly in (F.20), allow us to determine 
the functions hjip) completely. 

Having computed the Fourier transform of Q, we can now take the product with 
the operator P according to (F.IO). We begin with the simplest case where we take 
for P one massive Dirac sea in the vacuum (2.6.13) with m > 0. In this case, P{k) 
is supported inside the lower mass cone C^. According to Lemma F.2, the operators 
Qj^ vanish identically inside the lower mass cone. Hence the supports of Qj^ and P 
do not intersect, and it follows immediately that 

Q^f p = 0. (F.22) 

This means that after multiplying by P, all the information contained in the expansion 
of Lemma F.l is lost. Wc refer to this difficulty as the problem of disjoint supports. 
Using the notion introduced after (F.8), it is impossible to detect the poles of Q by 
testing with the negative-energy solutions of the free Dirac equation. This situation 
is indeed quite similar to the example (F.6, F.8) for / according to (F.7) and rational 
test functions with support in the upper half plane, in particular since the Fourier 
transform r]{k) = 7?(x) exp(— z/cx) of such a test function is supported in the half 
line A; < 0, and can thus be regarded as a one-dimensional analogue of the negative- 
energy solutions of the free Dirac equation. 

It is an instructive cross-check to see how the problem of disjoint supports comes 
about if instead of analyzing the behavior of Q in position space (F.4) and then trans- 
forming to momentum space, we work exclusively in momentum space. For simplicity 
we give this qualitative argument, which will not enter the subsequent analysis, only 
in the special case of a monomial, i.e. instead of (F.l) for an expression of the form 

f{x, y) ri"^^ • • • ri"''^ T^^'^ ■ ■ ■ ri''^^ . (F.23) 

In this case, we can, instead of taking the product of the factors Ti"^ and Tq^^ in position 
space, also compute their convolution in momentum space. As explained in §4.2, the 
singular behavior of the fermionic projector on the light cone is determined by states 
near the lower mass cone. More precisely, the main contribution to P{x,y) comes 
from states close to the hypersurface TC = {k \ = 0} , which for ^ on the light cone is 
tangential to the mass cone C = {/c^ = 0}, so that the singularity on the light cone can 
be associated to the states in a neighborhood of the straight line HnC. For objects 
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derived from the fermionic projector like the regularized distributions To and 
this qualitative picture applies just as well. The Fourier transforms of the factors ri"-'' 
and T^^^ are supported in the the interior of the lower mass cone. Thus when forming 
their convolution, 

h := Ti-) * ... * fi»^) , ^2 := ^ * ... * fi^^) , (F.24) 

the resulting convolution integrals are all finite integrals over a compact domain (e.g., 
the integrand in (Ti"^-* * ri"^^)(fc) = J d^q To"'^\q) To"'^\k — q) vanishes unless q lies 
in the "diamond" {q^ > 0, g° > 0} n {{q - kf > 0, (g - kf < 0}). Moreover, the 
supports of gi and 52 are again inside the lower mass cone. Exactly as described for 
the fermionic projector in Section C.3, the behavior of gi and ^2 near the lower mass 
cone determines the well known singularities of gi and g2 on the light cone, whereas 
the form of g\ and 52 in the high-energy region away from the mass cone depends 
essentially on the details of the regularization and is thus unknown. Using gi and 52 > 
we can write the Fourier transform of the monomial as 

M^^(fe) := j (fi fi"^) * ... * Ti"''^ fi^^^ * ... * f^^"^ e-'^^ (F.25) 

d'^qgMUQ-k). (F.26) 



In (F.l) the monomial is multiplied by the smooth function /. Thus the corresponding 
contribution to Q is obtained by taking the convolution of / with M^^, 

Q X f*M^^ . (F.27) 

Since / is a macroscopic function, its Fourier transform f[q) is localized in a neigh- 
borhood of the origin, i.e. in the region \q\ ~ ^macro. '^^^ Taylor expansion (F.3) 
corresponds to expanding / in terms of distributions supported at the origin, more 
precisely 

f = J2fj ^it^ fj^^) = i27r)^ f J {idkys\k), (F.28) 
J 

and substituting this expansion into (F.27) yields the expansion of Lemma F.l, 

^ fj*M^9_ (F_29) 

Since the distributions fj{q) are supported at g = 0, the support of Qj^ coincides with 
that of M^9. Hence to discuss the problem of disjoint supports, we must consider M^^ 
as given by the integral in (F.26). Note that this integral differs from a convolution in 
that the argument of §2 has the opposite sign; this accounts for the complex conjugation 
in (F.25). As a consequence, the integration domain is now not compact, and the 
integral is finite only due to the regularization. More precisely, the integration range 
is the intersection of two cones, as shown in Figure F.l in a typical example. We 
have information on the integrand only when both q and q — k are close to the lower 
mass cone, i.e. when q lies in the intersection of the dark shaded regions in Figure 4.1. 
Outside of this so-called intersection region, however, the integrand depends on the 
unknown high-energy behavior of gi or ^2- Since the intersection region becomes large 
when k comes close to the mass cone and does not depend smoothly in k for k on the 
mass cone, the contribution of the intersection region to the integral is localized in a 
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neighborhood of and is not regular on the mass cone {/c^ = 0}. The contribution of 
the high-energy regions to the integral, on the other hand, is regular in k and thus 
well-approximated by a polynomial in k. This qualitative argument illustrates why in 
Lemma F.l, Q is determined only modulo a polynomial, and why the singular behavior 
on the light cone (F.4) is in momentum space encoded near the lower mass cone (see 
Lemma F.2 and the discussion thereafter). An important conclusion from Figure F.l 
is that Q{k) is in general not zero in the interior of the lower mass cone, and even the 
intersection region gives a contribution there. Thus for generic regular izations or for 
simple regularizations like mollifying with a smooth function, the supports of Q and 
P will have a non-empty intersection, and even the singularities of Q on the light cone 
will contribute to the product QP. On the other hand, it seems possible that there 
are special regularizations for which the contributions from the high energy regions 
and the intersection region compensate each other in the integral (F.26) in such a way 
that M^^ indeed vanishes inside the lower mass cone. We call such a regularization an 
optimal regularization. According to the method of variable regularization (see §4.1), 
we want to keep the regularization as general as possible. Therefore, we must allow 
for the possibility that the regularization is optimal, and this leads to the problem of 
disjoint supports. 

The above consideration in momentum space gives a hint on how to resolve the 

problem of disjoint supports. Namely, let us assume for a moment that the macroscopic 
function / has nice decay properties at infinity. Then its Fourier transform / is a 
regular function. As a consequence, the convolution (F.27) mollifies M^^ on the scale 
^macro) ^^'^ the support of Q will be larger than that of M^^. Clearly, /macro is very 
small on the Planck scale, but since the mass shell {A;^ = m^} and the mass cone 
{k^ = 0} come asymptotically close as the energy \k^\ gets large, mollifying even on a 
small scale leads to an overlap of the supports of Q and P. This is an effect which is 
not apparent in the expansion of Lemma F.l because by expanding / in a Taylor series 
around ^ = 0, we did not use the decay properties of / at infinity, and thus we did not 
see the smoothing in momentum space (cf. also (F.28) and (F.29)). More generally, 
the above "mollifying argument" shows that the supports of Q and P should overlap 
if we take into account the macroscopic perturbations of P and Q more carefully. Thus 
in order to solve the problem of disjoint supports, we shall now compute the product 
QP in the case of general interaction, without assuming that P or Q are homogeneous. 
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Our key result will be an expansion of the operator product around the light cone (see 
Theorem F.5 below). 

Let us specify our assumptions on P and Q. For Q wc merely assume that the 
weak evaluation formula (F.4) holds. For P, on the other hand, we work with the 
formulas of the light-cone expansion, which are of the general form 

oo 

P{x,y) = ^ 5(p(a;,i/) r^)(a;,?/) + (smooth contributions) . (F.30) 

p=-i 

Here the gp are smooth functions involving the bosonic potentials and fields, and the 
unspecified smooth contributions are composed of the fermionic wave functions as well 
as the non-causal contributions to the Dirac seas (see (2.3.19, 2.5.45) and Appendix B). 

For clarity, we shall consider the product of Q with each of the summands in (F.30) 
separately, i.e. we will for given p > —1 compute the product 

QR with R{x,y) = g{x,y)T^\x,y) (F.31) 

and a smooth function g. To avoid confusion, we recall that in the case p = —1, T^P^ 
is defined via a distributional derivative; more precisely, we assume in this case that 
g{x, y) has the form g = with smooth functions f^{x, y) and set similar to (2.5.21) 

5(x,2/)T(-i)(x,y) = -2/^(x,y)^r(0)(x,y). (F.32) 

For technical convenience, we assume furthermore that 5 is a Schwartz function, g G 
5(M^ xM^), but this assumption is not essential and could be relaxed by approximation 
(see the discussion after (F.81) below). 

The contributions to Q in (F.4) are supported on the light cone. Thus we can 
write them in the form 

Q{x,y) = h{x,y) Ka=o{x,y) (F.33) 

with 

Ka=o{x,y) = -^S{e)e{e) (F.34) 

and a function h(x,y), which in general will have a pole at the origin x = y. This 
representation is useful because Ka is a solution of the Klein-Gordon equation, namely 
in momentum space 

Ka{k) = 6{k^ - a) e{k°) , aG[0,oo). (F.35) 

In what follows, we will also need the Green's function Sa of the Klein-Gordon equation 
defined by 

PP 1 1 

As is immediately verified with contour integrals, this Green's function is for a > 
causal in the sense that Sa{x,y) vanishes for space-like On the contrary if a < 0, 
Sa{x,y) vanishes for time-like ^. More precisely, the Green's function can be written 
as 

Saix, y) = 5{e) + Q{ae) e{a) Ha{x, y) , (F.37) 
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where Ha is a smooth solution of the Klein- Gordon equation with power expansion 

It is convenient to also introduce the Green's function 

= Sa-Q{a)Ha, (F.39) 

which for all a G M vanishes in time-like directions. As one sees explicitly using (F.38) 
and (F.37), both Ha{x,y) and S'^{x,y) are analytic in a for all a G R. Similarly, a 
short explicit calculation shows that Ka{x,y) is analytic for a E [0, oo). We set 

^•"'-S3(^)"^- 

The following lemma gives the light-cone expansion for an operator product involving 

two factors K^'\ A major difference to Lemma 2.5.2 is that the expansion now contains 
unbounded line integrals; this also requires a different method of the proof. 

Lemma F.3. The operator product K^^^VK^^^ with l,r > and a scalar function 
V E S has the light-cone expansion 

(K« T/ii:W)(x,y) 

-1 oo _ „oo 

+{non-causal contributions, smooth for x ^ y). (F.40) 

We point out that, exactly as in §2.5, we do not study the convergence of the 
infinite scries in (F.40), which are merely a convenient notation for the approximation 
by the partial sums. 

Proof of Lemma F.3. We first consider the operator product KaVKh for a,b > 
in the case when F is a plane wave, V{x) = exp{—iqx). Then in momentum space 
(similar to [F6, eqn (3.9)]), the operator product takes the form 

{KaVK,){p+^,p^^) = Ka{p+l) K,(p-l) ■ (F.41) 

If < 0, we get contributions when either the two upper mass shells of the factors 
Ka and Ki, intersect, or the two lower mass shells. Conversely if > 0, we only get 
cross terms between the upper and lower mass shells. Thus setting 

a + b a — b /-.^ 

u = , V = — — , (F.42) 
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we have 

= -6(g2) 5[{{p - \f -h) + {2pq - 2v)) 6({p - |)2 - b) 



-e{q ) S{2pq — 2v) S[ p — pq + — — b 



= --^^{q )^ip(i-v)s\^p + — -u] . 

Hence we can write our operator product as 



KaVKh = ^Auv (F.43) 



where A^v is the operator 



A. (p + |, P - I) = ^ S{pq -v)e (p2 + ^ - u)^ . (F.44) 

Our strategy is to first derive an expansion for Auy. Then we will differentiate this 
expansion with respect to u and v and take the limits a,b \ to get the desired 
expansion for i^Wy/cW. 

The right side of (F.44) involves a product of the form 6{a) @{P)- This product 
can be transformed into a line integral as follows. Consider for £ > the function 

P) = - &ieP - a^) (eP - a^)"^ . (F.45) 

TT 

This function is zero unless /3 > and a G [—^/eP, As e \ 0, the size of this last 

interval tends to zero, and so a is confined to a smaller and smaller neighborhood of 
the origin. On the other hand, the integral over a stays bounded in this limit; namely, 



i: 



fe{a,(3)da = @{(3) for alU > 0. 
Prom this we conclude that 



/ 

J — ( 



hm /,(«,/?) = S{a)e{P) (F.46) 
with convergence as a distribution. Moreover, the formula 

r-oo pp 

^ dT = TT 6(7) 7-2 (F.47) 

/_oo + 7 

allows one to write (F.45) as a contour integral. Putting these relations together, we 
obtain for A^v, 

^ni; (p + |, ]3 - I) ^'^■^tf'^^^ i lim fs {pq - V, q^ (p^ + ^ - u] 

lijn ^ e f + ^ ~ ^ _ ipq-^f \ ( P^ + T^^ _ {PQ-vf \ ~' 
2tt e\o \eq'^\ I eq^ e^q^ j\ eq^ e^q"^ j 

lim r !^ -, W "r ■ (F.48) 



27r2 e\0j_^ 2 2, 2,9^ (pq - v)' 

eq T +p + — — u — 



4 £g2 
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After shifting the integration variable according to r — r + {pq — v)/{sq^), we can 
identify the integrand with the Green's function (F.36), 

- hn. r ^dr 



/oo 
S,ie,r){p + rq), (F.49) 
-oo 



27r2 

where 2; is the "mass function" 

z{£,t) = u + 2tv + {1- £^) 5^ - ^ . 

If we solve (F.39) for Sa and substitute into (F.49), we can take the limit e \ to 
obtain 

= ^ l°°JS^ + e{z)H,){p + Tq)dT (F.50) 

with 

z = u + 2tv + {t"^ - J) g2 . (F.51) 

The calculation so far was carried out for fixed momentum q of the potential. In 
order to describe the case of general V e S, we must integrate over q. Furthermore, 
we transform to position space by integrating over p (similar to [F5, eqn (3.10)]) and 
obtain 

A,A^,y) = /|^y„)/|?,A„(p + |.p-|)e-<*-.)e-'J(-.), 

where V is the Fourier transform of V. Substituting in (F.50) and pulling out the 
r-integral gives 



/ 



"^'^ {S^ + e(z) H,){p + Tq) e-^vi^-y) 



(2vr)4 

and, after shifting the integration variable p according to p + rg — p, we can carry out 
the Fourier integral, 

/ ('^-^ + ®(^)^-)(^^ + ^«)^"'''^""'^ = e^''"^"-^^('5r + e(z)J/.)(z,y), 

and thus obtain 

In this way, we have transformed the line integral, which appeared in (F.48) as a 
contour integral in momentum space, into an integral along the straight line (^ — 
r)x + (^ + T)y through the space-time points x and y. 
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The operator product K^^'^VK^'^^ is obtained from A^^ by differentiating with re- 
spect to u,v and setting u = Q = v. More precisely, using (F.42) and (F.43), 

(i.'" V K^-K^.y) .^{l.^^)\l-±yi ^„|_„ . ,F.53) 

According to (F.51), the factors S*^, 9(z), and in (F.52) depend imphcitly on 
u and V. Thus when substituting (F.52) into (F.53), we can carry out the partial 
derivatives with the sum. product, and chain rules. Let us first collect the terms for 
which all the derivatives act on the factors or Hz- This gives the contributions to 

i/>.(i + ry(i-.r/|?,V(,)etf)e-.i<-.) 

with z = (r^ — 3 ) • After expanding in powers of z and introducing the new integration 
variable a = r + i , we obtain 



2^2 E y_ da a' (1 - af {a - J ^ e{q^) {-^T V{q) 



X 



If V is supported outside the mass cone, supp V C {q^ < 0}, we can carry out the 
g- integral in (F.54) and obtain precisely the two series in (F.40). Conversely if V is 
supported inside the mass cone {q^ > 0}, we get 



1 °° 1 /■ 

^^q"! 7ir\[0,1] 
1 00 -. „oo 



27r2 

n=0 

This does not coincide with the series in (F.40). Using (F.39) and (F.37), the difference 
can be written as 



1 ^ J poo 

— — 77 — / da a'' (1 — aY (a — a 
27r2 ^n[J_^ 



2\n 

n=o-— (F.55) 
X (□'^F)|„,+(i„„), eif) H^'+'+^+'\x,y) . 

Since unbounded line integrals are involved, this expression is clearly non-causal. We 
shall now prove that (F.55) is smooth for x / y. Notice that the line integrals in (F.55) 
are supported on the hyperplane {q\q^ = 0}, e.g. 

/OO I- li 

V{ay + (1 - a)x) da = / V{q) 2TrS{q^) e~' 2 i^+v) . (F.56) 

-00 J (27r)'' 

For time-like ^, this hyperplane does not intersect the support of V, and thus (F.55) 
vanishes identically inside the light cone. Furthermore, (F.55) is clearly smooth in the 
region {^^ < 0} away from the light cone. Thus it remains to show that all the partial 
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derivatives of (F.55) vanish on the hght cone. The boundary values of these partial 
derivatives on the light cone {^^ = 0} involve integrals of the form 

/oo 
da a'+^^ (1 - af^^^ (a - a^)" (a^n"F)|„,+(i_„), ^ 
-oo 

with parameters ki,k2 > and multi-indices K, L. Similar to (F.56), these inte- 
grals are supported on the hypcrsurface {q \ q(^} = 0, and for ^ on the light cone 
(and C 7^ 0)) this hypersurfacc docs not intersect the support of V. We conclude 
that (F.54) coincides with (F.40), both in the case when supp V C {q^ < 0} and when 
supp V C {q^ > 0}. Linearity and an approximation argument near the light cone 
yield that (F.54) coincides with (F.40) for general V E S. 

It remains to consider the contributions when some of the derivatives in (F.53) act 
on the factor Q(z) in (F.52). The resulting expressions are of the form 

/yj4l POO 
t>(9) e-^ t M / dTe(g2)<5(")(z)ijf)(x,2/)P(r) e-^«^(^-^) (F.57) 
(2vr)* y_oo 

with integers a,/3 > and a polynomial P(r). Using the formula 

6{z) = - — : lim 



27ri e\o\z — ie z + ie ^ 
we can write the r-integral in terms of the complex integrals 

/oo -1 

Depending on the sign of q{y — x), the integration contour can be closed either in the 
upper or in the lower half plane, and the residue theorem yields expressions of the 
form 

lim . „ ,\ . i?W(a; - y) V{t) e-'^^^^'"^) with r = ±1 (F.58) 

e\0 [q lb l£y 2 

and j<2a + l, P<K<P + a. These expressions are well-defined distributions, 
and thus the g-integral in (F.57) is finite. Due to the powers of in (F.58), the 
resulting contributions to the operator product K^'-^VK^^'^ are non-causal. Since the 
factor H^''^ in (F.58) is a polynomial in ^ and V{q) in (F.57) has rapid decay, these 
contributions are also smooth. ■ 



The above lemma can be used to derive the light-cone expansion for the operator 
product ii:(OFr('"). 

Lemma F.4. For l,r >0 and V eS, 

(^(0 yTW)(x,y) 

X e(y° - x^) r('+''+'*+^)(x, y) + {contributions smooth for x / ?/) . (F.59) 
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Proof. Using (F.39) and the fact that Ha{x,y) is smooth in x and y according 
to (F.38), the hght-cone expansion (F.40) yields that for p,q>0, 

+ (contributions smooth for x ^ y) , (F.60) 

where S^""^ = hm^^O'^i"^- The main difference between (F.59) and (F.60) is that the 
factors K^'^^ and are replaced by corresponding factors The method 

of the proof is to realize these replacements by multiplying (F.60) with a suitable 
operator from the right. 

In preparation, we rewrite the operators S^'^ in (F.60) in terms of K^-^ as follows. 
Using that multiplication in position space corresponds to convolution in momentum 
space, we have for a > 0, 



/ 



K„(x,y)e(y°-x°)e-*^«d4C 
J- 



°° rf^ 2 ,r,2 . . . . PP 



-oo 



27r ^ ' ' ^ ^ ^ ^ ^ fcO _ 



PP 



ITT 2 |a;| k'^ — u 

and thus 



u>=-\/\k^+a ^ PP -j^ 

u,=-^/W+^ ZTT fc^ - a ITT 



Sa{x,y) = iTT Ka{x,y) e{y° - x°) . 
We differentiate with respect to a and let a \ to obtain 

5W(a;,y) = i7rir(")(a;,y) e(y°-a;°). 
Substituting into (F.60) gives 

(i(r(0 Y K^'^^){x, y) = (contributions smooth for x ^ y) 



1 oo ^ „Qo 

+ Tr- / ^^«e(y°-x°)a'(l-a)'-(a-a2)^ 

n=0 J-oo 



-oo 

X (□"y)|„,+(i_«), K('+'-+"+i)(x,y) . (F.61) 

The operator Ta, a > 0, is obtained from Ka by projecting on the negative-energy 
states, more precisely 

Ta = KaX, (F.62) 
where % is the multiplication operator in momentum space 

X{k) = -e(-A:°) . 

In position space, x has the kernel 

X{k) e-''^^-y^ = -— lim - 5Hy-x). (F.63) 

(ZTTJ^ zm e\o y" — x'^ — le 

If a is positive, the mass shell {/c^ = a} does not intersect the hyperplane {k^ = 0} 
where x is iiot smooth, and thus we may differentiate (F.62) with respect to a to obtain 

ri") = i^i") X , a > 0. (F.64) 
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However, difficulties arise in (F.64) in the limit o \ 0. Namely, the limit of the left side 
exist only after "regularizing" T^"^ by subtracting a polynomial in (y — x) (see (2.5.42, 
2.5.43). On the right side, the problem is that K^^\x,y) behaves polynomially at 
timelike infinity, whereas x(a;,y) decays for large (y" — x^) only as (y" — x^)~^, and 
so the product K^^'^x does not exist. To cure this problem, we insert into (F.63) an 
exponentially decaying factor by introducing for k > the kernel 

I g-K|j/0-a;0| 

XK{x,y) = -— lini-p -d (y-x). (F.65) 

2m £\o y" — x'^ — le 

The exponential factor changes the product with only by a contribution smooth 
in y — X, and thus 

T^\x,y) + (smooth contributions) 

= rH(x,y) = {k'<^^ x){x,y) = (K^'^) XK)(a^, 2/) + (smooth contributions) . 

The very left and right of this equation converge for a \ 0, and we conclude that 

(i^(") X«)(x,y) = T^''\x,y) + (smooth contributions). (F.66) 

We multiply the operator product on the left of (F.61) by the operator Xk- Apply- 
ing (F.66) and using that multiplying by a smooth operator gives something smooth, 

we get 

(if(') y i^W X«)(2^,y) = (K^'^ yTW)(x,2/) + (smooth contributions). (F.67) 

It remains to show that multiplying the right side of (F.61) by the operator Xk gives 
the right side of (F.59). When we multiply the summands on the right side of (F.61) 
by Xk, we get according to (F.65) a convolution in the time coordinate of the form 

/(^^^ . ^ I I 
dT F(x;y^-T, y) ^('+'■+"+1) (x; j/° - r, y) . (F.68) 
.oo T-l£ 

Here the function F stands for the line integral in (F.61); it is smooth unless y^—T = x^. 
We write the convolution integral (F.67) as 



^.^ r~ F{x-y^-r,y)-F{x,y) o _ ^ ^ 

£\0 7-oo T-l£ 



(F.69) 



/oo —K \t\ 
dT U,.^ yO _r,y). (F.70) 
-oo T -le 

The term (F.70) can be written as F{x,y) (_fc(^+''+"+i) XK)ix,y), and applying (F.66) 
gives precisely the summands on the right side of (F.59). Let us prove that (F.69) 
is smooth for x ^ y. Thus assume that x ^ y. Our above construction is Lorentz 
invariant in the sense that we may introduce the operator Xk in any reference frame, 
without influence on our operator products. Thus we can choose the reference frame 
such that X / y . The r-integral in (F.69) can be regarded as an integral along the 
straight line {{y^ — T^y)}- Since x 7^ y, this straight line does not intersect the point 
X. Due to causality of the integrand vanishes unless \{y^ — t) — x'^l > \y — x\. 
Thus if ^ is spacelike, i.e. ly'^ — < \y — x\, the integrand vanishes in a neighborhood 
of r = 0, and as a consequence the integral (F.69) is smooth in x and y. On the other 
hand if lies inside the cone \y^ — x^\ > ^ \y — x\, the straight line {(y" — r, y)} does 
for small enough r not intersect the hyperplane {z | = x^} where F is not smooth. 
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Thus the function F{x; y° — r, y) is smooth in a neighborhood of r = 0, and the mean 
value theorem yields that the bracket in (F.69) is smooth. This implies that (F.69) is 
again smooth in x and y. 

It remains to show that multiplying the contributions smooth for x ^ y, which 
are not specified in (F.61), by the operator Xki gives terms which are again smooth 
for X ^ y. Again using Lorentz invariance, we can choose a reference frame such that 
X ^ y. Then multiplying by Xk yields a convolution along the straight line {(y'^ — r, y)}, 
and this line integral does not intersect the point x. In this way, we can avoid inte- 
grating across the origin where the contributions in (F.61) may be singular. We get 
a convolution of Xk with a smooth function, and this is clearly finite and depends 
smoothly on x and y. ■ 



After these preliminaries, we are ready to compute the operator product QR in 
an expansion around the light cone. For the statement of the result, we need to 
analytically extend hj{^V} in (F.4) from a function on S'^ to a function on Minkowski 
space and also regularize it at the origin: As a smooth function on S"^, hj can be 
expanded in spherical harmonics. Since the spherical harmonics are the boundary 
values on S"^ of the harmonic polynomials on M^, we have the unique expansion 

oo 

hj{x) = Y.'^nms^s^ , (F.71) 
n=0 

where Vnix) are suitable harmonic polynomials of degree n. The smoothness of hj 
implies that the summands in (F.71) decay in n faster than any polynomial. As a 
consequence, the series in (F.71) converges absolutely for any x G M^, and we can even 
extend hj{x) to a unique function on C^. For £ > 0, we define the regularization /ij 
of hj by 

X 



m) = hj ( 1 . (F.72) 



Theorem F.5. Assume that the operator Q satisfies the weak evaluation for- 
mula (F.4)- Then the operator product QR with R according to (F.31) with p > —I 
and g G 5(M^ x M^) has the expansion 

Lmax 9max 

{QR)ix,y) = Ylog<^{Ep)E^-^Yl (^-73) 

L=2 g=0 J 

X {-2n) f; hm Hda (1 - a)^ {a - a'r ^1 {jW^^ 9{z. y)) 

n=0'^^ ~ J \z=ay+{l-a)x 

X e(y° - x^) r(P+"+i)(x,y) mod V^'^'^^'^ [d^) R{x,y) + o{Ep) (F.74) 
+ (contributions smooth for x ^ y) + (higher orders in {lmacroFp)~^ ) (F.75) 
with C = z — X. 

Proof. We first consider the case p > 0- We fix two space-time points xq, yo and 

set 

QiO ■■= Q{xo,xo + 0, RiO ■■= Riyo-Lyo)- 

We regard Q(^) and -R(^) as the integral kernels of corresponding homogeneous op- 
erators, which with slight abuse of notation we denote again by Q and R. Then 
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in momentum space (F.ll), the operator Q has the expansion of Lemma F.l and 
Lemma F.2. Let us compute the product Qj^R. According to its construction in 
Lemma F.l, Qj^{$,) behaves for ^ 7^ hke 

and furthermore its Fourier transform Qj^{k) vanishes inside the lower mass cone (cf. 
Lemma F. 2). The distribution 

lim K^^^ (x, y) [-Att'^i (y° - x° - ie)-^+^ hj{0 ^"^j (F.76) 

has these two properties, as one sees immediately from (F.34) and when computing the 
Fourier transform of (F.76) with contour integrals. Indeed, a short calculation shows 
that (F.76) even coincides with Q^/{£,). 
We introduce for e > the potential 

V,{z) = -47rH{z'>-x''o-ie)-'^+'h^j{z-xo){z-xoy g{z,yo), 

where hj is the rcgularization of hj (F.72). This potential is a Schwartz function, and 
thus Lemma F.4 yields that 

(i^(0) ye T^P)){^X,y) 

= da {1- aria- a^T (□"^.)|a.+(i-a). 

n=0 ■^-^ 

X e(y° - x°) r(P+"+^)(a;, y) + (contributions smooth for x ^ y) . (F.77) 

We now set x = xq, y = yo and take the limit e \ 0. On the left side of (F.77), we can 
use that (F.76) = Qj^iO to obtain the operator product {Q^^ R){xq, yo). Furthermore, 
due to our rcgularization of /ij, the factors {W^Ve){z) on the right side of (F.77) are 
of the form smooth function times [z^ — x^ — ie)^\ and thus the limit e \ exists in 
each line integral in (F.77). We conclude that {Qj^R){x,y) coincides precisely with 
the series (F.74). Since in (F.77) we integrate across the origin, the higher orders 
in ((r + \t\)Ep)~^ in (F.4) (or equivalently the higher orders in k/Ep in (F.16)) yield 
contributions of higher order in (/macroFp)"^. Finally, calculating modulo polynomials 
in (F.17) means in position space that Q{x,y) is determined only modulo partial 
derivatives of 6^{x — y), and this gives rise to the term mod 'P^''\'^^~^ (dx) R{x,y) 
in (F.74). This concludes the derivation of (F.73-F.75) in the case p> 0. 

In the above derivation we neglected the contributions smooth for x ^ y, which 
are not specified in (F.77), implicitly assuming that they remain smooth in the limit 
£ \ 0. This is justified as follows. As £ \ 0, the potential T4(z) becomes singular 
only at z = xq. Thus if the contributions smooth for x ^ y in (F.77) had a non- 
smooth limit, the resulting non-smooth contributions to (Q^/-^ R)(xo, yo) would depend 
on g(z,yo) and its partial derivatives only at z = xq; i.e. they would be of the form 

dig{xo,yo)Wj{xo,yo), (F.78) 

where Wj are distributions independent of g. Suppose that the Fourier transform g of 
g{xo + ., yo) is supported inside the upper mass shell {k\k'^ > 0, k^ > 0}. Then due to 

our (— i£)-regularization, also the support of is inside the upper mass shell. As a con- 
sequence, the Fourier transforms of the distributions K^'^^xo, ■) and (V^T^'^s (P))(., t/o) 
have disjoint supports, and the left side of (F.77) is zero. Furthermore for ^ close to 
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the hght cone, the unbounded hne integrals in (F.77) vanish (notice that they are sup- 
ported on the hypersurface {A; | A;^ = 0}). We conclude that if g is supported inside the 
upper mass shell, then the contributions smooth for x ^ y va (F.77) are zero. Taking 
the limit e \ yields that the contributions (F.78) vanish if g is supported inside the 
upper mass shell, i.e. 

'' d'^kk^ g{k)Wj{xo,yo) = (F.79) 



for all g with supp g d {k\k^ = {),k^ > 0}. Since the polynomials restricted to the 
upper mass cone are linearly independent, (F.79) implies that the contributions (F.78) 
are all identically zero. 

The extension to the case p = —\ follows exactly as in [F6, Lemma 2.2]: We pull 
one derivative out of the operator product, 

j d^zQ{xQ,z)g{z,yo)T^-^\z,y) 

^'^=^ -2 A ld^zQ{xo,z)fiz,yo)T(^\z,y)^y=y^, 

substitute for the integral the expansion (F.73-F.75) for p = 0, and differentiate 
through. ■ 



The above theorem gives the clue for understanding the operator product QP as 
well as the commutator [P, Q] , as we shall now explain. Using that the product of Q 
with a smooth operator is smooth, we can write the operator product QP according 
to (F.30) in the form 

oo 

{QP){x,y) = ^{Q PP){x,y) + (smooth contributions) (F.80) 

p=-i 

with 

P''{x,y) = gp{x,y)T^\x,y). (F.81) 
The summands in (F.80) are precisely of the form considered in Theorem F.5, with 
the only exception that the functions g-p in (F.81) in general have no rapid decay at 
infinity. Fortunately, the behavior of the functions gp at infinity is of no relevance to 
us, and we can apply Theorem F.5 to each summand in (F.80) using the following 
approximation argument. For fixed x and y, we choose a Schwartz function -q which 
is identically equal to one on a compact set K 3 x,y. Then the function (gpf]) has 
rapid decay at infinity, and Theorem F.5 applies to the operator product (F.31) with 
g := gpT]. In the discussion below leading to Corollary F.6, the function g{z) enters 
only for z in a neighborhood of x and y. In this neighborhood, g and gp coincide, and 
thus the dependence on r] drops out. This shows that the behavior of the function gp 
at infinity is indeed of no importance for Corollary F.6 below. 

Let us briefly discuss the expansion (F.73-F.75). First of all, we point out that 
we calculate modulo terms of the form d^R{x,y) with \K\ < L — |J| — 3. This 
corresponds to the fact that we have no information on the behavior of Q near the 
origin. For generic regularizations or simple regularizations like a cutoff in momentum 
space, the terms R{x, y) will not be zero. Thus in this case, the operator product 
QR does not vanish, even when we take for R a Dirac sea of the vacuum; this is 
in agreement with our consideration in momentum space after (F.29). However, the 
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situation is much different if we assume that we have a regularization where the terms 
d^R{x, y) all vanish. Namely, if we then take for R a Dirac sea of the vacuum, for 
example 

R{x,y) = '-^T(-'\x,y), 

the integrand of the line integral in (F.74) is a rational function with poles in the up- 
per half plane, and QR is zero (up to the contributions not specified in Theorem F.5). 
This corresponds to our observation made after (F.22) that the poles of Q cannot be 
detected when testing with solutions of the free Dirac equation. The regularizations 
for which the terms R{x,y) vanish are just the optimal regularizations introduced 
after (F.29). The main advantage compared to our earlier consideration in momentum 
space after (F.29) is that the expansion (F.73-F.75) tells how the macroscopic per- 
turbations of P and Q effect the operator product. In particular one sees that, when 
taking into account the macroscopic perturbations, the operator product QP does in 
general not vanish (even for optimal regularizations), and thus the problem of disjoint 
supports disappears. 

In an interacting system, the factor g{z, y) in (F.74) is composed of the bosonic po- 
tentials and fields. Thus in the generic situation, the line integrals in (F.74) will vanish 
only if the operator Q is identically zero. In order to make this argument clearer, it is 
useful not to think of P as a fixed object, but to consider small dynamic perturbations 
of P. More precisely, we consider perturbations of P induced by perturbations of the 
bosonic potentials of our physical system. In order not to disturb the Eulcr-Lagrange 
equations, these perturbations must not be arbitrary, but should satisfy the physical 
equations; a typical example are perturbations by an electromagnetic wave. Thus we 
consider variations of our system by small, physically admissible perturbations of P 
and study the effect on the operator product QP. We refer to this procedure for an- 
alyzing the operator product that we test with physical perturbations of P. Clearly, 
the requirement that the perturbation should satisfy the physical equations is a strong 
restriction (in particular, such perturbations are not dense in the Lp' topology). The 
reason why it is nevertheless a reasonable concept to use physical perturbations for 
testing is that these perturbations enter into (F.74) only along the one-dimensional 
line {\y + (1 — A)a;}. In the example of the perturbation by an electromagnetic wave, 
the electromagnetic field appears in the function g in (F.74), and by changing the 
location and amplitude of the wave, we can completely determine the function hj{p) 
as well as the order L — | J| — 1 of the pole of the integrand at the origin. Notice 
furthermore that the summands in (F.73) scale in the Planck energy exactly as the 
corresponding summands in the expansion (F.4). We conclude that by testing the 
operator product QP with physical perturbations of P, we can reconstruct the weak 
evaluation formula (F.4) completely. 

Next, it is instructive to consider the behavior of the summands in (F.74) as y — x 
gets small. If y — x is scaled like (y — x) — > A (y — x) with A > 0, the variable trans- 
formation a — \~^a shows that as A \ 0, the line integral blows up like X~i'^'^+P+'^) . 
On the other hand, the factor 7"(p+"+i) goes to zero like X^^+'^P in this limit. Thus 
each summand in (F.74) scales like A~^^^. It is remarkable that, no matter how large 
the order of the pole of Q at the origin is, the operator product QP has at the origin 
a pole of at most the order one. The reason is that in (F.74), we integrate over the 
pole of Q, and this regularizes the singularity at the origin. Since a pole of order one 
is integrable in three space dimensions, we do not need to study the operator product 
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{QP){x, y) at the origin x = y. According to the hght-cone expansion of Theorem F.5, 
the information contained in the weak evaluation formula (F.4) is retrieved in the op- 
erator product QP by considering the singularities on the light cone away from the 
origin. 

The expansion of Theorem F.5 immediately allows us to study also the commutator 

[P,Q]. Namely, by taking the adjoint, {QR){x,y)* = {RQ){y,x), the formula (F.73- 
F.75) applies to the operator product PQ as well, and we can take the difference 
[P,Q] = PQ — QP. The key observation is that in the product {QR){x,y), the pole 
of Q at the origin appears in (F.74) together with the factor g{z,y) and 2; x, 
whereas in the product {RQ){x,y), this pole is multiplied by g{x,z) and z y. Thus 
when testing [P,Q\{x,y) with perturbations of P, one can distinguish between the 
contributions from PQ and QP by considering perturbations which are localized near 
y and x, respectively. 

These results are summarized as follows^. 

Corollary F.6. To every order in Ep, the poles ofQ{x,y) at the origin x = y 
can be detected in the commutator [P,Q] by testing with physical perturbations of P. 

We close with a general comment on the significance of the constructions in this 
appendix. Due to the problem of disjoint supports, we could make sense out of the 
commutator [P, Q] only after taking into account the macroscopic perturbations of 
P and Q. As a consequence, the relevant contributions to the commutator [P, Q] 
are by several orders of {lmacroEp)~^ smaller than expected from a simple scaling 
argument. This can be interpreted that the causal structure of Minkowski space and 
the structure of the Dirac seas, which are the underlying reason for the problem of 
disjoint supports, have a tendency to making the commutator [P, Q] small. In this 
way, the causal structure and the structure of the Dirac seas seem to correspond nicely 
to Euler-Lagrange equations of the form [P, Q] = 0. 



Online version: A simpler and cleaner method to obtain this result is to use the so-called method 
of testing on null lines as worked out in §5.2 in the paper arXiv:0908.1542 [math-ph]. 
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Perturbation Calculation for the Spectral Decomposition 

of P{x,y) P{y,x) 

In this appendix we shall develop a convenient method for analyzing the eigen- 
values and spectral projectors of the matrix A^y = P{x,y) P{y,x) and compute all 
contributions to the eigenvalues needed for the derivation of the effective gauge group 
in Chapter 7. Our strategy is as follows. We decompose the fermionic projector as 

P = Po + AP 

with Pq according to (6.1.1). This gives rise to the decomposition of A 

A = Ao + AA (G.l) 

with 

^0 = Poix,y)Po{y,x) (G.2) 
= APix,y)Po{y,x) + Po{x,y)AP{y,x) + AP{x,y)AP{y,x). (G.3) 

The eigenvalues and spectral projectors of were computed explicitly in Chapter 6, 
see (6.1.9, 6.1.10). On the light cone, Po{x,y) has singularities of order 0{{y — x)^^), 
whereas AP{x,y) = 0{{y — x)~^). Likewise, AA is compared to Aq of lower degree 
on the light cone. For this reason, can be treated perturbatively in the sense 
that the eigenvalues and spectral projectors of A can be expressed to any given degree 
on the light cone by a finite order perturbation calculation. Apart from the purely 
computational aspects, the main difficulty is that Aq may have degenerate eigenvalues, 
and in this case we need to carefully analyze whether the degeneracy is removed by 
the perturbation. Our method is to first compute projectors on invariant subspaces of 
A (§G.l). Analyzing the perturbation on these invariant subspaces will then give the 
spectral decomposition of A (§G.4). 

G.l. Perturbation of Invariant Subspaces 

We write the spectral decomposition of Aq as 

K 

^0 = ^ Fk 

k=l 

with distinct eigenvalues Afe and corresponding spectral projectors F^. As in §5.3 we 
use the convention Ai = 0. Clearly, the F^. are the sum of the spectral projectors 
counting multiplicities, 

n,c,s with Anc3=Afc 

with Xjics and Fncs according to (6.1.10). Since the perturbation AA will in general 
split up the degenerate eigenvalues, we cannot expect that by perturbing Fj. we obtain 
spectral projectors of the matrix A. But we can form projectors on the space 
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spanned by all eigenvectors of A whose eigenvalues are sufHciently close to Aj^. The Gk 
are most conveniently introduced using contour integrals. We choose £ > such that 

|Ai — Ajl < 2e for alli,j = l,...,K and i ^ j. 

Then we set 

Gk = ^ (f {z- A)-^ dz , (G.2) 



27ri J\z-\k\=e 

The Cauchy integral formula shows that Gk is indeed a projector on the desired sub- 
space. 

The integral formula (G.2) is very useful for a perturbation expansion. To this 
end, we substitute (G.l) into (G.2) and compute the inverse with the Neumann series, 

Gk = {z-Aq- AA)-i dz 



— i (1 - (z - Aq)-^ AA) ^ {z - Aq)-^ dz 

J\z-\^\=e 

oo 

-i J^((z-Ao)-'a^)" (^-^o)-'c^^• 
'''"^ J\z-\k\=£ „-n 



Interchanging the integral with the infinite sum gives the perturbation expansion, 

oo » 



— 27ri J\z^x,\=e 

where n is the order in perturbation theory. After substituting in the spectral repre- 
sentation for {z — Aq)~^ , 

{z-Ao)-' = E-^, (G.4) 

Z — A/ 

1=1 ' 

the contour integral in (G.3) can be carried out with residues. For example, we obtain 
to second order, 

Gk = Fk + i;i-^(^feA^^/ + Fi^AFk) + 0((Ayl)3) 



X {Fk AA Fi AA Fm + Fi AA Fk AA F^ + A^ F^ AA Fk) 

X {Fk AA Fk AA Fi + Fk AA Fi AA Fk + Fi AA Fk AA Fk) . (G.5) 

To order n > 2, the corresponding formulas are clearly more complicated, but even 
then they involve matrix products which are all of the form 

Fk, AA Fk, AA - - Ffc„ A^ Fk„^, . (G.6) 

Substituting in (G.l) and expanding, we can just as well consider matrix products of 
the form (G.6) with the factors Fk replaced by Fncs- Furthermore, for the computation 
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of the eigenvalues we need to take the expectation values of Gk with certain matrix 
elements of AA. This leads us to traces of matrix products of the form 

Tr (F^icisi AAi Fn^c2S2 • • • Fnidsi ^Ai) (G.7) 

with I = n+1. We refer to a trace of the form (G.7) as a matrix trace. Our first task 
is to develop an efficient method for computing matrix traces (§G.2 and §G.3); after 
that we will proceed with the calculation of the eigenvalues of A (§G.4). 



G.2. Factorization of Matrix Traces 

If one attempts to calculate a matrix trace (G.7) directly by substituting in the 
formulas of the light-cone expansion, the resulting expressions become so complicated 
and involve so many Dirac matrices that they arc almost impossible to handle. We 
shall now simplify the situation by giving a procedure which allows us to factor matrix 
traces into a product of so-called elementary traces, which are much easier to compute. 
According to (G.3), we can assume that each factor AAj in (G.7) is the product of a 
contribution to P{x,y) with a contribution to P{y,x). Denoting the contributions to 
P{x, y) by Bj and using that the corresponding contributions to P(y, x) are obtained 
by taking the adjoint with respect to the spin scalar product, we can write each /^Aj 
in the form 



Inserting the completeness relation 



and expanding gives for (G.7) a sum of terms of the form 

Tr (-Fnicisi Bi Fn2C2S2 B2 ■ ■ ■ Fuk-iCk-iSk-i Bk-l PukCkSk B)^) (G.l) 

with indices (rij, cj) (which arc in general different from those in (G.7)) and k = 21. 

In order to handle the sector indices in (G.l), we introduce operators Kn^^n2 which 
act on the sector index and map sector 77,2 to sector ni, i.e. in components 

{Knin2)n' = ^ni ^n'n2 ■ (G.2) 

Then 

Fncs = Knl Fics Kin ■ (G.3) 

If we substitute this relation into (G.l) and combine the operators K. and Bj to "new" 
operators Bj, we obtain a matrix trace again of the form (G.l), but with all indices 
rij equal to one. Therefore in what follows we can restrict attention to the case of 
one sector and omit the sector indices. The generalization to several sectors will be 
straightforward by inserting operators K. into the end formulas. 

We choose a space-like unit vector u which is orthogonal to ^ and ^. Then the 
imaginary vector v = iu satisfies the relations 

vj = = Vj^ , v'^ = 1, V = -V . (G.4) 

An explicit calculation using (6.1.10) yields that 

Fr+ = iFL+i, Fl- = -JiFL+ii, Fr^ = -JFl+^. (G.5) 
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Substituting these formulas into (G.l), we obtain an expression involving only the 
spectral projector -Fl+, namely 

(G.l) = TV {Fl+ Ci Fl+ C2 • • • Fl+ Ck) (G.6) 

with suitable matrices Cj. Since the Fl+ are projectors on one-dimensional subspaces, 

Fl+CFl+ = Tr{FL+C)FL+. 

By iteratively applying this relation in (G.6), we get the product of traces 

TV {Fl+ Ci) TV {Fl+ C2) • • • TV {Fl+ Ck) • 

If we express the matrices Cj explicitly in terms of Bj and Bj, we obtain the following 
factorization formTila, 



Tr (Fcisi Bi Fc2S2 P2 
= F^,ltl{B,)F^l^l{B*^] 



Bl) 



(Bk-i) F^^^liBl) , 



(G.7) 



where Fg-sj 



are the so-called elementary traces defined by 



F^^{B) = Tr{F+XLB) 

Ff^{B) = Tr{^F+ixLB) 

F^liB) = ^TriF+i^XLB) 

F^^iB) = ^Tr{^F+^XLB) 



F^^iB) = TriF+ixLB) 
F^?{B) = Tri^F+XLB) 

FL^{B) = -^Th:{F+ixLB) 

F^^{B) = - Ti{iF+iixLB) 
z 



(G.8) 



These formulas are also valid for the opposite chirality after the replacements L 
The elementary traces of B* are obtained by taking the complex conjugate, 

F£i(B*) 



R. 



F^^iB*) = F^^iB) 



F^^{B*) 



fWb) 



F^^{B*) = F^^{B) 
FLt{B*) = Ff^{B) 



F^^{B*) = F^fiB) 



(G.9) 



pLR(^B*) = F^^{B) . 

The relations (G.7-G.9) are verified by a straightforward calculation using (6.1.10, 
6.1.8, G.4). 

To summarize, the above procedure reduces the calculation of the matrix trace (G.7) 
to the computation of the elementary traces (G.8) for the contributions B to the light- 
cone expansion of P{x,y). Taking the complex conjugate (G.9), one obtains the ele- 
mentary traces of the corresponding contributions to P{y,x). By applying (G.7) and, 
in the case of several sectors, by suitably inserting the operators K., every matrix trace 
can be written as a linear combination of products of elementary traces. 

G.3. Calculation of the Matrix Traces 

We decompose AP{x,y) into its odd and even parts, denoted by Bo and Be, 

AP{x,y) = Bo{x,y) + Be{x,y) . 

Explicit formulas for the fermionic projector in the presence of chiral and scalar poten- 
tials are listed in Appendix B. For the purpose of this paper, only the contributions 
involving the mass matrices l^/zj and their derivatives are of importance. But for 
completeness and for later use, we will also compute the contributions which contain 
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the chiral field strength and the chiral currents. However, we will omit all contribu- 
tions quadratic in the field strength. Namely, these contributions are related to the 

energy-momentum tensor of the chiral fields, and it is therefore reasonable to post- 
pone their analysis until gravity is also taken into consideration. Thus the phase-free 
contributions relevant here are 

XlB, = ^XLmT^'\x,y)^ J' dz (DjYl) 

+XLmr(0)(x,y)yL(x) + ©(loglC^U^) 
ry 

XL Bo = |xLm2rW(a;,y)^ j cIzYlYr 

+iXL T(i) {x, y) r dz [0, 1 | 0] Yl 7^' {D^Yr) 

J X 

+iXL rr? t(i)(x, y) dz [0, 1 | 0] i\D^Yl) Yr 

J X 

-iXL rr? r(i)(x, y) Yl{x) dz i^D^Yr) 

J X 

+XLTW(x,y)r I' dz [0,l|0]Vi^,f 

J X 

-\xLT^'\x,y)iC [d^ [0,0|l]jf 

-XLr«(x,2/)r r dz [o,i|i] (a^) 

Jx 

-XLT^'Hx,y) r dz [0,2 lOjjtl" 

Jx 

+^o(r') + i'F^^eo{c') + o{Fi) -h o{\og\e\e)- 



A straightforward calculation yields for the elementary traces 

F^^Po) = (deg<l) = ix^(zr(-^)) (G.l) 

F^^iPo) = (deg <2) = '-Xl t{-'^ (G.2) 

F^i^iBe) = (deg < 1) = YLix) T[(°) + (deg < 1) (G.3) 

F^^{B,) = (deg<0) (G.4) 

F^^{Be) = (deg<l) (G.5) 

Fff:(Se) = (deg < 1) = YL{y) T^^^ + (deg < 1) (G.6) 

F^^iBo) = (deg<l) (G.7) 

= r dz [0, 1 I 0] F/, r|;) + (deg < 1) (G.8) 

J X 

+^ eijkl e f r dz [0, 1 I 0] Ft iU T^S^) (G.9) 
z z J J. 
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e'^^' i^j + Uf - Vk r dz CF^i T^J) (G.IO) 

F^^{Bo) = (deg<0) 

= J'dz Yl Yr {{z r(°)) + 4 T^])) + (deg < 0) (G.ll) 

-2i Yl{x) YR{y) T^l^ (G.12) 

[ [0' I 1] A {{z rg)) + 8 tI^) (G.13) 



F^^iBo) = (deg<l) 



2 

Ff^(5o) = (deg<l 



cLzYl Yr t[^^ + (deg < 1) (G.15) 
£dz[0,0|l]iiTg) (G.16) 

^^^■^'7^ /J^f (e'Tg)) (G.17) 



= I'dz [1,0 \ 0] F/, + (deg < 1) (G.18) 

J X 

+l'''''^^^^iykiTl^ (G.19) 
eijki e ? [0, 1 I 0] Fi- (e^ rg)) (G.20) 

J X 



Z — Z 



^ ^ e'^'' + - ^i^f ) £ dz CFni . (G.21) 



Here the totally anti-symmetric symbol Cijki appears because we applied the identity 
T^iXL/R i^H) = 2 {{ab){cd) + {da){bc) - {ac){bd)) T 2z eijM aVd'd^ . 

Therefore, the corresponding formulas for the opposite chirality are now obtained by 
the replacements 

L < — > R , eijki — > - eijki ■ (G.22) 

The elementary traces of the adjoints are computed via (G.9). All other elementary 
traces vanish. 

Applying (G.7) and the degree estimates for the elementary traces and omitting all 
terms quadratic in the field strength, we can factor and estimate the following matrix 
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traces, 



Tr{FL-AA) = Ff^(Po)FffTO + Ff:«(i?o)F_«f(Po*) 



+ F^^{Be) F^^{B:) + (deg<2) 



+ F^!l{Be) Ft_{Bl) + (deg < 2) 




(G.23) 



(G.24) 
(G.25) 
(G.26) 



T,{FlsAAFlsAA) = Tr{FLsAA)Tr{FLsAA) = (deg < 5) 
Tr{FLsAAFRsAA) = (deg < 5) 
Tr{FL+ AA Fr^ AA) = (F^^ (Se) F^^{P^) + F^^{Po) F^^iB^)) 




x{F^^{Be) F^^{P;) + F^t(Po) F^i^{B*)) + (deg < 5) 
TV(Fl_ AA Fr+ AA) = {F^^iBe) F^^iP^) + F^^{Po) F^^{Bl)) 



(G.27) 



x(Ff«(Be) F^t{P5) + Fff(Po) Ftt{Bl)) + (deg < 5) 
Tr(FL+ AA Fl- AA) = = Tr(FL_ AA Fl+ AA) . 



(G.28) 
(G.29) 



If we consider more generally the matrix trace of order Z, factorization gives a linear 
combination of products of elementary traces as in (G.7) (with k = 21). Let us estimate 
the degree of each of these products. Clearly, the number of factors equals the 
number of factors -Fl^, we denote the number of such pairs by p. Furthermore, let q be 
the number of factors F'.{AP') (where AP' stands for cither AP or AP*). According 
to (G.3), each AA contains at least one factor AP', hence q > I. The number of factors 

P++ and F^ is 2{l—p), and we saw above that each of these factors must involve AP', 

thus q > 2{l—p). Adding our two upper bounds for q gives the inequality q+p > 3//2. 

To estimate the degrees we first note that the degree of the pair P.^ (Pq) P1_|_(Po) is 

three, and is decreased at least by one each time a Pq is replaced by AP'. The total 
number of factors P_^_(AP') and Pl^(AP') is q — 2{l — p). On the other hand, the 
degree of each factor F'j^_^ and Pl_ is at most one. Hence the degree of the matrix is 
bounded from below by 2>p — {q — 2{l — p)) + 2{l — p) = Al — [q + p). Substituting in 
our above lower bound for g + p gives the degree estimate 



The above formulas are valid in the case = 1 of one sector. The generalization 
to several sectors is done by inserting suitable operators K. into the traces. This has 
no effect on the degree on the light cone, and thus the estimates of the matrix traces 
in (G.23-G.30) hold in the general case as well. We substitute the above results for the 
elementary traces (G.1-G.21) into (G.23-G.28) and insert the operators K. to obtain 
the following explicit formulas: 




for I > 3. (G.30) 



Tr(P„i+ AA) 



(deg < 2) 




(G.31) 



(G.32) 



(G.33) 



(G.34) 



288 



G. PERTURBATION CALCULATION FOR THE SPECTRAL DECOMPOSITION 



+- r dz [0, I 1] TV5 [In jl Xr] {{z rg)) + 8 T«) r(y^) (G.35) 
r dz [0, I 1] Tr5 {In Xl fn] {z t(-^)) ijf (G.36) 

+ J 6,,-,z C/ TV. {/n i^f } {z T^-^)) ) (G.38) 
Tr(F„i_A^) = (dcg<2) 

+Trs {/„ yL(y) M^)} t[^^ (G.39) 

+^ j^^ dz TV5 {/n Yr Xr] rg^ Tjl^j'^) (G.40) 

+^ j^" TV5 [in Xl Yr Yl] T^-^^^ {{z t[^^) + 4 Tj^^]^) (G.41) 
-IVs {/„ yR(y) Yl{x)] T^-^^ t[^^ (G.42) 
-k. /"^ dz [0, I 1] TV5 {/„ fR) t[-'^ {iJWh + 8 (G.43) 

+^ £ [0, I 1] Tts [In jl Xr] r(°) (zT^-^)) (G.44) 

TV(F„L+ AA) = (deg<5) 

X /„, Xr - Xr n{y) t([') I]f ) } (G.47) 

TV(F„l_ AA A^) = (deg<5) 

-\ [in (yUv) Xl ^-Xl Uv) t{-;^ ^) 

X In' {yr{x) Xr T[i°) {zT{-^'^) - Xr Yl{x) {z T^"^)) T^^ | (G.48) 

G.4. Perturbation of the Non-Zero Eigenvalues 

In §7 we calculated the eigenvalues Xncs of A in the presence of chiral and scalar 
potentials to the leading degree 3, (6.1.10). Now we shall compute the contributions 
to the non-zero eigenvalues of degree two, denoted by AA„cs5 = 1) • • • ) 7 (the kernel 
of A will be considered in §G.5). To this end, we need to analyze the matrix A on the 
invariant subspaces ImG^. First, we choose for fixed A; > 1 a convenient basis of ImG^ 
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as follows. The degeneracy of the unperturbed eigenspace ImF^ can be described by 
the index set /, 

/ = {(ncs) with Xncs = Afe} . (G.l) 
Note that, according to (6.1.10), s is the same for all elements (ncs) G /, provided 
that the eigenvalue is non-zero. The index c, however, may take both values L and R, 
giving rise to the partition of / into II and Ir, 

Il/r = {{ncs) G / with c = L/R} . 
The set / can be used to index a basis of F^; namely we choose 

(0ncs)(ncs)e7 with 7^ 0„cs € Im F„cs • (G.2) 

It is convenient to assume that the basis vectors are related to each other by 

4'n'cs ~ ^n'n 4'ncs > 'Pn'cs ~ ^n'n ^ 4'ncs ! (G^-S) 

this can clearly arranged according to (G.3-G.5). Since Fk projects onto a null 
space, the inner product of any two basis vectors 0„cs vanishes. Thus in order to 
be able to evaluate vectors in ImF^ using the scalar product, we choose a "dual basis" 
('^'"'')(nc.)G7 of ImF* given by 

</.--eImF4,, r'"' = K^'n^"', r"'' = K^'ni . (G.4) 

The basis vectors and their duals are orthogonal in the sense that for (ncs) / (n'c's), 

<</'"^' I 0n'c's> = <Kcs^''^"\F'n'c's(l)n'c',> = <4>'"'' \ Fncs Fn' d s (l^n' c' s> = 0. 

We normalize the basis vectors such that 

I (l>n'c's> = C for all (ncs), (n'c's) G /. (G.5) 

Next we introduce a basis (V'ncs)(ncs)G/ of the invariant subspace ImGjfc by applying 
the projector Pj. to the ^ncs^ 

i^ncs = Gk (pncs ■ (G.6) 

Finally, we introduce a basis (V'"''^*)(rtcs)e/ which is dual to (ipncs)- We must be care- 
ful because projecting on Im(G/j) and Im(Gfc), respectively, does not preserve the 
orthonormality; more precisely, 

Sn^'s = <Glcf>^'^'\^r^'c's> = <Glcf>^'^'\G''cf>n>c's> 

in general 

= < </>"^^ I G^^ I </)„v.> 7^ S:i,S',,. (G.7) 
But S is a perturbation of the identity, and thus it can be inverted within the perturba- 
tion expansion by a Neumann series. This makes it possible to introduce ('0"''^*)(ncs)e-f 

by 

^ncs ^ ^ {S-'y-^,^Gir'''''- (G.8) 

A short calculation shows that this basis of ImGj^ is indeed dual to {ipncs) in the sense 
that 

<V''"'' I ^n'c's> = ^n' for all (ncs), (n'c's) G /. (G.9) 

Using the basis {ipncs) and its dual (■i/'"'^*), we can write down matrix elements of 

A, 

Al?l,^ = Kij'^"' I A I Vn'c's> for (ncs), (n'c's) G I. (G.IO) 
From the orthonormality (G.9) one sees that ^J^f^/^ is indeed a matrix representation 
for A in the basis (V'ncs)) and thus the eigenvalues of A on the invariant subspace 
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ImGjfc are obtained simply by diagonalizing this matrix. In the unperturbed case (i.e. 
if = 0), the matrix ^^f^'s simplifies to 

K'ds = <0"'^^ I ^0 I </>nVs> = <</>"'^' I ^0 -FnVs </>n'c'.> 
= Afe <0"^^ I <j)n'c's> = Afe 51, 51, , 

in agreement with the fact that ImF^ is an eigenspace of Aq corresponding to the 
eigenvalue A^. Thus we see that the matrix elements A^f^'s leading order on the 

light cone of degree 3. In the following theorem we compute the matrix elements up 
to contributions of degree < 2. 

Theorem G.4.1. We consider the fermionic projector in the presence of chiral and 

scalar potentials (6.0.5) and in composite expressions disregard all terms quadratic in 
the field strength. Then for all k = 2, . . . , K and (ncs), (n'c's) G I, 

K'!'s = Xk5:,6',, + 5',,Tr{F,,,...AAK,,r,^ 

J2 ^ (^ncs AA Fi AA + (deg < 2). (G.ll) 

Proof. We begin by computing the matrix S, (G.7), and its inverse. This calcu- 
lation will also illustrate how the relations (G.3) and (G.4) make it possible to rewrite 
expectation values as matrix traces and thus to apply the results of §G.2 and §G.3. In 
the case c = c', we obtain from (G.7) and (G.3), 

'S'n'cs = Kcf/^'^^ I Gk I ^n'cs> = <(p'^^^ \ Gk | K^'n 4>ncs> 

(G.2,G.4) 



(*) 



Kn'n I ^ ncs 'VncS' 
'^4' I ^ncs G}~ Kn'n Pjics \ 

Tr {Fncs Gk Kn'n) <(l>^'^ \ F^cs | 'i>ncs> 

Tr {Fncs Gk Kn'n) <4'^'^ \ (l>ncs> ^ = ^ {Pncs Gk Kn'n) , 



where in (*) we used that Fncs projects on a one-dimensional subspace. If we sub- 
stitute the perturbation expansion for Gk, (G.3), into the obtained matrix trace, the 
estimate (G.30) shows that the orders n > 2 yield contributions to S of degree < — 1. 
Thus it suffices to consider for Gk the second order expansion (G.5). This gives 

5n'c. = <^n' - E n i^ncs AA Kn'n) + (deg < -1) . (G.12) 

Note that of the matrix trace appearing here we need to take into account only the 
leading contributions of degree 5; these are easily obtained from (G.47) and (G.48). 
In the case c 7^ c', wc obtain similarly 

•^n'!. = <0"" I Gk Kn'n i Kcs> = Tr {Fncs Gk Kn'n i) • 

We again substitute in the expansion for Gk (G.5). As a consequence of the additional 
factor y, the contribution to zeroth order in A A now drops out. The first order 
contribution to Sn?^g is 

, \ {Fncs ^AFi Kn'n i) = T {Fncs AA Fn'cs Kn'n i) 

l^^l^ ~ M ^ncs — ^n'cs 

^ {Fncs AA Kn'n i) = (deg < -1) , 



Ancs Ajj'cs 
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because according to (G.8) and (G.7) the last matrix trace has degree < 1. Here we 
imphcitly assumed that A„cs 7^ -^n'csi because otherwise we clearly get zero. A straight- 
forward calculation using the factorization formula (G.7) as well as the estimates for 
the elementary traces following (G.l) shows that the second order contribution to Sl^f^^^ 
also is of degree < — 1. We conclude that 

= (deg < -1) . (G.13) 

Now we can take the inverse of the expansions (G.12) and (G.13). This gives 

{S-')!;:^^,, = 6'^,6:^,+S'^,Y,^—^Tv{Fncs^AFiAAKn'n)+{<i^g<-l). (G.14) 
We next compute the expectation values 

«f)'"'' \AGk\ (t>n'c's> 

up to contributions of degree < 2. The method is the same as for the above calculation 
of the matrix S. In the case c = c', we obtain the following matrix trace, 

<(/)"^" \AGk\ (t)n'cs> = I A Gk Kn'n \ (t>ncs> 

_ ^^ncs I p^^^^ Kn'n F-ncs \ <t}ncs> = Tr {Fncs A Gk Kn'n) ■ 

Substituting in (G.l) and (G.3), the estimate (G.30) shows that it suffices to take into 
account G^ to second order (G.5). Wc get 

<r''' \AGk\ <Pn'cs> = Afc 6l[, 51, + Tr {Fncs Kn'n) 
+ V Tr {Fncs Fi AA Kn'n) 

- E l^ ^^^2 ^ i^ncs AA Fi AA Kn'n) + (deg < 2). (G.15) 
In the case c ^ c', we can rewrite the expectation value as follows, 

<(l/"^ \AGk\ (t>n'cs> = <0"''^ I A Gk Kn'n i \ <t>ncs> = Tr {Fncs A Gk Kn'n i) ■ 

If we substitute in (G.l) and (G.5), factor the resulting matrix traces and use the 
estimates of the elementary traces of §G.3, we obtain that 

K^"' \AGk\ 4>n'cs> = (deg < 2) . (G.16) 

In order to bring the matrix elements (G.IO) into a suitable form, we substitute 
the definitions (G.6) and (G.8) into (G.IO) to obtain 

Aye's = E is-')irs<Gir'^\A\Gkci)n'c's> 

(ncs)G/ 

= E is-')'^rs<r'''\AGk\cPn'c's>, 

{ncs)£l 

where in the last step we used that Gk commutes with A (as the projector on an 

invariant subspace). Putting in the expansions (G.14) and (G.15, G.16) gives the 

result. ■ 

If there are no degeneracies, the above theorem reduces to the well-known formula 
of second order perturbation theory. The important result is that to the considered 
degree on the light cone, the matrix elements A^,^ are all zero if c 7^ c'. In other 
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words, the left- and right-handed components are invariant subspaces of A. This fact 
immediately gives the following corollary. 

Corollary G.4.2. Consider the fermionic projector in the presence of chiral and 

scalar potentials (6.0.5), were in composite expressions we disregard all terms quadratic 
in the field strength. Suppose that the matrix A^^,^, (G.ll), is diagonal in the sector 
indices n, n' for all k = 2,. . . ,K. Then for n = 1, . . . , 7, the contributions to the 
eigenvalues of degree two are 

8 1 

AA„L+ = TV(F„L+ AA) + V ^ TV(F„L+ F,,,r_ A^) (G.17) 

Kl+ - K'R- 

8 1 

AA„i_ = TV(F„i_ A^) + V T,{F^L- AA F^m^ AA) . (G.18) 

The traces appearing here are given explicitly by (G.31-G.48), where the line inte- 
grals are in phase-free form. The corresponding formulas for the opposite chirality are 
obtained by the replacements (G.22). 

Proof. The result is an immediate consequence of Theorem G.4.1 and the esti- 
mates (G.23-G.29). ■ 



G.5. Perturbation of the Kernel 

The results of the previous section do not apply to the kernel of A. The reason is 
that for A; = 1, the index set /, (G.l), is 

/ = {(ncs) with n = 8, c = L/i?, s = it} , 

and this index set contains both elements with s = + and s = — , giving rise to 
different types of matrix elements. On the other hand, the situation for the kernel is 
easier because the unperturbed spectral projector on the kernel satisfies the relations 

X*FiX = (G.l) 

XrFiX = = X*FiXL, (G.2) 

and furthermore we can simplify the calculations using that Ai = 0. Using these 
relations, it follows that, neglecting all contributions of degree < 2, the dimension of 
the kernel is not affected by the perturbation. 

Theorem G.5.1. Consider the fermionic projector in the presence of chiral and 
scalar potentials (6.0.5) and assume that the fermionic projector is weakly causality 
compatible (see Def. 7.1.1). Suppose that in composite expressions all terms quadratic 
in the field strength are discarded. Then 

AGx = (deg<2). 
Proof. Using the definition (G.2), 

AGi = — / Aiz- A)-^ dz = — I iz{z-A)~^-l)dz 
^ </ z{z- A)"^ dz . 



27ri 71^1 
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Performing the perturbation expansion gives, similar to (G.3), 



oo 

AGi = ^ — f z{iz- Ao)-' A^)" {z - Ao)-' dz . (G.3) 



n=0 "^Fl-e 

When we substitute in (G.4) and carry out the contour integral with residues, we get 
zero unless the factor z is differentiated. For this to occur, the pole at 2; = must 
be at least of order two, and thus wc need to take into account only the orders in 
perturbation theory n > 2. If n > 2, we can as in the previous section transform the 
matrix products into matrix traces, and the estimate (G.30) yields that the resulting 
contributions to AGi are of degree < 2. Thus it suffices to consider the second order 
in perturbation theory. 



^Gi = — 



J- f z(z- Ao)-^ AA (z - Ao)-^ AA (z - Aq)"^ dz + (deg < 2) 

2m J\z\=^ 



= - V ^ (i^z Ayl Fi AA Fi + Fi AA Fi AA Fi + Fi AA Fi AA Fi) 
^ A/ 

1=2 ' 

+ (deg < 2) (G.4) 

The weak causality compatibility condition implies that 

XP{x,y) = P{x,y) = P{x,y)X*, (G.5) 

and similarly for composite expressions in the fermionic projector. As a consequence, 
the first two matrix products in (G.4) vanish; namely, 

A A Fi AA = {AA X*) Fi {X AA) = A A {X* Fi X) AA ^ . 

In the last matrix product in (G.4) we can apply (G.2), 

Fi AA Fi AA Fi = Fi {X AA) Fi {AA X*) Fi = xl Fi AA Fi AA Fi xr ■ (G.6) 

Next we substitute (G.l), rewrite the resulting operator products as matrix traces, 
factor these matrix traces into elementary traces, and apply the estimates of §G.3. 
This straightforward calculation shows that the matrix product (G.6) is of degree < 5 
on the light cone. From (G.4) we conclude that AGi is of degree < 2. ■ 
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